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FOREWORD 


This  report  documents  the  effort  performed  under  contract  number  N 00024- 
80-C-6292  for  the  Naval  Sea  Systems  Command  Code  63R ,  and  covers  the 
period  from  October  1980  to  September  1981.-  The  study  has  considered  the 
use  of  LtVlS  adaptive  cancellation  to  suppress  convergence  zone  surface  rever¬ 
beration  in  a  bistatic  sonar.  The  work  was  directed  by  Mr.  Daniel  Porter 
of  NAVSEA  Code  63R . 


ABSTRACT 


•  -J 

This  report  considers  the  cancellation  of  surface  reverberation  from  the 
convergence  zone  (CZ)  using  an  L\1S  multiple  canceller  structure.  In  CZ 
operations,  this  surface  reverberation  is  often  the  limiting  factor  in  the  de¬ 
tection  of  low  doppler  targets.  The  multiple  reference  canceller  approach  is 
bistatic  in  that  hydrophones  spatially  separated  from  the  active  sonar  trans¬ 
mit/receive  array  are  used  as  references  to  cancel  the  CZ  surface  reverber¬ 
ation.  When  the  target  is  submerged  below  the  CZ,  providing  spatial  separa¬ 
tion  of  the  target  and  surface,  it  may  be  possible  to  spatially  reject  the  CZ 
surface  reverberation  but  not  the  target,  allowing  detection. 

The  first  phase  of  the  study,  reported  in  this  document,  concentrates 
on  the  ability  of  the  canceller  to  reject  the  reverberation  in  the  absense  of 
signal.  The  CZ  surface  reverberation  is  modeled  as  an  extended  source  with 
range  extent  limited  by  the  width  of  the  CZ  and  horizontal  extent  limited  by 
the  horizontal  directivity  of  the  transmitter.  When  viewed  from  typical  con¬ 
vergence  zone  ranges,  this  source  appears  to  have  narrow  angular  extent. 
Using  the  narrow  extended  source  model,  expressions  are  developed  for  the 
second  order  statistics  of  the  canceller  output.  It  is  shown  that  the  solution 

A  .  , 

to  the  spatial  cancellation  problem  with  a  narrow,  extended  source  is  func¬ 
tionally  equivalent  to  the  linear  prediction  of  a  bandlimited,  temporal  sequence. 
This  equivalence  not  only  allows  use  of  the  results  of  linear  prediction  theory 
in  the  spatial  case,  but  provides  an  intuitive  feel  for  the  parametric  behavior 
of  the  spatial  canceller. 

Using  the  analytical  results,  the  cancellation  achieved  is  examined  as  a 
function  of  the  reverberation  and  canceller  parameters.  It  is  shown  that  as 
the  number  of  references,  K,  is  increased,  the  cancellation  improves  rapidly 


r 


up  to  some  threshold  value,  say,  K-Ko,  the  improvements  diminish  rapidly. 
Given  the  computational  cost  of  adding  a  reference,  the  use  of  K>Ko  refer¬ 
ences  may  not  be  justified.  The  value  of  Ko  is  determined  approximately, 
and  provides  a  design  guide  for  the  selection  of  the  number  of  references 
in  any  given  sonar  situation.  Significantly,  effective  cancellation  of  the 
reverberation  can  be  achieved  using  less  than  four  references.  Extensive 
curves  showing  the  parametric  dependencies  of  the  cancellation  on  the  de¬ 
sign  parameters  are  given  as  design  tools.  The  report  also  describes 
several  configurations  of  the  LMS  canceller  structure  for  use  in  the  non¬ 
stationary  reverberation  environment  and  discusses  their  performance. 

The  second  phase  of  this  study  will  focus  on  two  areas  r  the  incorpora¬ 
tion  of  a  plane  wave  signal  in  the  model;  and  the  verification  of  the  extended 
source  model  as  representative  of  the  CZ  surface  reverberation.  This  report 
on  the  first  phase  of  the  study  includes  some  initial  determination  of  the 
response  of  the  canceller  to  the  signal.  These  results  must  be  extended  to 
include  the  effect  of  the  primary  array  and  the  detailed  characteristics  due 
to  separation  of  the  target  and  the  reverberating  surface.  The  performance 
of  the  canceller  will  ultimately  be  evaluated  in  terms  of  the  improvement  in 
detection  performance  provided  by  the  canceller  over  that  achieved  without 
the  canceller.  Again,  this  will  be  examined  as  a  function  of  the  canceller 
design  parameters,  reverberation  extent,  and  the  spatial  separation  between 
surface  and  target. 

To  assess  the  validity  of  the  analytical  models,  a  computer  simulation 


of  the  multiple  canceller  structure  will  be  developed.  This  will  first  be 
tested  against  computer  generated  data  simulating  the  active  sonar  in  the 
CZ  detection  mode.  More  importantly,  the  simulation  will  be  run  on  actual 


data  from  sea  tapes  to  be  provided  by  the  Navy.  Any  changes  to  the  models 
necessitated  by  the  simulation  runs  will  be  made  in  the  analysis  and  the 
effect  of  these  changes  on  canceller  performance  examined. 


I 
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1.0  INTRODUCTION 


This  is  the  final  report  on  the  first  phase  of  a  study  of  a  technique  for 
processing  against  surface  reverberation  from  the  convergence  zone  in  an 
active  sonar.  When  attempting  to  detect  a  low  doppler,  submerged  conver¬ 
gence  zone  (CZ)  target,  this  backscatter  of  the  transmitted  waveform  from 
the  surface  above  the  target  can  be  the  limiting  background  noise.  The 
concept  involves  bistatic  or  multistatic  operation,  in  which  reference 
hydrophones  spatially  separated  from  the  primary  transmit /receive  array 
of  the  sonar  are  used  to  suppress  the  reverberation  in  the  primary  array 
output.  This  is  done  using  a  Least  Mean  Square  (LMS)  adaptive  multiple 
canceller  structure,  which  has  been  used  in  the  past  to  suppress  point 
interference  in  the  output  of  a  beamformer^  .  The  CZ  surface  reverber¬ 
ation,  however,  appears  as  a  narrow  extended  acoustic  interference,  so 
this  study  addresses  the  performance  of  the  multiple  canceller  concept 
with  extended  sources. 

The  most  common  (and  most  practical)  algorithm  for  the  implementation 
of  spatial  noise  cancellation  is  the  Least  Mean  Square  (LMS)  algorithm  de¬ 
scribed  in  [2] .  Define  the  inputs  to  the  algorithms  as 

d(n)  =  desired  input  at  n  a  sample  time 

x(n)  =  reference  input  at  n  n  sample  time. 

The  LMS  algorithm  stores  the  M- dimensional  data  vector 

X(n)  =  [x(n),  x(n-l) ,  - x(n-M+l)]T  (1-1) 

% 

and  computes  the  filter  output, 

y(n)  =  WT  (n)  X  (n)  ( 1-2) 

with 


1 


W(n)  =  [w0(n),  W^n),  •■•,WM.1(n)J 


the  weight  vector  of  length  M.  The  weights  are  updated  recursively  to 
minimize  the  mean  squared  value  of  the  error  signal,  defined  as 

e(n)  =  d(n)  -  y(n)  ( 1-3) 

This  is  done  by  means  of  a  gradient  descent  approach  in  which  the  actual 
gradient  of  the  mean  squared  error  is  replaced  by  an  estimate  of  the  gradi¬ 
ent  extracted  from  the  data  vector,  X(n).  This  yields  the  very  simple 
weight  update  equation 

W(n+1)  =  W(n)  +  fi€  (n)  X(n)  (1-4) 

where  m-  is  a  weight  update  coefficient  that  controls  the  rate  of  convergence 
and  stability  of  the  algorithm. 

When  the  inputs  to  the  LMS  adaptive  structure  are  zero  mean  random 

sequences  that  are  at  least  wide  sense  stationary,  the  weights  given  by 

(1-1)  can  be  shown  to  converge  in  the  mean  to  the  discrete  Wiener  filter 

r  31 

provided  h  is  sufficiently  small  to  assure  stability  .  That  is,  if  E[x]  de¬ 
notes  the  expected  value  of  x  and 


E[d(n)]  =  E[x(n)]  =  0 


xx 


=  E[X(n)  X  (n)] 


k  =  0,  1,  ...  K-l 


(1-5) 


and 


^-dx 


=  E[d(n)  X(n)] 


then 


lim  E[W(n)]  =  fT1  r^ 
n—  ® 


(1-6) 
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under  the  appropriate  conditions  on  y.  Reference  [3]  shows  that  (1-6) 
occurs  if 


,<l‘<  l'W  <1'7> 

where  ^max  is  the  largest  eigenvalue  of  the  covariance  matrix,  Rxx- 
Further,  it  is  shown  that  the  mean  of  the  weights  of  the  adaptive  filter 
approach  the  Wiener  solution,  (1-6),  approximately  exponentially  with  the 
time  constant 


TI 


„_J _ 

^\nin 


(1-8) 


where  r,  is  expressed  in  iterations  of  the  filter  and  X  .  is  the  smallest 
i  mm 

eigenvalue  of  the  covariance  matrix,  Rxx  • 


Because  the  weights  are  computed  from  the  random  input  data  in  (1-4), 
they  are  themselves  random,  and  exhibit  fluctuations  about  the  discrete 
Wiener  solution,  (1-6).  Therefore,  the  power  in  the  error,  e(n),  is  greater 
than  the  error  that  would  be  achieved  by  the  Wiener  weights.  This  increase 
in  error  power  due  to  algorithm  noise  is  termed  the  misadjustment  [3]  of 
the  algorithm.  In  general,  due  to  the  recursive  nature  of  (1-4),  even  the 
second  order  statistics  of  the  weights  cannot  be  determined  analytically  in 
closed  form.  As  will  be  discussed  below,  one  way  of  approximating  the 
weight  variance  is  to  develop  a  frequency  domain  model  for  the  adaptive 
filter,  as  in  [4],  and  to  then  use  [5]. 


Now,  consider  using  this  LMS  adaptive  filter  to  adaptively  cancel  a 
single  plane  wave  interference  from  the  output  of  a  beam  steered  in  the  di¬ 
rection  of  a  signal  of  interest.  A  single  hydrophone  is  chosen  as  a  refer¬ 
ence,  providing  the  reference  input  to  the  LMS  algorithm,  x(n).  The  de¬ 
sired  input,  d(n),  is  the  output  of  the  beam  steered  toward  the  signal 


which  is  assumed  to  arrive  from  a  different  angle  than  the  interference. 
Because  the  interference  dominates  both  the  reference  hydrophone  and 
beam  outputs,  the  adaptive  filter  minimizes  the  mean  squared  value,  or 
power,  in  e(n)  primarily  by  eliminating  the  interference  component.  It 
does  this  by  attempting  to  insert  a  delay  in  the  reference  input  equal  to 
the  propagation  delay  between  the  phase  center  of  the  primary  array  and  the 
reference  hydrophone,  plus  a  gain  equal  to  the  primary  array  gain.  The  in¬ 
terference  component  of  the  filter  output  is  then  nearly  coherent  with  the 
primary  array  output,  so  that  the  interference  is  cancelled  in  the  error  wave¬ 
form,  which  serves  as  the  canceller  output.  In  order  to  cancel  the  signal,  the 
filter  would  have  to  adapt  to  the  propagation  delay  and  array  gain  for  the 
signal,  which  differ  from  those  for  the  interferences  because  the  arrival 
angles  differ.  In  doing  this,  the  structure  would  pass  the  much  stronger 
interference,  increasing  the  mean  squared  error.  The  adaptive  algorithm 
therefore  converges  to  a  solution  which  cancels  the  interference  but  passes 
the  signal.  In  a  spatial  sense,  this  can  be  viewed  as  steering  a  null  in  the 
beam  response  in  the  direction  of  the  plane  wave  interference.  Generally,  the 
signal  will  be  somewhat  attenuated  by  the  cancellation  process,  but  this  will 
be  tolerable  in  view  of  the  much  greater  rejection  of  the  interference. 

If  the  interference  must  be  modelled  as  more  than  one,  say  K,  plane 
wave  sources,  then,  in  general,  K  spatially  separated  references  will  be 
needed  to  allow  spatial  nulls  to  be  steered  in  the  direction  of  each.  The  LMS 
algorithm  can  be  configured  to  use  K  reference  inputs  as  follows.  Let 


xk(n)  =  kth  reference  input  at  n**1  sample  time  (k  =  0,  1,  ...  K-l) 


and  let  the  algorithm  store  the  K  data  vectors  of  length  M , 

X.  (n)  =  [x.  (n) ,  x.  (n-1),  ...  x.  (n-M+l)]T,  k  =  0,  1, 


K-l 

(1-9) 


The  LMS  algorithm  computes 

yk(n)  =  wj(n)  Xk(n),  k  =  0,  1,  ...  K-l  (1-10) 

where  Wk(n)  is  the  k**1  weight  vector,  with  the  error  given  by 

K-l 

e(n)  =  d(n)  -  £  yk(n),  k  =  0,  1.  ...  K-l  (1-11) 

k  =  0 

The  weights  are  updated  according  to 

Wk(n+1)  =  Wk(n)  +  e(n)  Xk(n),  k  =  0,  1,  ...  K-l  (1-12) 

Note  that  each  implementation  of  (1-10)  is  a  non-recursive  digital  filter 
with  time  varying  weights,  so  the  multiple  reference  LMS  algorithm  can  be 
represented  as  shown  in  Figure  1-1.  Each  reference  input  is  passed  through 
a  non-recursive  digital  filter,  the  filter  outputs  used  to  form  the  error.  c(n) 
by  (1-13),  then  the  error  fed  back  to  update  the  filter  weights  according  to 
(1-4). 


In  many  practiced  situations,  the  computational  requirements  of  implement¬ 
ing  (1-16),  (1-17),  and  (1-18)  are  prohibitive,  and  the  cancellation  must  be 
done  using  the  so-called  frequency  domain  LMS  algorithm  [4,  6],  The  fre¬ 
quency  domain  algorithm  first  Fast  Fourier  Transforms  (FFT's)  the  reference 
and  primary  inputs,  forming 
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REFERENCE 
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Figure  1-1.  Multiple  Reference  LMS  Algorithm 


D.(n) 


=  i**1  FFT  coefficient  at  time  n 


M-l  ,  2n 

^  d(m  +  nM)  e 
m=0 


]  77-  mi 
1  M 


(1-13) 


and 


X^(n)  =  ^  FF^  coefficient  of  k**1  reference  at  time  n 


M-l  .  2t t  . 

V  '1  \T  mi 

=  2  xk(m  +  nM)  e  (i  =  0,  1,  ...M-l)  (1-14) 


m=0 


where  M  is  again  the  dimension  of  the  X^(n)  data  vectors.  A  single  tap  com¬ 
plex  LMS  algorithm  is  then  performed  independently  in  each  frequency  bin. 
Let 


F^fnl  =  the  frequency  domain  weight  in  the  l  1  FFT  bin  for  the 
til 

k  reference  at  time  n 

Then  the  LMS  filter  output  in  the  i**1  bin  for  the  kth  reference  at  time  n  is 

Yik(n)  =  Fik(n)  Xik(n)  (1'15) 

and  the  error  in  the  i**1  FFT  bin  (used  to  update  all  K  filters)  is 

K  - 1 

E.(n)  =  D.(n)  -  ^  (n)  (1-16) 

k=0 

The  complex  LMS  algorithm  [7]  gives  the  updated  weights  as 

Fik(n+1)  =  Fik(n)  +  Ei(n)  Xik(n)  (1-17) 


7 


If  required,  the  time  domain  canceller  output  can  be  obtained  by  inverse 
FFTing  the  E.(n) . 

While  this  algorithm  gives  significant  computational  savings  over  the 
time  domain  approach,  it  is  not  an  exact  implementation  of  the  time  domain 
algorithm  and  can  give  different  results,  even  in  relatively  simple  cases  [4]  . 
A  number  of  practical  cancellation  problems  have  been  solved  using  the  fre¬ 
quency  domain  approach,  however  [8,  9].  When  it  can  be  used,  it  has  some 
analytical  advantages  in  that  it  may  be  possible  to  determine  the  variance 
of  the  frequency  domain  weights  as  well  as  the  mean  using  [8]  .  This  allows 
the  effects  of  algorithm  noise  to  be  considered  analytically. 

The  CZ  surface  reverberation  differs  from  interferences  usually  associ¬ 
ated  with  LMS  spatial  cancellation  in  three  ways: 

( 1)  Extent ;  The  first  CZ  will  nominally  appear  at  a  range  of  30  miles 
and  have  a  width  of  approximately  2  miles,  so  that  it  appears  to 
have  a  narrow  vertical  angular  extent.  In  addition,  most  modern 
active  sonars  have  horizontal  directivity,  so  that  only  a  hori¬ 
zontal  sector  of  the  CZ  is  illuminated.  The  reverberating  CZ  sur¬ 
face  therefore  appears  as  an  extended  source  in  both  horizontal 
and  vertical  angles.  Within  this  sector,  the  reverberation  is 
made  up  of  a  multitude  of  individual  reflections,  often  regarded 
as  uncorrelated  [  10] .  The  surface  may  therefore  be  regarded  as 
narrow  extended  source  exhibiting  no  correlation  from  angle  to 
angle.  This  study  therefore  addresses  the  performance  of  the 
canceller  with  such  a  source. 
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(2)  Non -stationarity ;  In  the  active  sonar  environment,  the  input 
statistics  to  the  canceller  are  markedly  non -stationary ,  consist¬ 
ing  of  a  transmit  period  followed  by  periods  of  reverberation 
from  various  sources,  intervals  of  ambient  noise  only,  and  signal 
returns.  Although  the  derivation  of  the  LMS  adaptive  canceller 
assumes  stationarity  of  the  inputs,  it  has  often  been  applied  in 
non -stationary  situations  by  making  its  convergence  time  short 
enough  to  track  the  non -stationarities .  An  approach  to  the  ap¬ 
plication  of  the  LMS  algorithm  in  the  reverberation  environment 
is  also  developed  in  this  study  and  its  performance  analyzed. 

(3)  Noise  Related  to  Signal:  In  most  applications  of  the  adaptive 
canceller,  the  signal  and  noise  are  derived  from  different 
sources,  and  are  therefore  uncorrelated.  In  the  detection  of  a 
target  against  the  background  of  CZ  surface  reverberation,  the 
noise  consists  of  many  reflections  of  the  transmitted  waveform. 
Since  reflections  of  the  waveform  also  make  up  the  signal,  the 
signal  and  noise  are  related.  This  will  not  be  of  central  interest 
here,  since  this  report  concentrates  on  the  signal  absent  case. 

Suppose  that  instead  of  consisting  of  one  or  more  discrete  plane  wave 
sources,  the  interference  is  truly  extended  in  that  it  is  distributed  in 
angle  over  some  finite  sector.  Given  the  argument  that  the  adaptive  can¬ 
celler  is  steering  nulls  in  the  overall  response  pattern  in  the  direction  of  the 
interference,  it  would  seem  that  an  extended  source  could  never  be  completely 
cancelled.  However,  as  more  and  more  references  are  added  (K  increased), 
steering  an  increasing  number  of  nulls  to  the  interference,  it  seems  clear 
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that  the  amount  of  additional  cancellation  achieved  per  added  reference 
would  eventually  become  quite  small.  Since  the  computational  cost  of  adding 
a  reference  is  significant  (requiring  an  additional  LMS  adaptive  filter),  at 
some  point  the  improvement  will  not  be  worth  the  cost.  This  study  examines 
the  performance  of  the  multiple  canceller  approach  with  extended  sources 
using  a  Wiener  filter  model.  The  amount  of  cancellation  achieved  is  examined 
as  a  function  of  the  source  extent,  the  number  of  references,  and  their 
location  with  respect  to  the  primary  array.  Of  particular  interest  is  the 
following 

a.  The  degree  of  cancellation  theoretically  achieveable  for  an 
extended  source 

b.  The  selection  of  the  number  of  inferences  for  a  given  source  ex¬ 
tent,  taking  into  account  both  cancellation  performance  and 
computational  cost. 

The  Final  Report  is  divided  into  three  sections,  the  first  of  which  is  this 
Introduction.  Section  2  summarizes  the  results  of  the  study,  describes  the 
general  approach  used,  and  discusses  the  models  of  the  adaptive  canceller 
and  extended  interference  source.  The  detailed  development  of  the  models 
and  the  derivations  of  the  results  given  in  Section  2  are  included  as  Appen¬ 
dices  A  through  K.  Section  3  presents  the  general  conclusions  of  the  study 
as  to  the  applicability  of  the  LMS  cancellation  technique  to  the  rejection  of 
surface  reverberation  from  the  convergence  zone.  It  also  discusses  areas 
where  additional  study  is  needed . 
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2.0  SUMMARY  OF  RESULTS 
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2.0  SUMMARY  OF  RESULTS 


This  section  presents  the  results  of  the  study  in  summary  form  without 
detailed  derivations,  which  are  included  in  the  Appendices.  Where  results 
are  best  illustrated  by  the  use  of  graphs,  typical  cases  are  shown  here  and 
the  complete  set  of  curves  given  in  the  Appendices.  Section  2.1  discusses 
the  application  of  the  LMS  algorithm  in  the  non -stationary  reverberation 
environment.  Section  2.2  develops  the  model  for  the  extended  interference 
source  and  relates  it  to  surface  reverberation  from  the  convergence  zone. 

In  Section  2.3,  the  frequency  domain  Wiener  filter  model  for  the  adaptive 
canceller  is  derived.  Section  2.4  shows  that  there  exists  an  equivalence  be¬ 
tween  the  linear  prediction  of  a  bandlimited  temporal  random  process  and  the 
cancellation  of  an  extended  source.  This  equivalence  is  used  to  apply  a 
number  of  results  from  linear  prediction  theory  to  the  problem  of  interest 
here.  The  general  results  of  Section  2.4  are  then  specialized  to  a  source 
that  is  uniformly  distributed  in  angle  and  spatially  uncorrelated  in  Section 
2.5.  This  model  is  shown  to  be  a  reasonable  model  for  CZ  surface  reverbera¬ 
tion.  Performance  predictions  for  the  multiple  canceller  in  this  enviornment 
are  presented,  and  design  criteria  given  for  selection  of  hydrophone  spac¬ 
ing  and  number  of  references.  Section  2.6  discusses  the  performance  of  the 
canceller  with  a  source  producing  an  exponential  CSD  at  the  hydrophone 
outputs,  and  compares  this  performance  to  that  achieved  with  the  uniform 
source  of  Section  2.5.  Section  2.6  investigates  at  the  effects  of  the  non- 
stationarity  of  the  reverberation  on  the  performance  of  the  adaptive  cancel¬ 
ler.  Finally,  Section  2.7  takes  a  preliminary  look  at  the  response  of  the 
array  to  a  plane  wave  signal  in  the  presence  of  an  extended  source. 
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2 . 1  Canceller  Structure  for  Application  in  the  Non -Stationary  Reverber¬ 
ation  Environment.  As  pointed  out  in  the  Introduction,  the  development  of 
the  LMS  canceller  structure  is  based  upon  the  stationarity  of  the  inputs, 
although  the  adaptive  algorithm  itself  has  been  applied  to  non -stationary 
problems  [9,  11]  .  This  is  done  by  making  the  time  constant  of  the  adaptation 
process  short  enough  that  the  weights  can  "track"  the  non-stationarities. 

This  requires  use  of  value  of  n  larger  than  would  be  needed  in  a  stationary 
environment  with  a  resulting  increase  in  algorithm  noise.  There  is  there¬ 
fore  a  trade-off  between  the  ability  of  the  algorithm  to  follow  input  dynamics 
and  the  algorithm  noise. 

The  application  of  the  LMS  canceller  approach  being  considered  here 
departs  from  this  usual  approach  to  non -stationary  inputs.  The  active  sonar 
environment,  taken  as  a  whole,  is  markedly  non -stationary ,  consisting  of 
periods  during  which  reverberation  from  various  sources  (surface,  bottom, 
convergence  zone,  etc)  dominates  the  noise  field.  Between  these  periods  of 
reverberation  the  ambient  noise  field  will  be  the  main  limitation  on  sonar  per¬ 
formance.  However,  within  a  period  of  reverberation  from  a  particular  source, 
the  reverberation  is  often  regarded  as  a  stationary  random  process.  Hence, 
it  may  be  possible  to  avoid  the  need  for  increasing  y  to  track  input  dynamics 
by  adapting  the  weights  only  when  the  convergence  zone  (CZ)  reverberation 
is  present.  During  these  periods,  the  reverberation  dominates  the  ambient 
noise  component  and  the  signal,  so  it  should  be  possible  to  detect  its  presence 
reliability . 

These  considerations  suggest  the  following  canceller  structure.  Because 
the  propagation  distance  to  the  CZ  is  so  much  greater  than  that  associated 
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with  other  forms  of  reverberation,  it  is  possible  to  exclude  all  other  rever¬ 
beration  from  influencing  the  canceller  by  freezing  the  weights  for  Tw  sec¬ 
onds  after  each  transmission.  After  seconds,  it  can  be  assumed  that  the 
input  to  the  canceller  is  either  ambient  noise  or  CZ  surface  reverberation 
(possibly  with  a  target  present).  The  presence  of  the  reverberation  is  then 
determined  by  a  simple  energy  detection  in  the  reverberation  band  of  the 
primary  array  output,  and  the  weights  adapted  only  during  its  presence. 

This  structure  is  shown  schematically  in  Figure  2-1.  If  the  detector  can  de¬ 
termine  the  presence  of  the  CZ  reverberation  perfectly,  the  weights  only 
adapt  in  the  presence  of  the  reverberation.  Given  that  this  reverberation  is 
stationary,  as  assumed  above,  the  inputs  as  seen  by  the  weight  adaptation 
process  are  stationary  at  all  times,  so  that  they  converge  to  the  optimal 
weights  that  would  occur  if  the  reverberation  were  always  present.  The 
selection  of  u  can  be  made  solely  on  the  basis  of  acceptable  algorithms  mis- 
adjustment.  False  alarms  in  the  detector  will  cause  the  weights  to  adapt 
in  an  ambient  noise  background,  which  will  degrade  the  effectiveness  of  the 
canceller.  False  dismissals  in  the  detector  will  result  in  the  weights  being 
frozen  during  the  presence  of  the  CZ  reverberation,  so  that  the  algorithm 
will  take  longer  to  change.  It  is  anticipated,  though,  that  the  high  rever¬ 
beration  to  ambient  noise  ratio  associated  with  this  problem  will  allow  detector 
operation  at  a  low  probability  of  error. 

An  alternative  approach  can  be  based  on  the  well  known  fact  [  10]  that 
the  LMS  adaptive  algorithm  adapts  much  more  rapidly  in  an  environment 
producing  correlation  between  its  inputs  (such  as  reverberation  from  any 
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Figure  2-1.  Multiple  Canceller  Structure  for  Application 
to  CZ  Surface  Reverberation 


source)  than  to  one  producing  uncorrelated  inputs  (such  as  ambient  noise 
field  would).  That  is,  the  filter  "learns"  the  correlation  properties  of  an  input 
rapidly  in  comparison  to  the  way  it  "forgets"  those  properties  once  the  cor¬ 
relation  disappears.  Therefore,  in  order  to  maintain  the  properties  of  the 
convergence  zone  reverberation  in  the  filter  weights,  it  is  first  essential 
that  all  other  forms  of  reverberation  be  excluded  from  the  adaptation  proc¬ 
ess  by  inhibiting  the  adaptation  of  the  weights  for  T^  seconds  after  each 
transmission.  From  the  point  of  view  of  the  adaptive  weights,  then,  the 
input  consists  of  the  intervals  fnT  +  T  ,  (n+l)T]  for  n  =  1,  2,  ...  con¬ 
catenated  together  to  form  a  continuous  input  sequence.  This  concatenated 
input  is  then  sampled  every  T  seconds  to  produce  the  inputs  x(n)  and  d(n). 
The  question  is  how  much  the  adaptation  in  ambient  noise  only,  albeit  slow, 
degrades  the  operation  of  the  canceller  when  the  reverberation  actually  com¬ 
mences.  This  is  considered  in  Section  2.7. 

2.2  Model  for  the  Acoustic  Field.  The  majority  of  the  results  of  this 
study  are  based  upon  a  far  field  extended  source  model  developed  in 
Appendix  A  for  the  particular  array  geometry  shown  in  Figure  2-2.  In  this 
geometry,  the  K  reference  hydrophones  are  configured  in  a  line  array  with 
an  inter-hydrophone  spacing  of  d  feet.  These  K  references  are  to  be  used 
in  an  LMS  canceller  to  suppress  an  extended  interference  from  the  output  of 
a  single  omnidirectional  hydrophone  located  md  feet  from  the  nearest  reference 
and  colinear  with  the  reference  line  array.  In  this  model,  the  transmit/ 
receive  array  of  the  sonar  has  been  replaced  by  a  single  omnidirectional  hy¬ 
drophone.  The  impact  of  this  simplification  is  discussed  later. 
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The  extended  interference  is  assumed  to  be  distributed  in  a  horizontal 
plane  that  contains  the  line  array,  as  shown  in  Figure  2-2.  The  acoustic 
energy  arriving  at  the  array  from  the  interference  at  any  angle,  <J>  ,  is 
assumed  to  produce  a  plane  wave  across  the  entire  array  (including  refer¬ 
ences  and  primary).  This  requires  that  the  distance  to  any  point  on  the 
interference  be  large  in  comparison  to  the  array  dimensions  so  that  the 
wavefront  curvature  is  negligible.  As  is  usually  done  with  plane  wave 
models,  it  is  also  assumed  that  the  propagation  losses  from  a  point  on  the 
source  to  any  hydrophone  is  the  same,  so  that  the  difference  in  propagation 
to  different  hydrophones  is  characterized  in  terms  of  delay  only.  The  com¬ 
mon  propagation  losses  can  be  lumped  with  the  source  characteristics. 

Statistically,  it  is  assumed  that  the  acoustic  energy  arriving  at  each 
angle.  4>  ,  from  the  interference  is  a  zero  mean,  stationary  random  process. 

In  the  general  model,  the  arrivals  from  two  angles,  <J>^,  and  are  allowed  to 
be  correlated,  producing  a  cross -spectral  density  (single  frequency  correla¬ 
tion)  between  the  arrivals  of  Sj(<J>j,  <f>  ^ ,  w) .  Thus  the  spectral  density  of  the 
intererence  at  any  arrival  angle,  4>  ,  is  Sj(<t> ,  <t>  ,  w) .  In  addition  to  the  inter¬ 
ference,  the  acoustic  field  is  assumed  to  have  an  isotropic  background,  pro¬ 
ducing  zero  mean  noise  at  the  output  of  each  hydrophone.  This  noise  is 

assumed  to  be  uncorrelated  from  hydrophone  to  hydrophone  and  to  have 
2 

spectral  density  <rn(u)  • 

Next,  in  Appendix  A.l  the  model  is  restricted  to  narrow  extended  inter¬ 
ferences  that  are  distributed  over  a  narrow  angular  sector  fag  -#,  +  $] . 

This  allows  an  approximation  leading  to  a  two-dimensional  Fourier  transform 
relationship  between  the  source  angular  distribution  and  the  cross -spectral 
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density  (CSD)  of  the  hydrophone  outputs.  Appendix  A.  2  adds  the  additional 
assumption  of  a  spatially  uncorrelated  interference,  for  which 


Sj(<j>r  4>2,  W)  =  SjC^j.w)  6  («*> x  -  <t>2)  (2-1) 

In  this  case,  the  CSD  between  the  outputs  of  the  p**1  and  qth  reference 
hydrophones  is  given  by 


Although  it  is  not  immediately  apparent  that  this  simple,  far  field  model 
can  be  used  to  model  surface  reverberation  from  the  convergence  zone. 
Appendix  B  considers  two  source  geometries  representative  of  this  reverbera¬ 
tion  and  shows  that  they  reduce  to  the  form  of  the  far  field  model  under  very 
reasonable  assumptions.  Suppose  that  the  sonar  transmitter  has  a  horizontal 
beamwidth  of  29j  degrees,  so  that  it  only  illuminates  a  sector  of  the  annulus 
comprising  the  convergence  zone,  which  has  inside  radius  Rj  and  outside 
radius  Rq.  The  illuminated  sector  is  assumed  to  be  centered  about  the  angle 
eQ,  and  the  transmitter /receiver  is  located  D  feet  below  the  surface,  as 
shown  in  Figure  2-3.  The  reference  hydrophones  are  placed  d  feet  apart 
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reference  hydrophone 

hydrophone 


Figure  2-3.  Geometry  for  Surface  Model  with  Straight  Line  Propagation 


along  the  negative  y  axis,  with  the  distance  between  the  primary  and  the 
closest  reference  md  feet.  The  reverberation  from  each  point  on  the  sector 
of  the  annulus  is  assumed  to  propagate  via  a  straight  line  (direct)  path  to 

the  hydrophones,  and  can  be  regarded  as  a  plane  wave  across  the  array. 

2 

The  ambient  noise  is  assumed  to  produce  noise  with  spectral  density  crn(uO 
which  is  uncorrelated  from  hydrophone  to  hydrophone.  This  can  be  viewed 
as  a  model  of  the  convergence  zone  surface  reverberation,  except  that  the 
propagation  paths  are  direct  rather  than  along  the  refracted  paths  associated 
with  the  convergence  zone.  It  is  shown  in  Appendix  B  that  when  the  rever¬ 
beration  is  spatially  uncorrelated,  the  shape  of  the  propagation  paths  does 
not  affect  the  hydrophone  CSD  except  in  terms  of  the  overall  interference 
strength  (this  is  incorporated  in  the  interference  spectrum).  The  assumption 
that  the  reverberation  is  uncorrelated  from  point  to  point  is  commonly  made 
in  active  sonar  (see,  for  example,  Van  Trees  [9])  and  will  be  made  here. 

Now.  let  Sj(p,  0 ,  w)  be  the  spectral  density  of  the  reverberation  from 
the  point  on  the  surface  at  range  P  and  angle  0,  with  the  propagation  losses 
incorporated  in  the  source.  It  is  assumed  that  this  spectral  density  is 
separable  in  range, 

Sr(p,  6,  ui)  =  SIp(P,w)  SIa(0,  w)  (2-3) 

Then  under  the  additional  assumptions  that  0,  is  small  enough  that 
cos  e ^  =  1 

sin  0][  »  0X  (2-4) 

and  that  D  «P,  the  CSD  between  the  pt^  and  qth  hydrophones  in  the  array 
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Sir(P.^)  P “  dp 
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Comparing  (2-5)  with  (2-2)  shows  that  the  model  illustrated  in  Figure  2-3 
produces  the  same  form  for  the  hydrophone  output  CSD  as  the  far  field 
source  with  Sj(<t>+<^,cj)  replaced  by  the  angular  source  spectral  density 
Sia(0+0o,w)and  with  the  additional  factor  G  (w).  This  additonal  factor  can 
just  be  regarded  as  part  of  the  source  spectral  density  in  the  far  field 
model  of  Appendix  A.  Hence,  the  far  field  model  developed  in  Appendix  A 
can  be  applied  to  the  CZ  surface  reverberation  problem. 

Appendix  B  also  develops  a  second  model  for  the  CZ  reverberation  that 
utilizes  the  assumption  of  plane  wave  arrivals  at  the  array  from  each  point 
on  the  source,  but  distributes  these  arrivals  on  the  horizontal  sector,  $ 
e  ®  +  9  il  ’  and  the  vertical  sector, -$  .  0  Q  +*]  as  shown  in 

Figure  2-4.  When  the  reflections  from  each  point  on  the  surface  are  uncor¬ 
related,  the  arrivals  from  different  angles  will  also  be  uncorrelated.  As  be¬ 
fore,  the  horizontal  extent  is  limited  by  the  horizontal  directivity  of  the 
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transmitter.  The  vertical  sector  comprises  those  angles  producing  CZ  propa¬ 


gation.  Let  Sj(0,o.~)  be  the  spectral  density  of  the  arrival  at  horizontal 
angle,  e  ,  and  vertical  angle,  <t>,  and  assume  that 

Sj(0. *>.«-)  =  Sjh(0,  <*0  Sjv(4>,  )  (2-7) 

That  is,  S j (9  4» ,  u  )  is  separable  into  a  horizontal  density,  Sj^(6,u),  and  a 
vertical  density,  Sjv(<t>,u,).  It  is  assumed  that  is  small  in  the  sense  of 
(  2-  4) ,  and  that 

max  [<t>0  0Q  +*]  <<  1  (2-8) 


so  that  cosImaxf'pQ  Q ^  +  $]]  »  1.  This  requires  that  the  vertical  arrivals 

be  clustered  near  horizontal,  and  is  consistent  with  CZ  arrivals  in  a  surface 
ship  sonar.  Then  the  CSD  of  the  outputs  of  the  pt}l  and  q**1  hydrophone  is 
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Again,  comparison  of  (2-9)  with  (2-4)  reveals  that  the  model  of  Figure  2-6 
produces  the  same  hydrophone  CSD  as  the  far  field  model  of  Appendix  A, 
so  that  the  latter  can  be  used  to  study  the  CZ  surface  reverberation  problem 
An  important  special  case  of  the  far  field  model  is  that  of  uniform  source 
spectral  density 


Sj($,w)  = 
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The  hydrophone  output  CSD  then  becomes 


sin  (p-q)*sin0Q 
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Note  that  for  the  surface  model  (Figure  2-3)  and  the  vertical  arrival  model 
(Figure  2-4),  (2-10)  is  replaced  by 


OTj  (^),  0e[0Q  -  »v  0Q  + 9 


SIa(0,u»  =  SIv(fl,u >)  = 


0,  otherwise 

This  requires  that  the  reverberation  intensity  does  not  vary  as  a  function 
of  bearing  within  the  transmitter  main  lobe. 

It  is  important  to  note  the  behavior  of  the  hydrophone  output  CSD, 
s  (p.q.w),  as  the  extended  source  approaches  end  fire.pg  =  0.  Whenp^  =  0, 
(2-2)  reduces  to 


se(p,q,w;  =^n(^)5pq  + 
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which  is  exactly  the  CSD  produced  by  a  single  end  fire  plane  wave  with 

i 

spectral  density 
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+  0 ,  d0 


Similarly,  when  0q  =  0,  (2-11)  becomes 
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(2-13) 
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se(p.q.j)  =  ‘’’n^pq  +  2#<Tl2(a,)  e 
which  is  again  the  hydrophone  CSD  of  an  end  fire  plane  wave.  Further, 

(2-11)  indicates  that  as  the  source  approaches  endfire,  its  hydrophone  out¬ 
put  CSD  broadens  in  terms  of  the  difference  (p-q),  so  that  it  looks  like  the 
CSD  is  that  of  a  narrower  source.  It  is  as  if  the  effective  extent  of  the 
source  is  2$sin0g  rather  than  the  actual  source  extent  of  2$.  It  is  important 
to  show  that  this  apparent  narrowing  of  the  extended  source  near  end  fire 
is  a  characteristic  of  the  physical  problem  of  cancelling  an  extended  source 
using  a  line  array  of  references,  and  not  just  an  anomaly  introduced  by  the 
model  or  by  approximations  within  the  model. 

Although  analysis  of  the  general  case  has  not  been  completed,  a  number 

of  numerical  evaluations  of  the  exact  hydrophone  output  CSD  have  been 

made  for  the  special  case  of  a  uniform  source.  These  agree  very  well  with 

the  CSD  derived  using  the  narrow  source  approximation,  given  by  (2-11)  in 
2 

the  case  where  on  (w)  =  0.  This  tends  to  support  the  validity  of  the  model 
near  end-fire  in  the  uniform  source  case,  and  suggest  that  the  apparent 
source  comparession  actually  occurs. 

In  most  cases  of  interest,  the  primary  sensor  will  not  be  a  single  omnidi¬ 
rectional  hydrophone  but  an  array  of  hydrophones  steered  in  some  particular 
direction.  Appendix  I  investigates  the  applicability  of  the  single  hydrophone 
model  to  this  more  realistic  case.  It  considers  a  primary  input  derived  from 
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a  line  array  of  N  hydrophones  spaced  every  d^  feet  and  located  a  distance 
Ld  feet  from  the  reference  array  as  shown  in  Figure  2-5.  The  results  show 
that  the  equivalent  model  of  Figure  2-6,  in  which  the  primary  array  is  re¬ 
placed  by  a  single  directional  hydrophone,  produces  the  same  statistics  for 
the  reference  and  primary  inputs  if 

(a)  the  single  primary  hydrophone  is  placed  Ld  feet  from  the  ref¬ 
erence  array 

(b)  the  primary  hydrophone  has  the  same  directional  response,  H(d) 
as  the  primary  array 

(c)  the  ambient  noise  component  of  the  primary  hydrophone  output 
has  spectral  density 
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where  an  is  the  shading  applied  to  the  nul  hydrophone  ando^(u))  is  the 
spectral  density  of  the  output  of  a  single  hydrophone  on  in  the  primary  array 
being  replaced.  Hence,  in  the  computation  of  the  spectral  density  of  the 
primary  hydrophone  output  and  the  cross -spectral  density  between  the  pri¬ 
mary  and  the  references  (as  required  to  evaluate  the  canceller  performance. 
Section  2.3),  the  source  CSD  function  Sj(d,  ui)  i.e.  modified  by  the  response, 
H(d). 

2.3  Model  for  the  Adaptive  Canceller.  In  this  study,  the  performance 
of  the  adaptive  canceller  is  analyzed  under  the  assumption  that  the  inputs 
to  the  canceller  are  stationary  random  processes.  This  requires  that  the 


26 


pi  lllkif  y  dl  I  ay 

steering  angle 


°N-I 


pr imary 
input 


Figure  2-5.  Use  of  a  Line  Array  as  the  Primary  Input 


Figure  2-6.  Equivalent  Model  for  Primary  Hydrophone  Array 
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weights  only  be  adapted  during  the  time  when  the  reverberation  is  present. 
Section  2.1  described  a  canceller  configuration  in  which  the  weights  are 
frozen  except  when  a  detector  indicates  that  the  CZ  surface  reverberation 
is  present.  The  analysis  which  is  discussed  here  therefore  applies  to  this 
canceller  structure  when  the  detector  determines  the  presence  of  the  re¬ 
verberation  with  no  errors.  Also  described  in  Section  2.1  is  a  second  con¬ 
figuration  of  the  adaptive  canceller  to  reverberation  in  which  the  weights 
continue  to  adapt  when  only  ambient  noise  is  present.  The  performance  of 
this  canceller  configuration  will  therefore  be  degraded  somewhat  from  the 
performance  predicted  using  the  stationary  input  canceller  model.  If  the 
adaptation  in  ambient  noise  is  very  slow  and  the  length  of  time  during  which 
only  ambient  noise  is  present  is  relatively  short,  then  the  results  using  the 
stationary  model  will  be  a  good  approximation  to  the  actual  performance  of 
the  second  canceller  configuration .  The  degree  to  which  the  performance  of 
these  two  canceller  configurations  is  degraded  from  the  stationary  model  is 
discussed  in  Section  2.7  and  analyzed  in  detail  in  Appendix  J. 

The  multiple  LMS  adaptive  canceller  structure  of  Figure  1-1  is  analyzed 
by  replacing  each  of  the  LMS  adaptive  filters  with  a  continuous  Wiener  filter. 
In  a  stationary  environment,  the  LMS  adaptive  filter  converges  in  the  mean 
to  the  discrete  Wiener  filter.  In  turn,  the  discrete  Wiener  filter  will  be  a 
good  approximation  to  the  continuous  Wiener  filter  if  the  sampling  rate  high 
enough  and  the  filter  long  enough  for  the  discrete  impulse  response  to  ap¬ 
proximate  that  of  the  continuous  filter.  This  will  be  the  case  in  any  well 
designed  adaptive  canceller.  Therefore,  the  mean  of  the  converged  LMS 
adaptive  filter  weights  (which  comprise  its  impulse  response)  will  be  well 
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approximated  by  the  continuous  Wiener  filter.  The  actual  weights  (impulse 
response)  of  the  converged  LMS  adaptive  canceller  randomly  fluctuate  about 
the  mean  weights,  and  these  fluctuations  increase  the  canceller  output  power 
above  that  of  the  discrete  Wiener  filter.  However,  these  fluctuations  can  be 
made  arbitrarily  small  in  a  stationary  environment  by  making  the  feedback 
coefficient,  n ,  small.  The  use  of  the  continuous  Wiener  filter  to  model  the 
LMS  adaptive  canceller  can  therefore  be  regarded  as  neglecting  the  effects 
of  sampling,  finite  filter  length,  and  algorithm  noise,  all  of  which  will  be 
small  in  a  well-designed  canceller.  Obviously,  the  Wiener  filter  analysis  also 
applies  only  to  the  converged  LMS  canceller,  and  provides  no  information  as 
to  its  transient  response. 

Appendix  C  derives  the  transfer  functions  of  the  Wiener  canceller  illus¬ 
trated  in  Figure  2-7  in  terms  of  the  second  order  statistics  of  the  hydro¬ 
phone  outputs.  These  are  characterized  in  the  frequency  domain  by  S1;L(w) 
the  power  spectral  density  of  the  primary  hydrophone  output,  Se(u>),  the 
cross -spectral  density  (CSD)  matrix  of  the  K  reference  hydrophone  outputs, 
defined  as 

S0(w)  =  (se(p,q,u)j  p  =  0,  1,  . . .  K  - 1  (2-15) 

q  =  0,  1,  ...  K-l 

it.  it. 

with  se(p,q,w)  the  CSD  between  the  outputs  of  the  p  1  and  q  hydrophones 
as  given  in  Section  2.2,  and  by  s^u),  the  CSD  vector  between  the  primary 
and  the  references.  With  the  primary  hydrophone  is  md  feet  from  the  near¬ 
est  reference,  as  shown  in  Figure  2-2, 

Su(w)  =  sg(K  +  m-l,  K  +  m-1,  w)  (2-16) 
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and 


Sj(w)  =  i->eiO  ,  K+m  - 1) ,  se(  l,K+m-l),  ....  se(K  -  l,K+m- 1)  ]T 

(2-17) 

In  this  notation,  the  transfer  functions  of  the  K  filters  in  Figure  2-12  are 
given  by 

H°(w)  =  [H0(u.),H1(w),  ....  Hk_1(to)]T  =  Se1(u)s1(w) 

(2-18) 


and  the  power  spectral  density  of  the  canceller  output  is 

EK(m,w)  =  su(w)  -  s^(cj)  Sg1(uJ)  s^w)  (2-19) 

where  the  subscript  K  on  E^m.o;)  indicates  that  K  references  are  used  in 
the  canceller  and  m  is  the  distance  from  the  reference  array  to  the  canceller 
in  units  of  d  feet.  Since  the  models  described  in  Section  2.2  provide 
se(P-q.  w),  (2-15),  (2-16),  and  (2-17)  can  be  used  to  evaluate  the  cancel¬ 
ler  output  spectrum  (2-19).  This  is  the  primary  vehicle  for  the  analyses  re¬ 
ported  here. 

2.4  Results  From  Linear  Prediction  Theory.  Proceeding  from  the 
general  models  developed  in  Appendices  A  and  B ,  and  from  the  canceller 
model  derived  in  Appendix  C ,  Appendix  D  shows  that  the  Least  Mean  Square 
cancellation  of  a  spatially  uncorrelated,  narrow  extended  source  is  equiva¬ 
lent  to  the  linear  minimum  mean  square  error  prediction  of  a  bandlimited 
random  process  in  noise.  This  equivalence  allows  direct  application  of  known 
results  from  linear  prediction  theory  to  the  extended  source  cancellation 
problem  of  interest  in  this  report.  The  most  useful  of  these  results  involve 
the  irreducible  canceller  output  spectrum,  defined  as 
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E^tm.w)  =  lim  EK(m,u>)  (2-20) 

K  -o» 

Recall  from  Section  2.1  that  the  acoustic  field  model  requires  that  the  plane 
wave  assumption  be  valid  across  the  entire  receiving  array,  which  will  be 
violated  as  K—<»( since  the  length  of  the  reference  array  becomes  infinite). 
Therefore,  Eootm.u’)  must  be  regarded  as  a  lower  bound  to  the  canceller 
output  spectrum. 

When  m=l,  so  that  the  entire  array  is  uniformly  spaced  with  spacing,  d, 
the  irreducible  canceller  output  spectrum  is  shown  to  be  given  by 


V1*") 


/ 


ir 

log  He(a,u ))  da 


ir 


(2-21) 


where  He(a,ui)  is  the  wave  number- frequency  spectrum  (WNFS)  of  the  hy¬ 
drophone  outputs.  For  a  narrow  source  of  angular  extent,  2$  ,  this  is  related 
directly  to  the  source  spectral  density,  S ^ ( <J> ,  u> ) »  as  follows 

OP 

He(o,w)  =  <r^(u)  -t-  ^  i  Sj  (<5^L  +P0>W)  Rect 
k=-« 

where 

V=  —  sin^»  and  W  = 

A  ^ 

When  w<l/2,  (2-22)  reduces  to 

R“*[fW*]  <2'23> 

\ 

The  condition  that  W<  1/2  therefore  eliminates  spatial  aliases  from  the 
WNFS,  and,  as  such,  can  be  regarded  as  a  sampling  cirterion. 


a  +  2irk 


2ir 


vv 


ae[  -ir.ir] 
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An  important  footnote  to  equation  (2-21)  is  that  when  log  [H^a.w)]  is 
not  integrable,  the  right  hand  side  must  be  replaced  by  zero.  Appendix  D 
shows  that  if  there  is  no  ambient  noise  and  the  criterion  y*<l/2  is  met. 
then  Egotl.w)  =  0,  so  that  at  least  asymptotically,  the  extended  source  can 
be  cancelled  to  zero.  However,  when  ambient  noise  is  present. 

EooC l,w)  >ct^(w)  (2-24) 

so  that  the  extended  source  can  never  be  cancelled  to  the  ambient  noise 
floor. 

When  there  is  no  ambient  noise  present,  not  only  can  the  irrducible  can¬ 
celler  output  power,  E^d.u),  be  shown  to  be  zero,  but  the  canceller  output 
spectrum,  Ek(1,cj),  can  be  shown  to  approach  zero  asymptotically  in  K,  the 
number  of  references.  That  is,  Ek(1,cj)  is  asymptotic  to 

L(K)  =  G (uj) [ sin (7ry*) ]  2K+ 1  (2-25) 

where  G(w)  depends  upon  the  source  extent,*,  the  sampling  parameter,  y, 
and  upon  the  WNFS  of  the  hydrophone  output,  but  not  upon  K  (G(<->)  is 
given  by  equation  D-61  of  Appendix  D).  The  rate  at  which  this  asymptotic 
approaches  zero  does  not  depend  upon  the  shape  of  the  source,  as  can  be 
seen  from  (2-25).  It  is  pointed  out  in  Section  2.5  that  even  for  relatively 
small  K,  L(K)  is  a  good  approximation  to  E^d.w).  Further,  L(K)  will  be 
useful  in  predicting  the  value  of  E^(l,w)  when  ambient  noise  present,  as 
developed  in  Appendix  K . 

The  more  general  case  of  m  >  1  is  also  considered  in  Appendix  D.  It  is 
shown  that  if  the  irreducible  canceller  output  spectrum  for  m=l,  that  is 
Eood.w),  is  zero,  then 
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Ea;(  m .  w) 


0 


m  >  1 


(2-26) 


This  means  that  if  the  output  spectrum  gives  asymptotically  to  zero  with  K 
for  m=l,  then  it  will  do  so  for  m  >  1,  too.  Although  there  is  no  result 
analogous  to  (2-25)  for  m  >  1  indicating  exponential  behavior,  numerical 
results  discussed  in  Section  2.5  suggest  that  E^(m,ui)  also  goes  to  zero  ex¬ 
ponentially  but  at  a  slower  rate.  As  in  the  case  of  m=l,  the  irreducible  out¬ 
put  spectrum  is  always  non- zero  when  ambient  noise  is  present,  so  this 
result  applies  only  to  the  noise  free  case.  It  is  also  shown  that  E<n(m,uj)  is 
a  non -decreasing  sequence  in  m,  so  that  the  irreducible  output  spectrum, 
which  places  a  floor  on  the  cancellation  performance,  cannot  decrease  as  the 
references  are  moved  farther  from  the  primary  (m  increases).  In  fact,  for 
the  cases  of  interest  here,  E«o(m,u)  always  increases  with  m. 

Appendix  D  also  develops  an  implicit  expression  for  the  cancellation 
floor,  E^tm.w),  in  terms  of  the  source  spatial  distribution,  S ^ (<J> . t*>) .  It  is 
shown  that 


m- 1 


E<o(  m ,  w) 


■  s 


k=0 


Ck(w) 


(2-27) 


where 


Co(w)  =  exp^ 


±  f 

2  rr  J 


log  [H  (a.u;)]  da 


(2-28) 


and  the  remaining  C^'s  satisfy 

o°  fb  (< jj )  °° 

£  C*(w)  Zn  =  expj— +  £  bn(«)  Zn 
n=0  I  n=l 


(2-29) 
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with 


bn(w) 


2/r 


I 


7 r 

log  [He(a,aj)] 


e‘jn® 


do 


(2-30) 


-  7 r 

Again,  He(a,w)  is  the  hydrophone  output  WNFS  given  by  (2-22)  or  (2-23) 
in  terms  of  Sj(0,oj).  The  implicit  relationship,  (2-29),  can  be  solved  explicitly 
for  small  m  (see  Appendix  D,  supplement  II),  but  the  expressions  are  quite 
complicated  and  give  little  insight  into  the  canceller  performance. 

2.5  Results  for  a  Uniform  Spatially  Uncorrelated  Interference.  Suppose 
that  the  model  shown  in  Figure  2-3  is  used  to  represent  the  surface  rever¬ 
beration  from  the  CZ.  A  reasonable  model  for  the  reverberation  is  that  the 
intensity  of  the  reverberation  does  not  vary  with  bearing  in  the  mainlobe 
of  the  transmitted  signal,  so  that,  using  (2-3) 


SIa(fl,cj)  = 


<^(“0,  9e[0o  -  ev  eQ  +  ej 

0,  otherwise 

V* 

where  6  is  assumed  small  in  the  sense  of  (2-4).  Incorporating  two  way 
cylindrical  spreading  loss  into  the  source  characteristics  suggests  that 


(2-31) 


l 


SIr(P’w)  =  2p 


(2-32) 

so  that  the  reverberation  falls  off  linearly  with  range.  Then  from  (2-6), 


Gr(u>) 


(2-33) 


and  the  hydrophone  output  CSD  is 
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cjd 

O  -j-e-  (p-q)  cosd0 

Se(p,q,u;)  -  V",4pq  +  e  4 


R2-R2 
R0  RI 


^j2(w) 


sin  ^  (p-q)  3 2  sin  sq 


?f  (p'q)  sin*o 


(2-34) 


which  is  just  a  straight  forward  extension  of  (2-11)  using  (2-5).  Hence 
the  uniform  source  model  developed  in  Appendix  A  can  be  used  to  investi¬ 
gate  the  performance  of  the  canceller  with  CZ  surface  reverberation  if  the 


2  2 
R0RI 


/  4. 


source  spectrum ,  or  j  (w) ,  in  (2-10)  is  replaced  byff"(u>) 

Appendix  F  uses  uniform  source  model  and  the  results  of  Appendix  D 
to  determine  the  cancellation  floor,  E^,(l,m),  for  the  CZ  surface  reverbera¬ 
tion  modelled  as  in  (2-34).  A  significant  result  is  the  existence  of  a  spatial 
sampling  criterion  relating  the  hydrophone  spacing 


Y=  xsin% 

and  the  source  extent,  0^  It  is  shown  that  if 

yei*l 

then  only  very  slight  cancellation  is  possible.  One  the  other  hand,  if 

ye  <1 

then  the  cancellation  floor  is  given  by 

2W 


(2-35) 


(2-36) 


(2-37) 


Eocd.w)  =  o*(u0 

2VV 


1  + 


INR 


,  ye1<2 


(2-38) 


where 


36 


1NR 


(2-39) 


is  the  interference  to  noise  ratio  at  the  array,  and 

w  =  y  e{  (2-40) 

2 

Using  (2-38),  the  cancellation  floor  relative  to  the  input  power,  cr^f  ^)(  1+INR) . 
is 


E,*,  (  1  ,tj)  -if  INRl  1 

~n -  -  [1+INR]  1  +  -W  ,w<4  (2-41) 

yu)(  1+INR) ,  2W  "  -*  “ 


This  is  plotted  in  Figure  2-8  for  values  of  W  between  0  and  0.5,  and  for 

INR's  of  40,  30,  and  20  dB  .  It  can  be  seen  that  only  when  Wis  small  is 

cancellation  to  near  the  ambient  noise  floor  possible. 

A  second  useful  result  derived  in  Appendix  F  for  the  uniform  source 

case  is  an  asymptote  in  K  for  E^(l,  j)  when  the  ambient  noise  is  not  present. 

This  is  developed  using  the  form  of  the  Strong  Szego  Limit  Theorem  [12] 

2 

discussed  in  Appendix  D.  It  is  shown  that  for  o  (  w)  =  0  and  W<l/2,  E,  (1,  ,  ) 

n  k 

is  asymptotic  in  K  to 


A(K) 


20  a1  (w) 
e 

2W 


[sinirW]  2K  +  1 


(2-42) 


that  is  asymptotically,  E^O,  w)  decreases  asymptotically  with  K.  Although 
the  noise  free  case  is  not  of  practical  interest,  it  will  be  shown  later  that 
A(K)  is  useful  in  describing  the  canceller  behavior  in  the  noise  present  case. 

The  results  of  Figure  2-8  indicate  that  the  cancellation  can  be  made 
arbitrarily  close  to  the  noise  floor  by  making  d  small,  which  moves  the 
references  closer  together.  However,  it  should  be  noted  that  as  d  gets 
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smaller,  the  references  are  also  moving  closer  to  the  primary  in  this  model, 
accounting  for  some  of  the  cancellation  improvement.  A  meaningful  examination 
of  the  performance  therefore  requires  evaluation  of  E^(m,0,  where  the  refer¬ 
ence  are  md  feet  from  the  primary .  This  would  allow  determination  of  the  cancel¬ 
lation  floor  as  the  references  are  moved  closer  together  but  maintained  a  con¬ 
stant  distance  from  the  primary.  Although  Appendix  D  developed  an  expres¬ 
sion,  (2-29),  which  allows  explicit  determination  of  Eao(m,cj),  the  expres¬ 
sions  are  very  complex  and  yield  no  insight. 

An  alternative  approach  is  to  evaluate  EK(m,w)  numerically,  using  the 
results  developed  in  Appendix  E.  Using  the  fact  that  the  eigenvalues  and 
eigenvectors  of  the  reference  hydrophone  output  CSD  matrix  can  be  ex¬ 
pressed  in  terms  of  the  Discrete  Prolate  Spheroidal  Sequences  (DPSS's)  and 
their  associated  eigenvalues  [13]  ,  it  is  shown  that 


EK(m,«)  =  +  29i<re(u;) 


K-ir 

■hi 


Ji 

i 

[VK+m-l  <K'W>] 


261Ak(K,W) 


2Wo^(w)  +  2  01o-g(o;)^(K,W) 


Ak(K,W) 


(2-43) 


where  k  is  the  number  of  references,  as  before,  and 


{  V^k)  (K,W)} 


=  kth  DPSS 


m= . . .  - 1 , 0 , 1 ,  . . 


Xk(K,W)  =  eigenvalue  associated  with  the  kth  DPSS. 

The  DPSS  and  their  eigenvalues  are  easily  evaluated  on  the  computer, 
which  allows  numerical  evaluation  of  E^Cm.w)  from  (2-42).  Appendix  K 


L 
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includes  extensive  plots  of  EK(m,'d  as  a  function  of  K  for  various  m  and  values 
of  W  from  10  ^  to  5  x  10  *.  Interference-to-noise  ratios  of  40  dB ,  30  dB  ,  and 
20  dB  are  given,  as  well  as  some  results  in  the  ambient  noise  free  case. 

In  the  ambient  noise  free  case,  plots  show  that  the  asymptote,  A(K),  given 
in  (2-42),  is  an  excellent  approximation  to  E^d.tu)  even  for  small  K,  that  is 


Ek(  l,uj) 


2  VgOd 

2W 


[sin 77  W]2K+1, 


on(w)  =  0 

,W<± 


(2-43) 


Recall  from  Section  2.4  that  an  extended  source  cannot  be  cancelled  to  the 

o 

ambient  noise  floor,  but  approaches  some  cancellation  floor,  E  (l,w)  >c  (jj). 

oo  n 

However,  it  seems  reasonable  that  the  noise  present  case  should  behave 
approximately  like  (2-43)  if  the  interference  to  noise  ratio  is  high  and  the 
cancellation  achieved,  Ek(1,oj)  is  well  above  the  cancellation  floor,  E^d.oj). 
This  conjecture  is  shown  to  be  valid  by  comparison  with  numerical  results 
later. 

The  numerical  results  from  Appendix  K  for  the  ambient  noise  present 
case  are  typified  by  Figures  2-9  through  2-11,  which  show  the  40  dB  INR 
for  W  =  10  ^ ,  10  2,  and  10  **  respectively,  and  by  Figures  2-12,  2-13.  and 
2-14  which  show  the  same  cases  for  INR  =  20  dB  .  Also  shown  on  these  fig¬ 
ures  (and  the  others  in  Appendix  K)  are  the  cancellation  floor,  E^d.w), 
computed  using  (2-41)  and  the  approximate  value  of  E^d.w)  in  the  noise 

free  case  from  (2-43).  Note  that  Ea(l,u)  is  indistinguishable  from  the 

_  o 

noise  floor  for  W  =  10  .It  can  be  seen  that  the  canceller  output  spectrum. 


Ek(1,w),  approaches  the  floor  quite  rapidly  at  first,  following  the  asymptote. 
(2-42).  However,  there  is  a  pronounced  inflection  point  above  which  the  can¬ 
cellation  floor  is  approached  slowly.  This  point  is  very  significant  from  the 
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mu 
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point  of  view  of  the  cost  of  cancellation  improvement  per  dB  .  Suppose  the 
inflection  occurs  at  K  =  Kg.  Then  there  is  a  definite  "law  of  diminishing 
returns"  in  using  more  than  Kq  references,  since  each  additional  reference 
incurs  the  computational  cost  of  another  LMS  adaptive  filter,  regardless  of 
the  additional  cancellation  it  provides.  For  example,  in  Figure  2-15,  use  of 
2  references  achieves  about  33  dB  of  cancellation,  while  use  of  6  additional 
references  only  improves  the  cancellation  3-4  dB . 

It  is  shown  in  Appendix  K  that  this  inflection  point  occurs  at 
approximately 

Kfl  =  [Kj]  (2-44) 


where  [x]  is  the  smallest  integer  greater  than  or  equal  to  x  and 

Kj  =  |  [log(sin  it  W)]  _1  [  (2W-1)  log  —.]  -1  ,  W<|  (2-45) 

The  parameter  Kj  is  just  the  intersection  of  the  asymptote.  (2-42)  and  the 
cancellation  floor  (2-41).  Choosing  K  =  Kq  (for  m=l)  does  not  mean  that 
further  cancellation  is  not  possible,  but  that  this  additional  cancellation  is 
costly . 

The  number,  Kj,  provides  insight  as  to  how  the  number  of  references 
required  in  a  given  situation  changes  with  W  and  with  INR.  From  (2-45), 


it  can  be  seen  that  Kj  is  linear  in  10  log  (INR),  which  is  the  interference- 
to-noise  ratio  in  dB .  The  slope  of  Kj  with  respect  to  10  log  (INR)  is 


(2-46) 
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so  that  approximately  ( 2VV- 1 )  /(  20  log  (simW))  references  are  required  per  dB 
of  INK.  Figure  2-15  shows  Kj  plotted  as  a  function  of  10  log  (INR)  for  vari¬ 
ous  W.  The  dependence  on  the  interference-to-noise  ratio  is  quite  weak 
when  W  <^0.01,  that  is,  when  the  interferences  are  very  narrow. 

The  dependence  of  Ij  on  W  is  more  complicated  as  shown  in  Figure  2-16, 

which  shows  Kj  plotted  as  a  function  of  W  for  INR  =  20,  30,  and  40  dB .  For 
-2 

small  W  (say,  W  <10  )  the  dependence  of  on  W  is  very  weak,  regardless 

-2 

of  INR.  As  W  increases  above  10  ,  the  value  of  Kj  begins  to  increase 

rapidly  with  W,  going  asymptotically  to  infinity  at  W  =  1/2.  Recalling  that 
no  cancellation  is  possible  when  W  =  1/2,  it  would  be  expected  that  -+00 
at  this  point. 

Now,  suppose  the  criterion  K  =  KQ  is  used  to  select  the  number  of 
references  used  in  a  given  situation.  Figure  2-17  shows  that  the  cancella¬ 
tion  achieved  is  within  several  dB  of  the  cancellation  floor,  (2-41)  for  this 
choice  of  K.  In  fact,  it  can  be  seen  that  the  value  ofE^  (l.o.)  is  very  nearly 
log  (INR)  dB  above  the  floor,  so  that 

2VV 

Er  (1,w)  *  (INR)0' 1  Ew(  !,<*>)  =  2ei(X“(u)(INR)'-9  [l  + 

(2-47) 

provided  that  (INR)0,1  ( 1 , io)  20j<j^  (ui)  +  (oj).  This  condition  just 

assures  that  the  approximation  is  not  used  when  it  would  produce  an  arrar- 
ent  increase  in  the  interference  power  relative  to  the  uncancelled  primary 
output.  If  INR  >>1  and  W  <  1,  as  it  will  be  in  all  cases  of  interest,  (2-47) 
can  be  approximated  as 

Ek  (1,u)  as  28.af  (<u)  ( INR ) 2W' '  9  (  2W)  "2W  (2-48) 

1  e 
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For  0<  IV  <  1,  (2W)  “  is  within  about  1.6  dB  of  unity,  so 


o  2VV-  9 

E  ( 1 ,  j )  «  26  -  U)  (INR)  • 

A0  1  e 


(  2-49) 


The  cancellation  achieved  is  then  approximately 


CK  (l,w) 
K0 


Ev  (1  ,<o) 

K0 


(INR) 


2W+0.1 


1  +  INR 


(INR) 


2W-.9 


for  INR  >>  1.  Expressed  in  dB ,  there  yields 

10  log  C„  (l.io)  a  -0.9  [10  log(INR)]  +  2W  [10  log  INR]  (2-51) 

0 

The  interference  rejection  in  dB  varies  approximately  linearly  with  both 
the  interference  to  noise  ration  in  dB  and  the  paramter  W.  This  result  is 
valid  over  the  range  10  ^  <  Ws  0.45  and  102s  INR  s.  10^. 

Examination  of  the  numerical  results  of  Appendix  K  shows  for  any  given 
selection  of  K  and  W,  the  curves  of  Ek(m,j)  for  m  >  1  have  roughly  the  same 
shape  as  that  of  Ek(l,a))  and,  in  particular,  have  the  same  inflection  point 
in  K.  Therefore,  the  value  of  Kq  given  above  is  suitable  as  a  guide  to  the 
selection  of  the  number  of  references,  K,  even  when  m  >  1.  For  some  values 
of  W,  however,  Ek  (m,uO,  m  >  1,  falls  off  somewhat  faster  with  K  >  KQ  than 
E^  (l.w)  does.  For  example,  in  Figure  2-10,  with  m=l  and  KQ=2  yielding 
about  33  dB  of  cancellation.  Six  additional  references  (K=8)  gives  only 
3-4  dB  more  cancellation,  so  choosing  K  =  Kq  seems  justified.  However, 
for  the  same  case  with  m=5,  K  =  Kq  =  2  provides  22  dB  of  cancellation,  while 
using  only  2  more  references  gives  6  dB  more  cancellation.  Thus.  Kq  appears 
to  be  a  rougher  guideline  when  it  is  used  for  m  >  1.  Of  course,  there  will 
always  be  some  ambiguity  as  to  how  much  rejection  must  be  provided  by  a 
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reference  to  justify  its  use.  Certainly  in  any  given  design  situation,  Kq 
can  be  used  as  a  "rule  of  thumb,"  with  final  choice  of  K  made  by  referral 
to  the  curves  given  in  Appendix  K . 

The  numerical  procedure  for  calculating  E^.  (m,-)  as  given  in  Appendix 
K  allow  assessment  of  the  effects  of  changing  the  hydrophone  spacing,  d, 
while  holding  the  distance  from  the  primary  to  the  closest  reference  (md) 
constant.  If  d  is  divided  by  4  (reducing  hydrophone  spacing),  then  m  is 
multiplied  by  A  to  make  md  constant.  The  results  of  this  are  shown  in  Fig¬ 
ures  2-18,  2-19,  2-20,  and  2-21.  These  curves  show  E^(m,uj)  as  a  function 
of  K  for  various  values  of  A,  where 
m  =  A 

and 

«yA) 

W  =  W  /A  =  — -  6L  sinfl  (2-52) 

O  A  X  O 

Hence,  the  curves  can  be  interpreted  as  varying  spacing  d,  while  holding 
md  constant.  Figures  2-18,  2-19,  2-20,  and  2-21  use  Wq  =  .1,  .05,  .005,  and 
.001  respectively,  and  vary  m  over  the  range  indicated  on  the  figures.  The 
plots  are  all  for  an  INR  of  40  dB .  Recall  that  as  m  was  increased  for  fixed  W 
and  K  (that  is,  as  the  reference  array  is  moved  farther  from  the  primary  with 
fixed  spacing  between  reference  hydrophones)  the  canceller  output  spectrum 
increased  monotonically .  The  situation  is  quite  different  when  the  reference 
spacing  is  changed  holding  the  distance  to  the  primary  constant.  The  be¬ 
havior  is  not  monotonic,  and  the  E^  (m,aj)  may  increase  or  decrease  depend- 
int  upon  WQ  and  K.  The  details  of  this  behavior  cannot  be  predicted  from 
the  analytical  results  at  present,  and  will  be  investigated  further  in  the  next 
phase  of  the  study. 
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Figure  2-20 

Canceller  Outputs  Spectrum  vs  K 
WQ  =  .005 

1  NR  =  1»0  dB 


Figure  2-20.  Canceller  Output  Spectrum  vs.  K  for  W  =  0.005,  I  NR  =  40  dB 
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Figure  2-2! 

Cancel ler  Output  Spectrum  vs  K 

W  =  .001 
o 

I  NR  =  40  dB 


Figure  2-21.  Canceller  Output  Spectrum  vs.  K  for  W  =  0.001,  INR 


2.6  Cancellation  for  an  Extended  Interference  Producing  an  Exponential 
Hydrophone  Output  CSD.  It  is  quite  common  in  the  study  of  bandlimited 
temporal  random  processes  to  model  the  process  as  having  an  exponential 
correlation  function. 

p(T)  =  or2  e“B*T*  (2-53) 

2 

where<r  is  the  power  of  the  process  and  where  B  is  the  bandwidth.  This 
correlation  function  has  the  advantage  of  being  representative  of  many 
physical  processes  while  having  nice  analytical  properties.  In  particular, 
given  a  sequence  of  N  samples,  uniformly  spaced  m  time,  of  such  a  process, 
the  correlation  matrix  of  the  N  samples  is  invertible  explicitly  (this  is  dis¬ 
cussed  thoroughly  in  Appendix  G).  It  therefore  seemed  reasonable  to  con¬ 
sider  a  source  producing  an  exponential  hydrophone  output  CSD  and  the 
performance  of  the  canceller  with  such  an  interference. 

Appendix  G  assumes  the  same  hydrophone  geometry  as  Figure  2-2,  but 
assumes  that  the  extended  interference  produces  the  hydrophone  output  CSD 

se(p,q,w)  =  <r2( «)  6pq  +<x2(&»)  e  b(w>Ip  qld  (2-54) 

iL  th 

between  the  outputs  of  the  p  and  q  1  references,  where  B(w)  is  a  spatial 
bandwidth.  Using  the  results  of  Appendix  D  and  the  special  properties  of 
the  CSD  matrix  produced  by  (2-54),  the  irreducible  canceller  output  power, 
E^m.oO,  can  be  determined  in  closed  form  for  arbitrary  m.  It  is  shown  that 


Eq,;  m ,  g>  )  =  <Tj(<a) 


/3  ( 1-  P2m)  +  (1 .-  P2m) 


2/3 2Pm  ’ 

)3(  l+p2)+(  1-p  2)+\  [)3(  1+P2)+(1->P2)]  2-(  2/3p ) 2 


■> 

► 


(2-55) 
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where 


p  =  e-B(6,)d 


(2-56) 


and 


(3  = 


0-^(0)) 


(2-57) 


In  the  interference  only  (ambient  noise  free)  case,  (3=0  and  (2-55)  re¬ 
duces  to 

E^m.oO  =  o-^a. )(1-P2m)  =  <rx2(M)(l  -  e‘2mB (<U)d) , <r^(«> )  =  0 

(2-58) 

From  Appendix  G,  however,  (2-58)  is  just  the  cancellation  that  would  be 
achieved  using  a  single  reference  placed  md  feet  from  the  primary.  There¬ 
fore,  use  of  more  than  one  reference  (K  >  1)  does  not  improve  canceller 
performance.  This  is  a  fundamental  difference  between  the  uniform  inter¬ 
ference  discussed  in  Section  2.5  and  the  one  producing  the  exponential  CSD 
of  (2-54).  Recall  that  for  the  uniformly  distributed  source 
lim  E„(m, ti>)  =  0 

i ,  iv 

K-*-oo 


The  difference  in  the  behavior  of  E^(m,w)  for  the  two  hydrophone  output 
CSD's  can  be  explained  in  terms  of  the  special  nature  of  the  exponential 
CSD.  It  can  be  shown  that  the  sequence  of  hydrophone  outputs,  taken 
at  a  single  frequency  (in  the  interference  only  case),  are  a  first  order 
Markov  sequence.  As  a  result,  the  optimal  canceller  uses  only  the  closest 
reference,  regardless  of  how  many  are  possible. 

Figures  2-22  and  2-23  plot  the  irreducible  canceller  output  power 
relative  to  the  input  power  given  by  (2-55)  as  a  function  of  p  for  values  of 
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A 


10  log  (1//5)  of  40  and  20  dB ,  respectively.  It  can  be  seen  that  cancellation 

2 

to  near  the  noise  floor ,  <rn  ( u ) ,  requires  that  the  primary  and  references  be 
very  highly  correlated.  This  is  again  due  to  the  Markov  nature  of  the  hydro¬ 
phone  outputs  under  the  exponential  CSD.  At  high  interference  to  noise 
ratio,  the  optimal  canceller  uses  only  the  reference  closest  to  the  primary  and 
the  output  power  is  proportional  to  (1  -p^m).  Additional  references  have  no 
effect,  so  the  canceller  depends  completely  on  the  proximity  of  the  closest 
references  to  the  primary  to  achieve  cancellation.  Hence,  p  must  be  very 
close  to  unity  to  produce  cancellation  to  near  the  noise  floor. 

The  canceller  performance  with  the  exponential  hydrophone  output  CSD 
can  be  compared  with  that  achieved  in  the  presence  of  a  uniform  interference 
(producing  a  sin  x/x  shaped  CSD)  by  noting  that  the  parameter  B(«)d  in 
the  exponential  case  plays  the  same  role  as 

W  =  j-  <t>  sin<£o 

in  the  uniform  case.  Appendix  G  makes  the  following  observations: 

(a)  If  W  11/2  in  the  sin  x/x  case,  virtually  no  cancellation  is 
achievable,  as  already  observed  in  Section  2.5.  No  such  sampl¬ 
ing  criterion  is  present  in  the  exponential  case,  but  cancella¬ 
tion  is  still  severely  limited  for  B(o>)d  >  1/2.  Either  CSD  requires 
a  value  of  W  much  smaller  than  1/2  if  the  canceller  is  to  produce 
cancellation  in  the  10  to  40  dB  range. 

(b)  The  optimal  canceller  produces  better  results  for  the  exponential 
CSD  for  W  larger  than  approximately  0.4.  However,  in  the  im¬ 
portant  region  of  small  W ,  producing  more  than  10  dB  of 
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cancellation,  the  canceller  is  better  able  to  reject  the  sin  x/x 
source.  This  is  due  to  the  Markov  property  of  the  exponential 
source,  as  discussed  above. 

Summarizing-,  the  exponential  CSD  produces  significantly  different  re¬ 
sults  than  the  sin  x/x  case,  as  noted  above.  The  differences  between  the  two 


reflect  a  Markov  property  inherent  in  the  double  exponential  CSD.  It  is 
likely  that  the  differences  also  result  from  the  fact  that  the  sin  x/x  CSD  is 
produced  by  an  angle  limited  source,  while  the  exponential  case  results  from 
a  source  distributed  on  (- n,it ).  Like  the  sin  x/x  CSD,  the  source  producing 
the  exponential  CSD  cannot  be  cancelled  to  an  ambient  noise  floor  for  m=l. 
Further,  the  minimum  canceller  output  power  is  shown  to  increase  with  m, 
just  as  the  numerical  results  indicated  for  the  sin  x/x.  However,  the  param¬ 
etric  behavior  of  the  canceller  with  the  two  CDS's  differs  sufficiently  that 
the  exponential  case  cannot  be  used  to  predict  the  canceller  performance  with 
the  sin  x/x  CSD. 

2.7  Effects  of  Reverberation  Non-stationarity  on  Canceller  Behavior. 

Section  2.1  considered  two  methods  of  applying  the  LMS  adaptive  canceller 
to  the  non-stationary  reverberation  environment.  Both  methods  held  the 
weights  fixed  for  an  interval  of  time  after  transmission  in  order  to  assure 
that  the  canceller  was  not  operating  during  periods  of  reverberation  from 
other  sources.  This  is  critical  because  the  highly  correlated  inputs  to  the 
canceller  during  these  periods  would  severely  degrade  cancellation  of  the 
CZ  reverberation  if  the  weights  were  not  frozen.  Once  the  initial  time 
interval  has  passed  (and  all  other  sources  of  reverberation  died  out),  the 
question  arises  as  to  the  effects  of  ambient  noise  on  canceller  performance. 


The  first  method  discussed  in  Section  2.1  seeks  to  avoid  canceller  adaptation 
in  ambient  noise  completely .  It  does  this  by  using  an  energy  detector  in  the 
reverberation  band  to  determine  when  the  CZ  reverberation  is  present,  and 
only  allowing  the  weights  to  adapt  during  this  time.  Given  the  fact  that  the 
reverberation  dominates  the  ambient  noise  (and  signal,  if  present),  this 
detection  should  be  very  reliable. 

However,  it  is  well  known  that  the  LMS  algorithm  adapts  much  more 
rapidly  in  an  environment  producing  correlation  between  its  inputs  (such  as 
reverberation)  than  to  one  producing  uncorrelated  inputs  (like  ambient  noise 
does).  That  is,  the  filter  "learns"  the  correlation  properties  of  an  input 
much  more  quickly  than  it  "forgets"  those  properties  once  the  correlation 
disappears.  The  second  approach  discussed  in  Section  2.1  takes  advantage 
of  this  by  operating  the  adaptive  canceller  continuously  after  the  initial 
time  interval  has  passed.  Thus,  the  weights  adapt  during  the  period  of 
ambient  noise  alone,  as  well  as  when  the  CZ  reverberation  is  present.  Since 
the  filter  "forgets"  slowly,  the  degradation  due  to  adaptation  in  the  presence 
of  the  ambient  noise  will  be  tolerable  if  the  period  of  ambient  noise  is  not  too 
long.  This  approach  eliminates  the  need  for  (and  cost  of)  the  reverberation 
detector . 

The  latter  approach  is  considered  in  Appendix  J.  A  pulsed  model  for 
the  adaptive  canceller  input  is  developed  then  used  to  derive  an  expression 
for  the  mean  weights  and  man  square  error  (canceller  output  power)  in  such 
an  acoustic  field  in  the  absence  of  signal.  It  is  shown  that  if  the  reference 
hydrophone  cross-spectral  density  matrix  due  to  interference  has  a  dominant 
eigen  value,  X  (as  it  will  for  a  plane  wave  or  narrow  interference  and 
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high  INR).  then  the  transient  response  of  the  system  is  determined  by 

y  ( P  +  \  D  ),  where  P  is  the  ambient  noise  power  and  where  D  is 
n  max  c  n  c 

the  duty  cycle  of  the  reverberation,  i.e.. 


time  CZ  reverberation  is  present  in  one  transmit  cycle 


time  canceller  is  adapting  during  one  transmit  cycle 


Hence,  the  transient  response  is  determined  by  the  noise  power  plus  the 
time  average  of  the  eigenvalue  over  the  adaptation  time. 

It  is  also  shown  that  the  pulsed  character  of  the  input  produces  an 
increase  in  the  steady  state  canceller  output  power  above  that  which  would 
occur  if  adaptation  only  occurred  in  the  presence  of  CZ  reverberation.  Note 
that  the  canceller  reaches  a  steady  state,  time  varying  solution  in  response 
to  the  non- stationary  input.  For  a  single  plane  wave  interference,  this 
maximum  increase  is  shown  to  be 

k-D  2 

[1  -  (l-yP)K  P] 

Y  =  - — I? - , -  (2-59) 

max  [1  -  (l-uPj)P  <i-uPn)k'p] 

where  Pj  is  the  interference  power,  is  is  the  total  number  of  samples  in  the 
period  during  which  filter  adapts  in  one  transmit  cycle,  and  p  is  the  number 
of  samples  within  k  that  the  reverberation  is  present.  This  has  been  shown 
numerically  to  be  a  monotone  increasing  function  of  y  ,  so  the  increase  in 
canceller  output  power  can  be  made  arbitrarily  small  by  decreasing  y . 
However,  it  will  generally  be  required  that  the  canceller  converge  within 
several  pings,  so  that  y  will  be  selected  to  give  the  desired  dynamic 
response.  In  that  case,  (2-59)  will  provide  an  indication  of  the  increase 
in  cancellation  above  that  predicted  in  Section  2.5  and  2.6.  Note  that 
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'max  *s  m*n*m*zed  by  choosing  p  as  close  to  k  as  possible,  i.e. ,  only  adapting 
in  the  presence  of  interference.  It  remains  to  develop  an  equivalent  expres¬ 
sion  to  (2-59)  when  the  source  is  extended,  and  to  use  these  results  to 
generate  design  guidelines  for  the  implementation  of  the  canceller  structure. 

2.8  Spatial  Response  of  the  LMS  Canceller  in  the  Presence  of  an 
Extended  Source.  Although  the  main  thrust  of  this  study  was  the  analysis 
of  the  effectiveness  of  the  canceller  in  rejecting  an  extended  interference 
without  considering  the  effect  of  the  canceller  on  some  plane  wave  signal 
of  interest ,  it  is  clear  that  the  detection  of  such  a  signal  is  the  eventual 
goal.  Appendix  H  therefore  considers  the  spatial  response  of  the  canceller, 
that  is,  its  response  to  a  plane  wave  from  some  direction,  <J>  ,  as  a  function 
of  that  direction.  This  response  is  just  the  power  that  would  appear  at  the 
canceller  output  due  to  a  plane  wave  signal  at  the  angle  <J>  .  In  the  case  of 
a  uniformly  distributed,  spatially  uncorrelated,  plane  wave  source,  this 
spatial  response,  B  (0  ),  is  developed  in  terms  of  the  DPSS's  and  their 
eigenvalues  and  evaluated  numerically  on  the  computer.  Figures  2-24 
through  2-26  show  the  response  of  the  canceller  in  the  presence  of  such  a 
source  centered  at  45°  relative  to  broadside  and  with  an  angular  extend 
of  20°.  The  interference  to  noise  ratio  is  40  dB ,  and  the  hydrophones  are 
uniformly  spaced  1/2  wavelength  apart.  The  Figures  show  the  results 
using  1,  4,  and  8  references,  respectively.  The  spatial  response  can  be 
seen  to  include  a  notch  in  the  direction  of  the  interference ,  as  one  would 
expect.  From  the  Figures  2-24  through  2-26  and  those  given  in  Appendix  H, 
it  can  be  seen  that  the  notch  broadens  and  deepens  as  references  are  added 
until  there  are  3-4  references,  at  which  point  additional  references  only 


66 


sharpen  the  sides  of  the  notch  and  improve  the  spatial  response  to  signals 
in  directions  other  than  that  of  the  interference.  The  results  of  Section  2.5 
indicate  that  for  narrow  sources,  more  use  of  more  than  3-4  hydrophones 
will  not  yield  sufficient  cancellation  improvement  to  justify  their  cost.  There¬ 
fore,  the  results  of  Appendix  H  suggest  that  the  reason  for  adding  more 
references  may  be  to  assure  adequate  signal  response  in  some  direction  of 
interest . 

It  was  shown  in  Section  2.5  that  as  the  distance  between  the  reference 
array  and  the  primary  increase,  the  cancellation  of  the  interference 
decreased.  Figures  2-27  and  2-28  show  the  spatial  response  of  the  canceller 
for  the  same  source  and  an  8  element  reference  array  as  above,  but  with  m  =  10 
and  m  =  20  (recall  that  the  reference  hydrophone  spacing  is  d  and  the  distance 
from  the  primary  to  the  closest  reference  is  md).  It  can  be  seen  that  the 
notch  depth  has  decreased,  indicating  a  degradation  in  the  cancellation, 
and  that  there  has  been  some  degradation  is  the  response  pattern  in  direc¬ 
tions  other  than  that  of  the  interference. 

When  interpretting  these  figures  it  should  be  kept  in  mind  that  the 
primary  consists  of  a  single  omnidirectional  hydrophone,  rather  than  an 
array,  as  would  be  used  in  most  sonar  situations.  As  indicated  in 
Appendix  L,  when  a  primary  array  is  used,  the  responses  like  those  of 
Figures  2-24  through  2-28  would  be  multiplied  by  the  response  pattern  of 
the  primary  array.  This  would  greatly  attenuate  (by  the  primary  array 
sidelobe  level)  those  regions  outside  the  mainlobe  of  the  primary  array. 
Assuming,  then,  that  the  reverberation  that  limits  the  active  detection 
performance  arrives  on  or  near  the  main  lobe,  the  signal  response  in  the 
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vicinity  of  the  notch  is  of  primary  interest.  The  reason  for  adding 
references  would  then  be  to  increase  the  slope  of  the  sides  of  the  notch. 

With  this  in  mind,  it  can  be  seen  that  in  Figures  2-27  and  2-28, 
increasing  the  distance  between  the  references  and  primary  did  not  signifi¬ 
cantly  degrade  the  spatial  response  in  the  vicinity  of  the  notch.  It  can  be 
conjectured  then  that  signal  attenuation  will  not  be  a  primary  consideration 
in  selecting  the  distance  between  the  primary  and  references  when  the 
primary  is  a  steered  array  and  the  interference  is  in  the  mainlobe. 

As  the  interference-to-noise-ratio  decreases,  there  is  a  reduction  in 
the  depth  of  the  null  at  the  location  of  the  interference.  Figure  2-29  shows 
the  spatial  response  for  the  cases  of  Figure  2-26  (m=l,  K=8,  $=10°)  but  the 
interference-to-noise  ratio  has  been  reduced  to  30  dB.  Figure  2-30  repeats 
the  same  case ,  but  with  an  interference-to-noise  ratio  of  20  dB .  It  is 
interesting  to  note  that  the  array  response  in  directions  other  than  the 
interference  direction  is  affected  very  little  by  the  change  in  interference- 
to-noise  ratio. 

A  major  task  of  the  next  phase  of  this  study  is  the  investigation  of  the 
effect  of  the  canceller  upon  a  plane  wave  signal,  and  the  quantification  of 
the  effectiveness  of  the  canceller  in  improving  the  detectability  of  such  a 
signal.  During  the  next  phase,  the  above  conjectures  will  be  examined  and 
quantified. 
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3 . 0  CONCLUSIONS  AND  RE  COMMENDATIONS 


3.  CONCLUSIONS  AND  RECOMMENDATIONS 

This  study  has  considered  the  LMS  adaptive  cancellation  of  CZ  surface 
reverberation  using  a  bistatic  approach,  where  passive  hydrophones  spatially 
separated  from  the  transmit /receive  array  are  used  as  references.  This  was 
done  by  modelling  the  surface  reverberation  as  a  narrow  extended  source  that 
is  spatially  uncorrelated  and  uniformly  distributed.  The  work  to  date  has 
concentrated  on  the  ability  of  the  LMS  structure  to  reject  the  reverberation 
in  the  absence  of  a  plane  wave  signal  of  interest.  It  was  shown  that  by  a 
judicious  choice  of  system  parameters,  such  as  hydrophone  spacing  and 
number  of  references,  the  reverberation  can  be  cancelled  arbitarily  close 
to  the  ambient  noise  flow.  Unlike  the  case  of  a  point  interference,  however, 
if  care  is  not  taken  in  the  placement  of  references,  little  or  no  cancellation 
may  be  possible  regardless  of  the  number  of  references.  Through  a  combina¬ 
tion  of  analytical  and  numerical  techniques,  the  parametric  sensitivities  of 
the  cancellation  performance  to  source  extent,  hydrophone  spacing,  distance 
between  references  and  primary  array,  and  number  of  references  has  been 
characterized.  These  results  comprise  the  basic  tools  for  the  design  of  a 
canceller  for  use  in  a  given  situation. 

An  important  result  was  that  for  narrow  extended  sources,  such  as  CZ 
reverberation,  there  is  a  law  of  diminishing  returns  with  respect  to  the  number 
of  references  to  be  employed.  That  is,  beyond  some  number  of  references, 
say  Kq,  more  references  do  not  provide  sufficient  additional  cancellation  to 
justify  the  computational  cost  of  including  them.  A  design  guideline  for 
selection  of  the  number  of  references  was  derived.  It  appears  that  for  the 
cancellation  of  CZ  surface  reverberation  in  a  sonar  with  a  horizontal 
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beamwidth  of  less  than  20°,  at  most  4-5  references  are  required.  Of  course, 
these  references  only  supply  adequate  cancellation  if  the  other  design  guide¬ 
lines  are  met. 

Since  this  phase  of  the  study  considered  the  behavior  of  the  canceller 
in  the  absence  of  signal,  the  results  presented  here  cannot  be  used  to  predict 
the  improvement  in  the  detection  of  a  plane  wave  signal  in  the  presence  of  the 
CZ  reverberation.  Some  of  the  analysis  of  the  response  of  the  canceller  to 
the  signal  has  been  initiated  in  the  work  reported  on  the  spatial  response 
of  the  canceller.  In  the  next  phase,  this  work  will  be  continued,  with  the 
goal  of  deriving  a  detection  performance  measure,  such  as  deflection,  for  the 
canceller  structure.  The  effect  of  signal  presence  during  adaptation  (to 
"signal  bias"  problem)  will  also  be  investigated. 

The  work  in  this  phase  of  the  study  has  been  based  upon  several 
theoretical  models  of  the  canceller  structure  and  the  CZ  reberberation .  During 
the  next  phase,  computer  simulations  will  be  run  to  validate  both  the  analysis 
and  the  mathametical  models  used.  This  will  include  development  of  a  com¬ 
puter  model  for  the  extended  source  for  use  with  existing  programs  for  the 
multiple  canceller  structure.  If  suitable  sea  tapes  can  be  obtained,  actual 
CZ  reverbation  will  be  compared  statistically  to  the  theoretical  models  used 
here,  and,  if  possible,  run  through  the  multiple  canceller  algorithm. 
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APPENDIX  A:  DERIVATION  OF  THE  FAR  FIELD  MODEL 


APPENDIX  A;  DERIVATION  OF  THE  FAR  FIELD  MODEL 


Consider  the  array  geometry  shown  in  Figure  A-l  located  in  a  noise 
field  consisting  of  an  ambient  noise  background  and  a  strong  extended  inter¬ 
ference  source.  The  array  consists  of  K+l  acoustic  hydrophones  arranged 
as  a  horizontal  line  array.  Of  the  K+l  hydrophones,  K  will  be  denoted 
reference  hydrophones,  as  shown,  and  used  to  cancel  in  the  LMS  sense,  the 
interference  from  the  output  of  the  primary  hydrophone.  The  reference 
hydrophones  are  uniformly  spaced  d  feet  apart ,  while  the  primary  hydro¬ 
phone  is  located  md  ft  from  the  nearest  reference.  The  reference  hydrophone 
outputs  are  denoted  e^(t)  for  k  =  0,  1....K-1,  while  the  primary  hydrophone 
output  is  ed(t) . 

It  will  be  assumed  in  the  development  of  the  noise  field  model  that  the 
ambient  noise  produces  zero  mean ,  wide  sense  stationary  noise  at  the  output 

of  each  hydrophone  with  nk(t)  the  ambient  noise  at  the  output  of  the  k— 

reference  (k  =  0,  1,  ...  K-l)  and  n,(t)  the  output  of  the  primary.  It  is 

d 

assumed  that 

E[nk(t)n1(t+x)]  =  Rn(x)  6kJ  k  =  0,  1,  ...K-l  (A-l) 

1  =  0,  1,  . .  .K-l 

E[nd(t)nd(t+T)]  =  Rd(x)  (A-2) 

and 

E[nd(t)nk(t+x)]  =0  k  =  0,  1,  ...  K-l  (A-3) 


A-l 


Figure  A-l.  Array  Geometry  for  Far  Field  Model 

where  i.  ,  is  the  Kronecker  delta  and  R  (t)  is  the  covariance  function  of 
kl  n 

the  N  ,(t)  and  R  ,(x)  is  the  covariance  function  of  n,(t).  That  is,  the 
d  d  d 

ambient  noise  at  the  output  of  each  hydrophone  is  uncorrelated  with  the 
ambient  noise  at  the  output  of  any  other  hydrophone.  This  assumption  will 
be  reasonable  for  realistic  ambient  noise  fields  if  the  hydrophone  spacing  is 
not  much  less  than  one  half  wavelength  at  the  lowest  frequency  of  interest. 

Let  the  extended  source  be  distributed  in  angle,  <J>,  as  shown  in 
Figure  A-l,  and  assume  that  the  arrival  from  each  angle  is  a  plane  wave 
across  the  length  of  the  array.  Let  the  arrival  from  the  angle  <J>  at  reference 
hydrophone  number  zero,  denoted  i(t,  <j>) ,  be  a  zero-mean,  wide  sense 
stationary  random  process  with 


A- 


E[i(t ,  <j>  )  i  (t+x,  <t>2)  1  =  Rj(4>  j_ . <J>  2.t) 


( A-4) 


If  the  speed  of  sound  in  the  water  is  given  by  c ,  then  the  propagation  delay 
from  reference  number  zero  to  a  point  kd  feet  along  the  line  array  for  the 
arrival  from  angle  <J>  is 


x  (<t>)  =  k  -  cos 


( A-5) 


‘k''  T'  ’  c 

Integrating  over  all  <f>  to  account  for  the  entire  extended  source  gives  the 
hydrophone  outputs  as 


/’■ 


ek(t)  =  n^lt)  +  I  i  (t  -  k-  cos  <j>,  <Jr)  d  d>,  k  =  0,1,  ...K-l  (A-6) 


and 


I; 


e  .(t)  =  n  ,(t)  +  /  i  (t  -  J  -  cos  0,  0)  d  0 


( A-7) 


where 


J  =  K  +  m  -1. 


Assuming  that  the  interference,  i  (t,  0),  is  uncorrelated  with  the 
ambient  noise,  that  is, 


E[i(t,  0)nk  (t  +  t)]  =  0 


and 


E[i(t,  <())nd  (t  +  I)]  =  0 


Vt,  0 

k  =  0,  1,  ...  K-l 


Vx,  0 


( A-8) 


(A- 9) 


A -3 


then  the  cross-covariance  between  the  output  of  the  p—  and  q—  reference 
hydrophones  is  given  by 


Re(p *  q.  x)  =  E[ep(t)  eq(t  +  t)I 
=  RnCr)  <$pq  +f  f 


(pcos  -  qcos<$2  ]  d$jd<J>2 

( A-10) 


i.  U 

Therefore,  the  cross-spectral  density  (CSD)  between  the  outputs  of  the  p — 
and  q—  hydrophones  is 


se(p,  q,  w)  =.>[re(p,  q,  t)] 


=  0n(uj)6pq  +f  J  Sl  (^i*  *2’  “>  e+j  c  tpCOS$l  qCOS<t>2]d4>1d((.2  (A-ll) 


where 

,^[x(t)I  =  J  x(r)e  ^  dt 
-00 

a2(u)  =  ^[Rn(T)] 

and 


Sj(w,  0lf  <J>2)  =t/[RI(x,  4^, 


f  A- 12'» 


(A- 13) 


is  the  spatial  cross -covariance  function  of  the  source  at  frequency  w. 


Similarly,  it  can  be  seen  that 


rx(p,  x)  =  E[ed(t)ep(t  -  x)]  p  =  0,  1.  ...  K-l 


i  /*TT  d 

=  J J  Rj[  x  -  2  (Jcos^j^  -  pcosi?2)]d^1d<j)2 

s1(p,  oj)  =^'tr1(p>  r) ] 


and 


rn(x)  =  E[ed(t)  ed  (t  -  t) ] 

r  r 

=  Rd(x)+^  J  4>2>  t) d <t> j' 


dc}>r 


S11(w)  =^'[r11(T)] 


_  2,  ,  f  f  \  ,  .  .  Jlcosij).  -  cos<j)  „ 

=  od(o))  +  J  j  s I c 4> »  <t>2»  w)  e  c  1  2 

*^7T  7T 


where 

crd(u))  =df[Rd(r)I 

Now,  let  e(t)  be  a  random  vector  of  the  reference  hydrophone 
outputs, 

e(t)  =  [eQ(t),  e^t),  ...  e^ft)]1 


(A  - 14) 


(A-15) 


(A-  16) 

(  A-17) 
1  d<js  ^4*2 


CA-18) 


( A- 19) 
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The  cross -spectral  density  (CSD)  matrix  of  the  reference  hydrophone  output 
is 


Sg(to)  =  '  E[e(  t)e+  ( t)  ]  j 


=  ise(p.  q.  w)\  p  =  o,  l,  ...  k- 1 


q=0,  1,  ...k-1 


(A  -  20) 


The  CSD  vector  between  the  reference  hydrophone  outputs  and  the  primary 
hydrophone  output  is 


SjU)  =  .jr'E[ed(t)e(t)]' 


=  [s^o,  to),  Sl  (l,w),  ..:  SjfK-l,  u) ) ] 


(A-21) 


The  power  spectral  density  of  the  primary  hydrophone  output  is  then  S-^tto) 

It  will  sometimes  be  convenient  to  explicitly  represent  the  number  of 
references,  K,  in  the  notation  for  the  reference  hydrophone  CSD  matrix  and 
the  CSD  vector  between  references  and  primary.  In  these  cases,  Se(K,  to) 
will  be  the  reference  hydrophone  CSD  with  K  references. 


Se(K, 


a))  Sg(u))  |j£  references 


( A-22) 


and  s  ^(K,  to)  will  be  the  CSD  vector  with  K  references 


s^(K,  to)  s^( ^  I k  references 


A -6 


( A-23) 


A .  1  Narrow  Source  Approximation 

When  the  extended  source  is  the  surface  reverberation  from  the 
convergence  zone,  the  extent  of  the  source  will  be  much  less  than  2^  radians. 
Suppose  that  the  source  is  centered  on  an  angle,  <p  ,  and  is  limited  to  the 
sector  [  .}>0  -  $  ,  0o  +  £  ]  ,  so 

SjC'aj,  =  0  (^i  ^  ^  +  1 

<  ( A  -  24) 


^2  i  Uo 


<5  .  <f»0  +  *  1 


Then  (A- 11)  can  be  written 


s  (p,  q,  ui)  =  a  (u)5 
e  n  pq 


+LLsi^  o+v  W)e 


pcos(4>0+4>1)-qcos(d>0+(f)2)j 

CL  -J 


(A -25) 


If  it  is  assumed  that  the  source  is  narrow  enough  that 
sin  <[>=$ 


and 


(A-26) 


cos  $=1 


then 


pcos  (4>0+4>1) -qcos(4>0+4>2)  =  (p-qjcos^-fp^-q^sind^ 


(A-27) 


Therefore 


2 

s„(p,  q,w)  =  6  (4c pq 

c  n 

+j~  (p-q)cos 
+  e 

(A-28) 


<6  /""/**  -j  — 

VV^  '  C 


(p<J>1-q4>2)sin*0 


d*1d«2 


It  can  be  seen  that  under  this  assumption,  there  is  a  two-dimensional  Fourier 
transform  relationship  between  the  hydrophone  CSD  and  the  source  angular 
distribution.  Similarly 


Sjfp.uj)  =  e 


-j^(J-P)cos<b0 


/•$  /•$ 


(J  <t>i-p«t>2)  sin4>0 


d*1d<(>2 


(A -29) 


and 


S11M  *  °'d  M  +  /  /  si(<^o  + 


$9)sin<t> 

si(<^o  +  +  ^2’  ^  e  "  0 


dcf)  1d4>2 


(A-30) 


Using  the  notation  given  in  (A-20)  and  (A-21),  these  expressions  yield  the 
following 


seM  =  <r*  H  I  +  G+(w)  S  (w)  G(w) 


(A -31) 


where 


Seo<“>  * 


?I<V*1’W>#  p-0.1. 


-S> -$> 


q  =  0,  1, 


..  K-l 
..  K-l 

(A -32) 


A -8 


and 


/  -j^-pcos*^ 
=  diag  I  e  f 

r»  '  / 


(A-  33) 


Also 


SjM  =  Ge 


_+-j^-Jcos*0 


5-0 


(A-34) 


where  the  p^-  element  of  is 


\{i/l(V<*>l’V<V~)e 


j  ^(J(t>1-P<i>2)sin40 


dijddig 


(A -35) 


A . 2  Spatially  Uncorrelated,  Narrow  Source 

It  will  be  shown  in  Appendix  B  that  for  reverberation  from  the 
convergence  zone,  the  extended  source  may  be  assumed  to  be  spatially 
uncorrelated,  that  is, 


Sj(w,  4^,  Ag)  _  6  <4>^ — <±>2) 


(A -36) 


where  6 (<t>)  is  the  Dirac  delta  function.  Under  these  assumptions,  (A-28), 
1 

(A-29),  and  (A -30)  reduce  to 


2  -j“(P-<l)cos<^  f*  +r£-(p-q)d>sin<i> 

se(P*  q.'*')  =  °-n  M&pq+  e  4SI<V6>u,)e  d6  (A-37) 

-rf-(J-P)cos  +j  ^  (J-p)<tsind)0 

s1(p,w)  =  e  -$SI(V6,w)e  c  dcp  (A-38) 


and 


Sll(w) 


2  r  <E> 

°d  M  +  Sj(®q  +<t>»  w)d6 


(A -39) 


A -9 


For  the  spatially  uncorrelated  source,  there  is  a  simple  Fourier  transform 
relationship  between  the  source  distribution  function,  and  the  CSD  of 

the  hydrophones. 


Of  particular  interest  will  be  the  case  when  the  source  is  not  only  narrow  and 
spatially  uncorrelated,  but  uniformly  distributed  on  the  sector  (oq-$,  +§, ) 


ST(w,<t>)  =  j<rI2(uj),  4>e  4>Q+  ®j 

0  ,  elsewhere 


I 


(A -40) 


Under  these  conditions,  the  integrals  in(A-37),  (A-38),  and  (A-39)  can  be 
evaluated,  giving 


se(P,  q,  w)  «  an2M6pq 

+  <r2(w)  sin[(P-q)~<Psin<t>Q|  ^(P-qJcos^ 
l[(p-q)^  sin<60] 

[(J-Pf^sinfcj  -:^(J-p)cos<^ 

(J-P)^sin<t>0J 


sin 


s0(P,w)  =  <rI  H 


and 


SnM  =  <rn2M  +  2  «Oj2(w) 


(A -41) 


(A -42) 


(A -43) 


Therefore,  in  the  notation  of  (A-33)  and  (A-35), 


sin  (p-q)  c  ®  sin<t0 

\  [(P-q)-^-ain4>0j 


h  =  o,  i, 

q  =  0,  1, 


...  K-l 

...  K-l 


(A -44) 
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APPENDIX  B;  APPLICATION  OF  THE  FAR  FIELD  MODEL  TO  THE 


SURFACE  REVERBERATION  FROM  THE  CONVERGENCE  ZONE 
B.l  Surface  model  with  straight  line  propagation.  Consider  the  geometry 
shown  in  Figure  B-l,  in  which  the  noise  enamates  from  a  sector  of  an  annulus 
on  the  ocean  surface.  The  annulus  is  centered  at  the  origin  of  the  cylindrical 
coordinate  system,  and  has  inside  and  outside  radii  Rj  and  R  ,  respectively. 
The  sector  is  assumed  to  be  of  angular  width  29^  centered  about  an  angle  0  . 
A  hydrophone  is  located  D  feet  below  the  origin  and  a  second  hydrophone  d 
feet  away  along  the  negative  y  axis,  as  shown.  This  can  be  viewed  as  a 
model  of  the  convergence  zone  surface  reverberation ,  except  that  the  propa¬ 
gation  paths  are  direct  rather  than  along  the  usual  refracted  paths  associated 
with  the  convergence  zone.  In  the  following  results,  it  will  be  shown  that  the 
propagation  paths  from  the  source  to  the  array  do  not  affect  the  hydrophone 
CSD  under  the  plane  wave  and  narrow  source  assumptions  if  the  source  is 
spatially  uncorrelated.  The  assumption  is  that  the  reverberating  surface  is 
uncorrelated  from  point  to  point  is  commonly  made  (see,  for  example,  Van 
Trees  The  hydrophone  below  the  origin  represents  the  transmit /receive 

array  of  the  active  sonar,  and  the  other  hydrophone  is  the  reference.  The 
reverberation  is  limited  to  the  angular  sector,  0  -  9,  to  0  +  0, ,  by  the 
horizontal  directivity  of  the  transmitter. 
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hyd rophone 


Figure  B-l.  Geometry  for  Surface  Model  with  Straight  Line  Propagation 


If  the  ambient  noise  is  uncorrelated  between  the  two  hydrophones 


spaced  d  feet  apart ,  the  CSD  of  their  outputs  can  be  written  as 


se(d,u>)  =  on  (a  )  f  (d) 


ro+ei  rRo[%+h  <*o 

I  I  I  Sl  (ol,02,yl’"2’ 

J  A  -3  •'R  •'A  -A  •'P 


6o  -01  'RI  '6o  -ei  RI 


(B-l) 


.e*\c  [ro(pl’9l}  "  rl(p2’62)] 


plp2dp2de2dold0I 


where 


on  (oj)  -  power  spectral  density  of  ambient  noise 


SI(pl’p2'ere2’  ^  =  CSD  ncdse  radiated  at  (pj.Bj)  and  C  P  2 » ®2> 


rQ(p,e)  =  range  from  point  at  (p,9)  to  hydrophone  zero 


r^p.0)  =  range  from  point  at  (p,9)  to  hydrophone  one 


c  =  speed  of  sound 


This  result  is  just  a  straightforward  extension  of  the  model  developed  in 
Appendix  A. 
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It  is  assumed  that  the  source  is  spatially  uncorrelated,  i.e.  , 


Sj  ( p^,  ^2»  (P^.Q^.w)  o  (p^  -  Og)  6  ( Qj  ”  ^ 

Then  (B-l)  reduces  to 

se(d,oj)  =  Pn2(^)  6  (d) 


( B  —  2 ) 


+j  w  [r  (o,0)  -  r.(p,6)]  „ 
Sj(p,0,w)e  c  p^dpd9 

(B-3) 

Now,  as  in  Appendix  A,  assume  that  the  wavefront  from  each  point  on  the 
annulus  arrives  at  the  hydrophones  as  a  plane,  and  that  the  propagation 
distance  can  be  approximated  as  the  sum  of  the  range  from  the  point  to  hydro¬ 
phone  zero  and  the  plane  wave  propagation  delay.  That  is, 


r..(p,s)»r  (p,s)+dcosa 
1  o 


(B-4) 


where  a  is  the  angle  the  arrival  direction  makes  with  the  axis  of  the  two 
hydrophones.  It  can  easily  be  shown  that 


cos  a  = 


p  cos  9 


(B-5) 


It  is  now  assumed  that  p»D,  which  is  certainly  true  for  the  convergence 
zone  problem,  so 


cos  a<*cos  6 


(B-6) 


] 
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and 


sg  ( d ,  a) )  =  oni'(aJ)6(d)  + 


/9„  +9i  ^  R 

o  1  r  o 

-9.  *'R, 

o  1  1 


+j  —  d  cos  6 

Sj  ( p  ,  0, u))e  p2dpd0 


(B-7) 


or,  making  the  change  of  variables  6'  =  0-9  , 


se(d’w)  =  CTn 


fdl  fRo 

(d)+/  /  SjCp, 

J-f I  -'r 


0'+0Q,tj)e 


+j-^  cos  ( e'+6  )  „ 
c  o  2 


p  dpd@' 
(B-8) 


Note  that  rQ(p,0),  representing  the  long  range  propagation  from  the 
reverberating  surface  to  the  array,  is  not  present  in  (B-8).  Therefore, 
the  use  of  straight  line  propagation  paths  instead  of  refracted  convergence 
zone  paths  should  have  minimum  effect  upon  the  results,  except  as  they  pro¬ 
duce  different  vertical  arrival  angles. 

It  is  now  assumed  that  is  small  enough  that  cos  0^<*  1,  sin  0^  9^, 


so 


cos  (0"  +  0  )  =  cos  0"  cos  9  -  sin  0'  sin  0  **  cos  9  -  0'  sin  0  (B-9) 

o  o  o  o  o 


and 


se(d,(jj)s*an  (go)  6  (d) 


+  e 


(B- 10) 


Sj  (  p ,  0,oo)e 


aid 

c 


0  sin  0 


p  dpd0 
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This  has  exactly  the  same  form  as  equation  (A-25)  of  Appendix  A,  which 
gave  the  CSD  for  a  uniformly  distributed,  uncorrelated  source  in  the  far 
field  model. 


Next,  it  is  assumed  that  the  source  CSD  function,  Sj  (p.e.w),  is  separable 
in  p  and  0,  that  is 

ST  (p,0,uj)  =  ST  (p,<jj)  ST  (0,w)  (B-ll) 

l  l  r  la 

so  that 


se(d,uj)  =  an2<w)  <5  (d) 


j  ^  cos  0 

J  c  o 

R 

6,  . 

f  °  2 

f  1  -37T  9  sin  9o 

+  e 

/  SIr(p,u))p  dp 

_RI 

/  SIa(p,w)e  c  d0 

(B-12) 

Now,  if  the  source  CSD  is  assumed  to  be  uniform  in  angle, 

Sia(9,w)  =  ox2  (w)  (B-13) 

then 


se(d,oj)  =  on  (oj)  6  (d) 


(B-14) 


„  sin  [  2tt01  t-  sin  0  J 

+  Oj  <u>)  Gr(w)  — g - — - —  e 

it  f-  sin  0 

A  O 


j  ^  COS  0 

c  o 
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where 


Gr<“>) 


Ir 


(p,ui) 


2, 

P  dp 


(B-  15) 


If  the  pair  of  hydrophones  on  the  x-axis  is  replaced  by  a  line  array  of 
uniform  spacing,  d,  then  the  CSD  between  the  outputs  of  the  pth  and 
hydrophones  is 

se(p,q,uO  =  on2  (w)  5  (B-16) 

sin  [2ir0  (p-q) -^-sin  6]  )  ~  (p-q)cos  e 

+  o2(u-)  G  (u>)  - - * - —  e  C  ° 

-rtf-  sin  6 

X  o 

where  <5pq  is  the  Kronecker  delta.  Note  that  the  plane  wave  assumption  must 
be  valid  over  the  entire  length  of  the  array. 

This  has  exactly  the  same  form  as  (A-41),  so  that  the  far  field  model  can 
be  applied  to  this  surface  source  under  the  above  assumptions. 

B.2  Vertical  arrival  model.  Consider  the  spherical  coordinate  system 
shown  in  Figure  B-2.  As  before  the  noise  field  is  assumed  to  consist  of 
spatially  uncorrelated  hydrophone  noise  plus  an  extended  source.  Here,  the 
extended  source  is  assumed  to  be  distributed  over  a  horizontal  sector,  0Q-0^ 
<0<0q  +  0^  and  a  vertical  angle,  <<J>  < <J>Q+$  ,  and 

Si  ( 0^,  02>  <P i>  )  =  CSD  between  arrivals  from  angles  (8j,<J>^)  and 

r.(0,<j))  =  propagation  distance  from  source  arriving  at  ( 0,  <J> )  to 


hydrophone  i  (i  =  0,1) 


It  is  assumed  that  the  source  is  spatially  uncorrelated, 


si( 0r 92’ *r ^ 2’^  =  Sj( ©i . i» 92>  °  ^ i_<^ 2^ 


(B-17) 


so  that  the  CSD  of  the  two  hydrophone  outputs  is 


se(d,w)  =  cn  5  (d) 


9^+9i  - 


7 

-9.  J( P 

o  1  Yc 


i|[ro(e,$)-r1(6,0))] 

sin  <j:d<t>do 


(B-18) 


Sj( 0, cf> ,u})  e 


As  in  case  B,1  above,  it  is  assumed  that  the  wavefront  from  each  point 
on  the  source  is  planar,  and  that 


r  9 , 4> )  =  ro(0,4>)  +  d  cos  a 


(B-19) 


In  this  case,  however, 


cos  a  =  cos  9  cos  $ 


(B-20) 


so 


se(d,w)  =  an  (<jj)  5  (d) 


e„+0i  /•^+$ 


°  1  f  o  ) 

I  Sj(  0,  <t>  ,co)e 
-a  7a  _* 


—  COS  0  COS  (J)  „ 

sinZifd<})d0  (B-21) 


'9  -0,  J<t>  -<t> 

o  1  ^o 
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Note  that  the  kernel,  exp  f-j  cos  9  cos  0]  depends  upon  both  variables 
of  integration  and  is  not  separable.  Because  of  this,  reduction  of  (B-21) 
to  a  line  source  requires  that  either  cos  9  1  or  cos  <j>  =  1.  This  means 

that  max  [6o-@,0o+9]  or  max  [  <{>  -$,  .  respectively,  be  small,  which  is 

a  much  stronger  condition  than  in  case  B.l,  which  only  required  that  $  be 
small . 

It  will  be  assumed  that  the  vertical  arrival  angles,  <p£  [  $  - $  >  t>Q+  *  1  >  are 
close  enough  to  the  horizontal  that  cos  4>  <*  1  and  sin  $(%<£,  which  means  that 
the  arrivals  must  be  within  approximately  ±10°  of  the  horizontal.  This  is 
consistent  with  the  arrival  angles  associated  with  the  convergence  zone  in 
surface  ship  hull  mounted  sonars.  Under  this  assumption 

se(d,ij)  *  on2(uj)  6  (d)  (B-22) 


9  +6.  <p  +<f 

<"  o  1  r  o 


f  7 


•  a 

j  —  cos  9  „ 

Sj  ( 0,  $  ,w)e  c  0  d<J>d0 


v  J0  -<t 

o  1  yo 


=  on  (w) 6 (d) 


(B-  23) 


f  01  A+*  c  i  it  GOS  <e+e  )  2 

+  f  /  Sj(  6+0q,  (t+^o,w)  e  q>  d$d0 

J-\  J0Q-0 


Proceeding  as  in  Case  B.l,  assume  that  is  small,  so 


cos  ( 0+9  )  s%  cos  0„  -  0  sin  0 
o  o  o 


(B-24) 
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and  assume  that  Sj  (0,$,^)  is  separable. 


SI( 


>  =  SIh  ( 


Slv( 


(B-25) 


so 


se(d.  v)  cn2(  .)  5  (d) 


+ 


sIv(  4> 


e 


SIh(e’") 


■i# 


1  sm  d 


de 


(B-26) 


When  the  source  CDS  is  assumed  uniform  in  S,  this  gives 


sg(d,  j) 


<*<j  2(jj)  6  (d) 
n 


(B-27) 


+  a , 


(!*>)  G  (ui) 
<P 


sin  [  2^0^  y  sin 


V 


d_ 

it  A  sin 


•  bud 

l  — cos 
J  c 


where 


<VW) 


Sj(  <p,ui) 


d0 


(B-28) 
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This  has  exactly  the  same  form  as  (A-41)  and  as  the  CSD  of  the  far  field 
model  with  a  uniform,  uncorrelated  source.  As  before,  for  a  uniformly 
spaced  line  array  along  the  y  axis  with  element  spacing  d.  the  CSD  between 
the  outputs  of  the  pt^1  and  q**1  hydrophones  is 


se(p,q..)*on  U)  cpq 


(B-29) 


+or,  Sin  (Hp'q)6lTSin  j 

I(co)  dr  ' - -= - e  J 


x)  d  /  v  ^ 

—  (p-q)  cos  eQ 


T(P-Q)  -  sin  eQ 


B .  3  Summary.  Two  source  models  that  can  be  used  as  approximate 
representations  of  surface  reverberation  from  the  convergence  zone  have 
been  developed.  Under  very  reasonable  assumptions,  it  has  been  shown  that 
these  models  produce  a  hydrophone  output  CSD  of  the  same  form  as  an 
uncorrelated  source,  which  was  considered  in  The  assumptions 

needed  for  use  of  the  model  are 

(a)  The  wavefront  emanating  from  each  point  on  the  reverberating 
surface  can  be  regarded  as  planar  across  the  array. 

(b)  The  reverberation  is  confined  to  an  azimuth  angular  sector  of 

width  20^,  with  9^  small  enough  that  sin  0^  P^  and  cos  0^  1. 

The  transmitting  sonar  must  therefore  provide  sufficient 
directivity  to  justify  this  assumption  if  the  model  is  to  be 
valid . 

(c)  The  reverberation  at  each  point  on  the  surface  is  uncorrelated 
with  that  at  any  other. 


B-12 


(d)  The  spectral  density  of  the  reverberation  over  the  two 


V 

r 


I 


dimensions  is  separable. 

(e)  The  spectral  density  of  the  reverberation  is  uniform  in  azimuth. 
In  addition,  the  vertical  arrival  angle  model  requires  that 

(f)  The  vertical  arrival  angles  be  concentrated  in  a  narrow  sector 
about  the  horizontal,  say  -<t>  <<j>  <  $  ,  such  that  sin  $  s% 

cos  <f  »  1. 

Under  these  assumptions,  the  results  for  an  uncorrelated  line  source 
can  be  directly  applied  to  the  study  of  the  cancellation  of  convergence  zone 
surface  reverberation.  These  results  also  suggest  that  the  line  source 
model  can  be  used  to  approximate  other  two  dimensional  sources  of  narrow 
angular  extent. 
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APPENDIX  C:  DERIVATION  OF  THE  OPTIMAL  CANCELLER 
AND  CANCELLER  OUTPUT  SPECTRUM 


APPENDIX  C;  DERIVATION  OF  THE  OPTIMAL  CANCELLER 
AND  CANCELLER  OUTPUT  SPECTRUM 

Consider  the  multiple  reference  canceller  configuration  shown  in  Figure 
C-l,  in  which  the  output  of  K  reference  hydrophones,  denoted  e.  (t)  with 
i  =  0,1  ...  K-l  are  to  be  used  to  cancel  an  interference  from  the  output  of  a 
primary  hydrophone,  ed(t).  In  keeping  with  the  discussion  of  Appendix  A, 
the  hydrophones  outputs  are  assumed  to  be  zero  mean  random  processes,  that 
are  at  least  wide  sense  stationary.  That  is, 


E[ep(t)  eq(t+T)]  =  r(p,q,x) 


p  =  0.1 _ K-l 

q  =  0,1, ...K-l 


where  E [  - ]  denotes  expectation.  The  filters  are  linear  and  time  invariant 

with  impulse  responses  h.(t)  for  i=0,...K-l.  The  impulse  responses  are  to  be 

2 

chosen  to  minimize  the  mean  square  value  of  the  error  output,  E[e  (t)] .  The 
derivation  here  closely  follows  that  of  reference  [14] 

From  Figure  C-l,  it  can  be  seen  that 


K-l 

£(t)  =  e^tt)  -  7^  j  h^fx)  ek(t-T)  dT 
k=0 


(C-l) 


Figure  C-l.  Multiple  Reference  Canceller 


so  that 


E[ e  (t)]  =  ru(0)  -  2 


K-l 

El 

k=0 


hk(r)  r1(k,  t)  dr 


(C-2j 


K-l  K-l 


E  E  L  L  \ 

r  r  a  1»  — A  A 


K^O  k2=0 


(x)  hk  (p)  re(k1>k2,i-p)  drdp 


with  r  (p,q,i)  the  cross-correlation  between  e  (t)  and  e  (t),  r,(p,T)  the 
e  p  q  l 

cross-correlation  between  e^t)  and  e^(t),  and  r^(i)  the  covariance  of 

ed(t). 

Let  the  optimal  impulse  responses  be  denoted  h°(t)  for  i  =  0,1,. .  .K-l, 
and  let 


hk(t)  =  hk(t)  +  7k  5hk(t) 


CC-3) 


where  6hk(t)  is  an  arbitrary  variation  about  hi^(t )  and  Yk  a  real  constant. 
Substituting  (C-3)  in  (C-2)  gives  the  mean  square  error  as 


E[E*(t)]  =  JQ  +  2 


k^O 


K-l  K_1  ff 

£  Yk1  2  /rhk(T)  hk2(p)  rl(k’T-p)  dTdP 


k2=° 


f6\ 


(t)  fgO^i’k^x)  dx 


K-l  K-l 

E  E 

kfO  k2=0 


ff  <5hk  (x)  5hk  (p)  r^k^kg.T-p)  didf 


Yk  Yk 
K1  K2 


(C-4) 
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i 


where 


J 

o 


ru(0)  -  2 


K-l 

2  / 

1 ,  A 


hkl(T^  ri(ki>T)  dT 


kx=0 


+ 


K-l  K-l 

2  S  /T  h°  (t)  h°  (p)  re(k1,k2>x-D)  dxdp 
k  =0  k  =0 


is  just  the  minimum  mean  square  error.  All  integrals  have  upper  and 
limits  of  °°  and  -°°,  respectively. 

If  a  minimum  exists,  it  will  occur  when 


hk(t)  =  hk(t)’  k  =  °’  l'  K_1 
A  necessary  condition  for  this  minimum  to  occur  is  that 


E[e2(t)] 


=  0,  k  =  0,  1, 


K-l 


(C-5) 


lower 


(C-6) 


(C-7) 
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Using:  (3-4) , 


3Y. 


E[e  (t)] 


=  2 


Z//ih„'-'  h° 

i=0  ^  J 


(p)  re(n,i,x-p; 


Y0  =  0 


K-l  =  0 


/6hn(T)  rl 


(n ,  t)  dx 


For  this  to  equal  zero  for  arbitrary  variation  of  6hn(t), 


K-l 


r1(n,x)  = 


■  Z  f< 

k=0  J 


k(p)  rg(n ,k, i-p)  dp 


This  must  be  satisfied  for  n  =  0,  1,  ...  K-l.  Recall  from  Section  2.1 


Sn(w)  =jr[ru(T)] 


Sjfp.w)  =  J^[r^(p ,  t)  ] 
Sg(p  .q .  w)  [re(p  ,q,x)] 


where  .&[•]  denotes  the  Fourier  transform  with  respect  to  x,  given  by 
Then  taking  the  Fourier  transform  of  (C-9)  gives 


Sjin.w)  = 


K-l 


k=0 


Hj(u)  Se(n,k,  «)  n  =  0.  1, 


,  K-l 


dxdp 


(C-8) 


(C-9) 


that 


(A- 12) . 


(C-10) 
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where  H°(oj)  is  the  transfer  function  of  the  kth  optimal  filter.  Now,  as 
in  Appendix  A,  let 


S^toj) 


Sj^O.w),  s1(l,w),  ...  s^(K-1,oj) 


T 


(C-U) 


and 


S  (oj)  =  (s  (p,q,o)))  ( C—  1 2 ) 

6  P=0, 1, . .  .K-l 

q-o,i,...K-i 


Then  (C-10)  can  be  written  as 


s^w)  =  Se(co)  il  (oj) 


( C—  1 3) 


where  H°(w)  is  a  vector  of  the  optimal  transfer  functions, 


pj 

H° (to)  =  ^Hj(w),  H°(w),  ....  H° ^(w)] 


(C- 14) 


Therefore, 

H°(uj)  =  S'^oj)  s^w) 


( C— 15 ) 


Fourier  transforming  (C-3)  to  yield  the  spectrum  of  the  irreducible 


error,  termed  the  irreducible  canceller  output  spectrum  here,  gives 


Su(w)  -  2 


Ek(co) 


K-l 

y:  s x ( k , u.' ) 

k=0 


+ 


E 

kr< 


K-l 

)  H°  Cuj)  h£*(u)  se(k1.k2.u) 
k2=0  1  2 


=  S11(oj)  -  2H°  (w)  s^w)  +  H°  (u)) 


S  (w)  H"  (u) 
e  — 


(C-16) 


Substituting  the  optimal  response  from  (C-13)  gives  the  canceller  output 
spectrum  as 


Ek(oj  ) 


s^(oj)  S{ 


(w) 


(C-17) 


Here,  the  subscript  K  in  E^(oj)  indicates  that  K  references  have  been  used  in 
the  canceller. 

It  will  be  useful  to  express  E^(w)  as  the  ratio  of  determinants, 


det[SH(w)] 
EK(w)  =  det[Se(cu)3 


(C-18) 


where  S„(cj)  is  a  (K+l)  x  (K+l)  matrix  defined  as 

H 


(C-19) 


This  is  easily  shown  using  the  following  matrix  identity. 


Sjj(od)  = 


Se(^ 

S^foo) 

S*(w) 

S11(w) 

C-6 


Theorem :  Determinant  of  a  Bordered  Matrix 


Let  the  (N+l)  x  (N+l)  square  matrix,  B,  be  defined  as 


B 


A  x 
+ 

I  c 


(C-20) 


where  A  is  an  arbitrary  N  x  N  matrix  x  is  a  N-dimensional  column  vector, 
y+  an  N-dimensional  row  vector,  and  c  a  scalar.  Then 

dettB]  =  c(det  [  A] )  -y+(adj[A])x  (C-21) 

where  det  [  •  ]  and  adj  [  •  ]  are  the  matrix  determinant  and  adjoint, 
respectively . 

Proof :  Let 


n 


*  =  [yi  •  •  •  yn  ] 


Using  the  Laplace  expansion  formula  with  the  elements  of  the  N+l  column 
of  B  as  coefficients  gives 


N+l 


i+N+1 


detlBl  =  2  bi,N+l  (- 


1) 


6i,N+l 


(C-22) 


i=l 


with  b.  M  the  ith  element  of  the  (N+l)th  column  and  6.  the  minor  of 
i,N+l  i,N+l 

b.  KTj_, .  Note  that 
l,N+l 


C-7 


bN+l,N+l  C 


and 


JN+1,N+1  =  (A) 


(C-23) 


Further, 


b.  KTj,  =  X.,  1  <  i  <  N 
i,N+l  1  -  - 


(C-24) 


Now ,  the  minor  of  b.  V7 , ,  for  1  <  i  <  N  is  the  determinant  of  the  matrix 
l,N+l  -  - 

obtained  by  deleting  the  i-th  row  and  the  (N+l)th  column,  denoted  C..  That 
is. 


C.  = 

l 


an  •  •  •  ■ 

'  al,N 

•i-i.r  •  ' 

■  ai- 1 ,  N 

ai+i,r  •  ' 

ai+l ,  N 

aN,i  •  • 

;  aN,N 

* 

Jk 

y  i  •  • 

•yN 

(C-25) 


.  .th 


where  a.  .  is  the  i , j  element  of  A.  To  obtain  8.  M  -  ,  reapply  the  Laplace 

1,]  1,N+1 


expansion  formula  with  the  y. 's  as  coefficients 

Si.NH  =  de«Ci)  "  E  y‘  (-I>  “i.i 


(C-26) 


with  a.  .,  determinant  of  the  matrix  remaining  when  the  Ntb  row  and  jtb  column 
t>] 

of  C.  are  deleted.  Inspection  of  (C-25),  however,  shows  that  a.  .  is  exactly 
x  i  >] 

the  minor  of  a.  .  in  the  matrix  A. 
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C-8 


Therefore , 


N 


N 


det [ B ]  =  e(det [ A] )  7  ,X.(-1) 

i=l 


Exi< 


i+N  +  l 


3=1 


N  N 

c(det [  A] )  EE 

i=l  3=1 


y.  (-l)1+]  a.  .  X. 
3  i.3  i 


The  quantity  a.  .is  the  cofactor  of  a.  .  in  the  matrix  A 

1.3  i.3 

Hence 


det[B]  =  c(det[A])  -  y 


A  A  A 

A11  12  •  •  •  a1N 


AN1  AN2  ‘  •  •  ANN 


or 

det [ B ]  =  c(det[A] )  -  }[*(Adj[A])X 
Using  this  theorem, 

det [ S ( oo) I  =  Su(oj)  det[Se(oo)]  -  s^(w)  Adj [ Se ( oo) ]  s ^ ( ca) 

where  Adj  [A]  in  the  adjoint  matrix  of  A.  Clearly  then, 
det[SH(uj)J  + 

det[Se(w)l  =  sn(a))  “  -i(w)  se  -i  (t0^  =  ek 
In  the  special  notation  of  (A- 22)  and  (A- 2 3) 


>j 


(C-  27) 


,  denoted  A.  .. 

i.) 


(C-28) 


(C-29) 


(C-30) 


(C-  31) 


C-9 


EK(jJ)  =  S11(jj)— l+(K,;i'+)Se1(K,u)  -1(K,U;) 


(C-32) 


det[SH(K,x')] 

det[Se(K,u))J 


where 


Sh(K,oj)  =  Sh(oj) 


'K  references 


(C-  33) 
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APPENDIX  D:  THE  EQUIVALENCE  OF  EXTENDED  SOURCE 
CANCELLATION  AND  THE  LINEAR  PREDICTION  OF  A 
BAND  LIMITED  TEMPORAL  SEQUENCE 


APPENDIX  D;  THE  EQUIVALENCE  OF  EXTENDED  SOURCE 
CANCELLATION  AND  THE  LINEAR  PREDICTION  OF  A 
BANDLIM1TED  TEMPORAL  SEQUENCE 

Consider  K+l  hydrophones  oriented  as  shown  in  Figure  D-l,  and  let 
the  output  of  the  pth  hydrophone  be  (t),  where  the  primary  hydrophone 
from  which  the  interference  is  to  be  cancelled,  takes  the  (J+l)1  position.  Let 
the  cross-correlation  function  between  the  pth  and  qth  hydrophone  outputs 
be  given  by 

R  (uj)  =  E[n  (t)  n  (t+T)]  (D-l) 

pq  P  q 

and  let  the  cross- spectral  density  between  the  pth  and  qth  outputs  be 

s  (p.q.co)  =,^[R  (t)]  CD-2) 

e  pq 

where,^[  ]  denotes  the  Fourier  transform  with  respect  to  the  variable  t. 

Now ,  denote 

Se(w)  =  (se(p,q,u»  p  =  0,  1....K-1 

q  =  0,  1....K-1  (D-3) 

SjCco)  =  [Si(0,id),Si(1,u>),  ....  s1(K-l,a))]T  (D-4) 

It  is  assumed  that  the  process  is  spatially  as  well  as  temporally  stationary, 
so  that 

Se(p.p+A  ,w)  =  se(q,q+A  ,oj)  (D-5) 


D-l 


n 


■(J-l)d 


-(K-l)d- 


|*-  d  •*]*-  d  4-d^) 


-ff- 


REFERENCE 

HYDROPHONES 


PRIMARY 

HYDROPHONE 


Figure  D-X :  Hydrophone  Geometry  for  Spatial  Cancellation 


/ 
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f 


for  any  p  and  q.  It  is  shown  in  Appendix  C  that  the  power  spectral  density 
of  the  canceller  output  for  K  references  is 


Ek(io)  =  SuU)  -  s1+(uj)  Sg  1(> 


(D-6) 


Next,  consider  the  temporal  problem  of  predicting  a  stationary  random 

sequence,  m  steps  ahead.  Let  a  continuous  random  process,  x(t),  with 

correlation  function  R  (t)  be  sampled  at  a  uniform  rate  of  l/T  samples  per 

X  s 

second,  so  that  the  random  sequence  is  jx(nT  )f  ,  „  ,  „  .  The 

I  s  )  n=.. ,-1,0, 1,2, . . . 

correlation  between  the  p  n  and  q  n  samples  will  be  R  [(p-q)T  ]  .  Given 

X  s 

the  finite  past  {x(-KT  ) , . . .  ,x(-T  )  },  the  goal  is  to  predict  in  the  minimum 

s  s 

th 

mean  square  error  sense  the  value  of  the  mu  1  element  of  the  sequence, 
x(mTs>.  It  can  easily  be  shown  that  the  prediction  error  is 


E(m)  =  Rx(0)  - 


r  ,  T  R 
-dx  xx 


-1 


-dx 


(D-7) 


where 


Rxx  =  (Rx[(p-q)Ts]).  P  =  -K,-K+l . -1 

q  =  -K  ,-K+l . -1  (D-8) 

=  (Rx[(p-q)Tsl),  p  =  0,1, . .  .K-l 
q  =  0,1,. .  .K-l 

and 

-dx  =  [Rx(m)’  Rx(m+1)’  •••  Rx(m+K)]T  (D-9) 
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Comparing  (D-6)  and  (D-7),  it  can  be  seen  that  the  temporal  linear 

prediction  problem  is  equivalent  to  the  spatial  LMS  cancellation  at  a  single 

frequency.  The  temporal  statistics,  R  ,  r  .  ,  and  R  (o)  correspond  to 

xx  — ax  x 

S  (u),  s^(a)),  and  S^(gj)  in  the  spatial  case.  It  is  convenient  to  note  that 
in  the  spatial  case,  using  the  notation  of  Appendix  A, 

S^Cp.u.')  =  se(J-l,p,w)  (D-10) 

and 

S11(w)  =  se(p,p,cj),  any  p  (D-ll) 

Therefore, 

s^uO  =  [se(J-l,0,oj),  s  (J-l.l.u),  ....  se(J-l,K-l,w)JT  (D-12) 

This  makes  the  equivalence  of  Rx[(p-q)Tg]  and  Se(p,q,w)  clear. 

The  analogy  can  be  extended  one  step  further  if  the  source  is  assumed 
to  be  a  narrow  plane  wave  source  in  the  same  plane  as  the  array ,  distributed 
in  angle  <t>,  as  shown  in  Figure  D-2.  Let  the  arrival  from  any  angle,  <J> ^ ,  be 
uncorrelated  with  arrivals  for  any  angle,  $2»  if  <j>^  4-  and  let  the  source 
be  distributed  in  angle  from  4>Q  -  <J>  to  <J>Q  +  $  with  angular  density  Sj  (ui,<J>) 
in  that  interval.  The  source  is  assumed  to  be  narrow  if 

sin  $  =  $ 

cos  <5=1  (D-10) 


D-4 


2 

Also  assume  that  ambient  noise  with  spectral  density,  (oo)  appears  at  the 
hydrophone  outputs,  and  is  uncorrelated  from  hydrophone  to  hydrophone. 
Then,  from  Appendix  A, 


se(p,q,w)«:  on  (w)  5pq 


(D-ll) 


+  e 


(p-q)cos  <j)0 


/\ 


+]'7r  (p-q)$sin$ft 

(0+$Q,aj)e  ud4> 


This  discussion  will  be  facilitated  by  demonstrating  that  the  value  of  the 
canceller  output  spectrum,  Ek(ij),  given  by  (D-6),  is  not  affected  by  the 
factor  exp  [-j^-  (P-q)  cos  0Q]  in  the  second  term  of  (D-ll).  From 
equation  (C-27) 


det  [S  (K,w)] 
EK(w)  =  det  [Se(K ,oj)] 


(D-12) 


where 


(D-13) 


Given  the  hydrophone  CSD  of  the  form  of  (D-ll),  it  is  possible  to  write 


SH(K  ,ui)  = 


Se(K,w)  s^fK.w) 


Sj  (K , oj)  S^Cw) 


SH  (K ,oj)  =  G*  (K)  SHo(K,w)  Gl  (K) 


(D-14) 


Sa  (K.w)  =  G+  (K)  S  (K.oj)  G(K) 
“  60 


(D- 15) 


where 


G(K)  =  diag 

p=0,  1 . K-l 


ojd 

+1P  —  cos  4>q 


(D- 16) 


G^K)  = 


G(K) 


-j  K—  cos  <p 
1  c  o 


with  0  a  K-dimensional  vector  of  zeros  and 


Seo  (K’U)  =  (seo  (p’q’w))p=0,l, 

q=0 , 1 , 


..K 

..K 


(D-17) 


with 


Seo(p’q’w)  =  °n  (w)  5pq  + 


f  $  2 

■  j  Oj  (U), 


-j  ^(p-q)  ~r-  sin  <Pn 
(J))e  d<p  (D-18) 


also 


SHo(K’w)= 


SeQ  (K  ,  a))  sQ  (K  ,u>) 


sQ  (K,w)  Sn(K) 


(D- 19) 


where 


3  (K ,w)  =  [s  (o,m,0))t  s  (1,  m,0)) . s  (K-l,  m,u>)]  (D-20) 

vi  C  w  Cv  «  (J 


D-6 


It  can  be  seen  that 


det  [  G  ( K )  ]  = 


K-l 


wd 

+  ]P  -  COS 

C 


0  =  (det  [G+(K)] )  * 


(D-21) 


and 


...cod 

-  jK —  cos  <j> 

det  IG1(K)]  =  e  c  0  det  [G(K)]  =  (det  [G1  (K)])  * 


(D-22) 


Then  using  the  fact  that 


det  [ABC]  =  det  [A]  det  [B]  det  [C] 


(D-12) ,  (D-.14) ,  (D-15) ,  (D-I7) ,  and  (D-18)  give 


det  [S  (K , w) ] 
EK(w)  =  det  [SeQ  (K  , oj) ] 


(D-23) 


Therefore,  (D-ll)  can  be  replaced  without  loss  of  generality  by 


se(p,q,w)  *  on  (uj)6 


(D-24) 


r  $  (p-q) 0  sin  4) 

+  /  S  j  ( 4>+<i>g » u)  )e  d  4> 


so  that  the  CSD,  s^p.q.oj),  is  the  Fourier  transform  of  the  source  density, 

Sj(  <!>+<}>  Q,oj) . 


D-7 


In  the  temporal  case,  when  uncorrelated  noise  with  variance 
added  to  each  time  sample , 


RXUP 


q>Ts]  = 


pq 


L 


S  (f)e 
xx 


.2irf(p-q)Ts 


df,  f  = 


0) 

2tt 


(D-25) 


where  Sxx(f)  is  the  power  spectral  density  of  the  process  x(t).  If  the 
process  is  bandlimited, 


S  (f)  = 
xx 


Sxx<f+V’  VBsfsf0+B 

0  ,  otherwise 


(D-26) 


Then  (E-14)  becomes 


Rx(P'q)Ts] 


=  a 


n 


pq 


/■ 

-B 


Sxx<f+f0)e 


j2irf  (p-q)T 


df 


(D-27) 


Comparing  (D-24)  and  (D-27),  it  can  be  seen  that  the  angular  density  of 
the  extended  source  is  analogous  to  the  spectral  density  in  the  temporal 
case.  Further  the  angular  extent  of  the  source,  corresponds  to  the 
bandwidth  of  the  temporal  process,  B,  and  the  temporal  sample  interval, 
Tg,  corresponds  to  the  spatial  sample  interval 


sin 


with  X  the  wavelength  of  the  frequency,  f. 


Now,  the  samples  of  the  random  process,  x(t),  comprise  a  random 


sequence  x  ,  n  ,  0  with 

M  n  n=. .,-1,0, 1,2,... 


xn  =x(nTs) 


Define 


Rk  ~  E  [xn  xn+k] 


and  note  that 


®k  =  Rx  tkTs] 


The  power  spectral  density  of  the  sampled  sequence,  jxnj  >  maY  be 
on  [-n,u)  as 


.ka 


>  =  2-,  VJ 

k=-°° 


where  the  inverse  transform  relationship  is 


Rk  = 


r  77 

Tn  [  H 

“  -  TT 


_.ka 

,(a)e  J  da 


Note  that  (D-27)  can  be  rewritten  as 


R  [  (p-q)T  ]  =  o  6  + 

x  v  M  s  n  pq 


.  2ttBT 

/  sxx  {wr +  fo)  ®3 


j(p-q)a  da 
2ttT 


(D-28) 

(D-29) 

(D-30) 

defined 

(D-31) 

(D-32) 

-  (D-33) 
s 


D-9 


If  it  is  assumed  that  BT  <  1/2,  then  direct  comparison  of  (D-32)  and  (D-33) 

s 

shows  that 


H 


x 


(a)  =  °n 


2 


+ 


Rect 


a 

2ttBT 

s 


(D-34) 


On  the  other  hand,  if  BTs  >  1/2,  the  aliases  of  the  sampled  power  spectral 
density  must  be  considered,  and  Hx(a)  must  be  replaced  by  H^x(a)  in 
(D-32)  where 


Htx(a) 


-°n2+  H 


k=-° 


(D-35) 


/a+2irk  f  \ 
xx  \2ttT  *0/ 


Rect 


k=-° 


a+2iTk 

2ttBT 


,  ae[-Tr,TT) 


Now,  consider  the  spatial  cancellation  problem  starting  with  equation 
(D-24).  Define  the  elements  of  the  hydrophone  cross- spectral  density 
matrix  as 

Sk  =  *\E[en  I1  len+k  (t+T)1|  (D_36) 
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Analogous  to  the  power  spectral  density  of  the  temporal  process  is  the  wave 
number-frequency  spectrum  (WNFS)  of  the  hydrophone  outputs,  given  by 


He(ct.u))  = 


Sk  e 


jka 


OC  £  [  —  TT  ,  7T  ) 


(D-  37) 


k=-° 


with  the  inverse  relationship 


Sk  = 


if 

•'-7T 


He(a,w)  e 


-jka  da 


(D-38) 


However,  (D-24)  gives  as 


Sk  =  °n  (w)  6  (k)  + 


L 


-]4>k^-  sin  0 

Sj  (w,  0+0Q)e  c  dd> 


(D-  39) 


Making  the  substitution 


.  ud  . 

a  =  4)  —  sin  <t>Q 


da  =  —  sin  tp  d<J> 
c  o  r 


(D-40) 


and  letting 


1  cod  >  ,  d  . 

Y  =  U  ~sm  (po=XSin  % 


D-ll 


gives 


(u>)  6  (k)  + 


/. 


2TTY0 


Sj  ( u3 , 


2  TTY 


V 


-jka  dot 


2tty 


(D-41) 


-2tty  $ 


Assume  for  the  moment  that  y^-1/2.  Then  comparing  (D-41)  with 
(D-36)  shows  that 

n  S  L,  J*-  +  <£  J 

He(a,a))  =  (oj)  +  - ill - Rect  [a/2Try4>]  (D-42) 


where 

•x]  _  j  1,  |x  I  s  X 

XJ  jO,  otherwise 

Next,  consider  the  case  when  y$  >  1/2,  so  that  the  noise  field  is  under 
sampled  spatially,  and  the  portion  of  the  spectrum  extending  beyond  [-tt.tt) 
will  fold  back  into  [  - tt,  tt).  Let  be  the  folded  spectrum 


Rect 


(w)  +  ^ 

k=-oo 


Hg(a  +  2-rrk  ,(*>) ,  ae  [-tt,  tt) 


(D-  43) 


where  He(a,u))  is  defined  by  (D-42). 
Then 


u> 


Hfe(a,w)  =  an  (to) 


a+2irk 

2tty 


+*o) 


Rect 


k=-» 


a+2irk 


2tty$ 


ae[-TT,  tt) 


(D-44) 
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D.l.  Results  from  linear  prediction  theory.  The  applicable  results  from 
linear  prediction  theory  fall  into  two  general  categories.  In  the  special 
case  of  m  =  1,  when  the  line  array  of  references  and  the  primary  constitute 
a  uniformly  spaced  line  array,  results  depend  upon  the  Toeplitz  properties 
of  the  hydrophone  cross- spectral  density  matrix.  The  second  class  of 
results  which  apply  to  m  >  1  are  from  the  more  general  theory  of  stationary 
random  functions.  Most  of  the  results  to  be  presented  concern  the  irre¬ 
ducible  prediction  error,  which,  for  the  temporal  prediction  problem,  is  the 
mean  square  error  in  predicting  the  desired  number  of  steps  ahead  using 
the  infinite  past.  If  E,.(m)  denotes  the  error  in  predicting  m  steps  ahead 
using  K  samples  of  the  past,  then 


ER(m)  =  min 


aK-l,aK-2,...,al’ao 


O 

x(m)  -  ^  aRX(k) 
k=-K+l 


(D-45) 


where  E(-)  denotes  expectation.  Then  the  irreducible  prediction  error  is 
E^Ui^lim  E„  (m )  (D-46) 

K-*oo  * 

In  the  context  of  the  cancellation  of  the  extended  source,  the  results 
apply  at  a  single  frequency,  so  the  notation  must  be  changed  to  E  (m,u)) 
and  E^Im.cj).  The  parameter  ER(m,u))  is  the  spectral  density  of  the  can¬ 
celler  output  at  the  frequency  u)  using  K  reference  hydrophones  and  with 
the  primary  md  feet  from  the  nearest  reference.  The  E^Im.w)  is  the  limit 
on  the  single  frequency  canceller  output  spectrum  as  the  array  becomes 
infinitely  long  (K-»-°°).  This  clearly  conflicts  with  the  assumption  of  a  plane 
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wave  source,  which  requires  that  wavefront  curvature  be  negligible  across 

the  extent  of  the  array.  The  irreducible  error,  E^Cm.-u),  is  therefore  best 

interpreted  as  a  lower  bound  on  the  canceller  output  spectrum  in  the  spatial 

case.  It  should  be  noted  that  in  numerical  evaluations,  given  in  Appendix  K, 

the  lower  bound  was  approached  quite  rapidly  as  K  increased,  so  that  K 

may  not  have  to  be  very  large  for  El/(m,cu)  to  be  near  E  (m,u)). 

A  basic  result  from  the  theory  of  Toeplitz  operators  in  the  Szego 

* 

r  i5| 

Theorem  1  ,  which  relates  the  irreducible  prediction  error,  E^Cl).  to  the 

spectral  density  of  the  random  sequence,  Hx(uj). 

Theorem  1:  Let  a  zero  mean  stationary  random  sequence,  jx(n)J  ,  have 
correlation  function 


Rk  =  E[x(n)x(n+k) I 


(D-47) 


and  power  spectral  density 


Hx(a)  =  2-r  V 


(D-48) 


Then 


E  ( 1)  =  lim  E  ( 1)  =  exp 


r 71 

2 -  /  tosH> 

-Tt 


(a)da 


(D-49) 


When  log  [Hx(a>]  is  not  inte  grable  on  the  right  hand  side  of  (D-49) 

is  replaced  by  zero.  The  theorem  as  given  here  is  somewhat  more  restrictive 


than  given  in  [15]  in  that  it  assumes  the  existence  of  the  spectral  density, 
rather  than  allowing  for  a  spectral  distribution  that  contains  jumps 
or  that  is  not  everywhere  differentiable. 

Note  that  under  the  last  stipulation  of  the  Theorem,  any  time  log  [Hx(a)j 
is  not  integrable,  the  irreducible  prediction  error  is  zero,  so  the  process 
can  be  predicted  perfectly  from  the  infinite  past.  Such  a  process  is  termed 
a  singular  stochastic  process  by  Lamperti  [16]  or  a  deterministic  stochastic 
process  by  Azego  [15].  The  condition  under  which  the  process  is  singular 
which  will  be  of  interest  here  is  as  follows. 

Theorem  2:  A  stationary  random  process  is  singular ,  i.e.,  it  can  be 
predicted  perfectly  in  the  mean  square  sense  from  its  infinite  past,  if 


log  Hx(a)  da  =  -» 


(D-50) 


This  will  occur  if  Hx(a)  is  zero  for  any  cce [ - rr , rr) .  Then  any  random 

sequence  whose  spectrum  is  zero  in  some  region  on  [ - rr , -t)  will  be  singular. 

2 

Note  that  if  an  ^  0  in  (D-34),  then  the  process  is  not  singular,  as  would  be 

2 

expected.  In  fact,  the  irreducible  error  will  always  be  greater  than  an  , 

2 

except  in  the  trivial  case  of  H  (a)  =  a 

x  n 

Theorem  3:  Let  the  stationary  random  process  x(n),  have  spectral  density 
of  the  form 


Hx(o)  =  on2*Hx0(“>'H*0<c‘>*° 


CD-51) 
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Then  E  (1)  •  o  .  This  can  be  seen  by  observing  that 
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since  the  argument  of  the  integral  is  always  positive  if  HxQ(a)  4  0. 

In  terms  of  the  cancellation  of  an  extended  source,  Theorems  1  through 

3  apply  to  the  power  spectrum  of  the  canceller  output,  E^d.u).  As  shown 

above,  the  wave  number- frequency  spectrum,  He(a,o)),  is  directly  related 

to  the  angular  distribution  of  a  narrow  extended  plane  wave  source.  These 

theorems  therefore  provide  conditions  on  the  source  distribution  under 

2 

which  complete  cancellation  is  possible  in  the  noise  free  (o  (w)  =  0)  case 
and  indicate  that  cancellation  to  the  noise  floor  is  not  possible  in  the 
presence  of  ambient  noise. 

Theorem  la:  Let  ^be  an  extended,  uncorrelated,  plane  wave  source  as 


defined  above ,  distributed  in  the  same  plane  as  the  reference  array  and 


r 

i 


primary  of  the  geometry  of  Figure  D-l.  In  the  case  whem  m  =  1,  the  wave 
number-frequency  spectrum  of  the  hydrophone  outputs  is 
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Under  the  condition  that  y$  s  1/2,  this  reduces  to 


«.«■•“>  -  -  7  SI  (liy  *  RMt  [n7T] 


(D-54) 
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Theorem  2a:  Let  (oj)  =  0  and  assume  that  y$  s  ?>•  Then  the  extended 
source  gives 


E.d.u)  =  0  (D-55) 

2 

Theorem  3a:  Assume  that  on  (or)  #0.  Then  by  Theorem  3,  the  extended 
source  cannot  be  cancelled  to  the  noise  floor. 

The  results  of  Theorem  2  can  easily  be  extended  to  m  >  1  by  observing 
the  following. 

Theorem  4:  Let  Eot  (m.w)  be  the  irreducible  error  in  predicting,  in  the  minimum 
mean  square  sense,  the  stationary  random  sequence,  x(n),  m  steps  ahead. 

If  Ea(  1)  =  0,  then 

E0(m)  =  0  m  =  1,2, -  (D-56) 
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This  can  be  proven  by  noting  that,  under  the  stationarity  assumptions, 
x(2)  can  be  predicted  perfectly  using  the  infinite  past,  J  ....  x(0),  x(l)j 
ifEoo(l)  =  0.  Given  j  ....  x(-l),  x(0)J  ,  it  is  possible  to  predict  x(  1) 
perfectly.  Hence,  x(2)  can  be  predicted  using  j  ....  x(-l),  x  ( 0)  J  with 
zero  mean  square  error,  £^(2)  =  0.  By  induction,  (D-56)  is  true. 

This  result  means  that  in  the  noise  free  case,  the  canceller  output 
power  for  an  extended  source  tends  asymptotically  (in  K)  to  zero,  regard¬ 
less  of  m,  that  is,  regardless  of  the  distance  between  the  primary  hydro¬ 
phone  and  the  reference  array. 

In  the  case  when  E^d)  =  0,  it  is  often  possible  to  determine  the 

asymptotic  behavior  of  E  (1)  for  large  K  using  the  Strong  Szego  Limit 

K 

r  19] 

Theorem  for  Circular  Arcs1  .  Applied  specifically  to  the  problem  being 
considered  here,  the  Theorem  is  as  follows: 

Theorem  5:  Let  x(n)  be  a  stationary  random  sequence  with  spectral 
density  H  (a)  and  correlation  R,  ,  k  =  ...-1,0,1,....  Let  H  (a)  be 
strictly  positive  on  -A  <  a  <  A  with  A  <  ir,  and  zero  otherwise.  Further 
let  dHx(a)/da  satisfy  a  Lipschitz  condition  with  exponent  greater  than  1/2 
on  -A  <  a  <  A.  Then  the  minimum  mean  square  error  in  predicting  x(n) 
given  |  x(n-K),  x(n-K+l),  ....  x(n-l)  j  using  a  linear  combination  of  the 
observations  is  asymptotic  in  K  to  L(K),  where 


L(K)  =  G 


2K+1 


(D-57) 
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and 
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The  theorem  requires  a  modification  of  the  version  presented  in  [  12] .  The 
modification  is  derived  in  the  Supplement  I  to  this  Appendix. 

Theorem  5  can  be  applied  to  the  spatial  cancellation  problem  by  replacing 
H  (a)  with  H  (a,u))  as  given  in  Theorem  la.  Inspection  of  equation  (D-42) 
shows  that,  due  to  the  aliasing  of  the  angular  distribution  of  the  source  when 
Y$  >  1/2,  the  Theorem  only  applies  when  yi>  <  1/2. 

Theorem  5a:  Let  y^  <  1/2,  and  let 


He(a,oj) 


1 

Y 


a 

2  tty 


+  $ 


0’ 


“) 


Rect 


fa] 


(D-59) 


where  H  (a.aj)  is  strictly  positive  on  -A  i  a  <  A,  and  where  9H  (a,cu)/3a 
e  e 

satisfies  a  Lipschitz  condition  in  a  with  exponent  greater  than  1/2  on  that 
interval.  Then  the  irreducible  canceller  output  spectrum,  E„(l,w),  of  an  LMS 

xV. 

canceller  using  a  line  array  of  K  references  spaced  d  feet  apart  to  cancel  the 
source  from  the  output  of  a  primary  in  the  same  line  as  the  references  and 
d  feet  away  satisfies 


E^(  1,cj)  'v  G(w) 
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where  ^  denotes  asymptotic  in  K  and  where  G(u>)  is  given  by 


G(w)  =  exp 
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From  this  Theorem,  it  can  be  seen  that  the  canceller  output  spectrum, 
E„(1,<jj),  goes  exponentially  to  zero  in  K  when  y$  <  1/2  in  the  ambient  noise 
free  case.  Somewhat  surprisingly,  the  rate  does  not  depend  upon  the  shape 
of  the  source  angular  distribution,  but  only  upon  y<*>.  In  evaluations  of 
(1,oj)  given  in  Appendix  K  for  uniform  angular  source,  (d-47)  appears  to 
be  a  good  approximation  to  E^(1,uj),  even  for  modest  K. 

Note  that  if  ambient  noise  is  present ,  so  that 

H  (ot.u))  =  a  2 
e  n 

(D-62) 
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Then  the  extent  of  He(a,oj)  is  [-ir,nj.  This  corresponds  to  A  =  it,  which 
violates  the  conditions  of  Theorem  5.  (With  noise  present,  (D-62)  may  also 
violate  the  Lipschitz  condition).  However,  Theorem  3a  has  already  shown 
that  Ek(l,u))  approaches  a  positive  floor  as  10°°. 

In  the  extended  source  cancellation  problem ,  the  noise  free  case  is  of 
somewhat  limited  interest.  It  would  be  desirable  to  characterize  E  (m,u)) 

Co 

for  m  >  1  in  the  noise  present  case ,  which  would  provide  a  lower  bound  on 
cancellation  as  the  reference  array  is  moved  farther  from  the  primary.  This 
can  be  done  using  some  results  from  the  general  theory  of  stationary  random 
processes  which  do  not  depend  upon  the  Toeplitz  properties  of  the  covariance 
matrix. 


D-20 


In  the  noise  present  temporal  case ,  it  can  be  seen  that 


log  [Hx(a)]da  >  -<*> 


(D-63) 


so  that  the  process,  x(n),  is  not  singular.  If  the  spectral  distribution  of 
the  process,  say  Hx(a),  does  not  contain  any  jumps  and  is  differentiable 
everywhere,  then  x(n)  is  a  regular  random  sequence,  and 


Hx(U>= 


H  (a)da 
x 


(D-64) 


A  regular  random  sequence  has  a  moving  average  representation 


x(n)  =  ck  C(n-k) 

k=0 


(D-65) 


where  £(n)  is  a  sequence  of  uncorrelated  random  variables.  It  can  readily 
be  shown  that  the  minimum  mean  square  error  predictor  of  x(m)  given 
!x(n)| 
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so  that  the  prediction  error  is 


E.O") 


(D-67) 


Therefore,  if  the  c^'s  in  (D-67)  can  be  found,  the  irreducible  prediction 
error  for  m  steps  ahead  is  known . 

The  restriction  of  the  process  to  the  class  of  regular  random  processes 
will  not  limit  the  usefulness  of  these  results  in  the  extended  source  problem. 
As  an  aside,  however,  this  restriction  need  not  be  made  at  all.  As  long  as 
the  process  is  non-singular, 


£n  Hx(a)da  >  -«° 


(D-68) 


then  it  has  a  unique  decomposition , 

x(n)  =  u(n)  +  v(n)  (D-69) 

where  u(n)  is  regular,  v(n)  is  singular,  and  u  is  orthogonal  to  v.  Since 
v(n)  is  singular,  it  can  be  predicted  perfectly  given  the  infinite  past,  and 
need  not  be  considered  in  determination  of  the  minimum  prediction  error 
for  x(n).  As  a  regular  process,  u(n)  has  the  moving  average  representation, 


U(n)  =  X]  dk  £(n~k) 

k=0 


(D-70) 
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where  £(n)  is  an  uncorrelated  random  sequence.  The  irreducible  m  step 
ahead  prediction  error  for  u(n)  (hence,  for  v(n),  taking  into  account  the 
singularity  of  v(n))  is  then 


E  (m)  = 


(D-71) 


These  results  are  rigorously  derived  in  Lamperti  [16]  or  Doob  [18]. 

Even  without  finding  the  c^'s  of  the  moving  average  representation  (D-65) 
provides  a  useful  result  for  the  cancellation  of  a  narrow  extended  source. 

By  (D-67),  the  minimum  canceller  output  power  in  the  noise  present  case, 
E^tm.w),  is  a  non-decreasing  sequence  in  m.  The  cancellation  floor  there¬ 
fore  cannot  become  smaller  as  the  reference  array  is  moved  away  from  the 
primary . 

The  c^'s  in  the  moving  average  representation  of  a  regular  process  can 
be  found  using  a  factorization  of  the  spectral  density  Hx(a),  given  by  the 
following  theorem  from  Doob^1®^. 


Theorem  6:  Let  jx(n)|be  a  stationary,  zero  mean  regular  random  sequence 
with  spectral  density,  Hx(a).  Then  the  coefficients  of  the  moving  average 
representation  are  uniquely  determined  by  the  following 
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where 

bn  =  I  log[Hx(a)]e  ~jnada  (D-74) 
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By  expanding  the  right  hand  side  of  (D-73)  in  a  power  series,  or  by  differ¬ 
entiating  (D-73)  and  letting  z+0,  it  is  possible  to  obtain  explicit  expressions 
for  the  c^'s,  and  hence,  for 

m-1  0 

«.<■>-£  ioki 

k=0 

However,  the  expressions  resulting  from  this  procedure  are  quite  complicated. 
Supplement  II  to  this  Appendix  derives  expressions  for  c^ ^.Cg.Cg.c^, 
and  c,.  as  a  function  of  the  bk's  defined  by  (D-74),  which  allows  evaluation 
of  EM(n)  for  m  =  1  through  m  =  5. 

As  before,  the  temporal  linear  prediction  results  given  here  can  be 
applied  to  spatial  LMS  cancellation  by  replacing  Hx  (a)  by  Hg  (a.w)  given 
in  Theorem  la.  This  yeilds  the  following  result; 

Theorem  7;  Let  ER  (m.w)  be  the  output  spectrum  of  an  LMS  spatial  canceller 

as  defined  above,  and  let  He(a,u>)  be  defined  as  in  Theorem  la.  Denote 
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Supplement  I:  Strong  Szego  Limit  For  Circular  Arcs.  The  Strong  Szero 
Limit  Theorem  for  Circular  Arcs  is  given  in  [12]  as  follows; 

Theorem  8:  Let  f  be  a  function  satifying  f  (  9)  =  f  ( 2v  -  e)  which  is 
supported  on  a  closed  arc  a  <  8  <_2tt  -  a  and  which  when  restricted  to  that 
arc  is  positive  and  satisfies  a  Lipschitz  condition  with  an  exponent  greater 
than  1/2.  Let  D^[f]  be  the  Toeplitz  determinant 


DR[fI  =  det 
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Then  as  ,  one  has  the  asymptotic  formula 


DR [f]  'v  21/12  e3 


^(_1)(sinjj  4  EtF]^  4G[F]K(cos  |  )  (D- 


86) 


D-26 


where  if 


log  f(  0)  =  ^  ffc  ejke 
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one  defines 

F  [  f  ]  =  exp(fg) 
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and 

F( 0)  =  f  (2  cos  _1  (cos  2  cos  (D-90) 

[  19] 

Here,  4'(x)  is  the  derivative  of  the  Riemann  zeta  function.  The 

proof  of  this  Theorem  is  given  in  [12] . 

Theorem  5  requires  the  limit  of  the  Toeplitz  determinant  when  the 
function  f  satisfies  all  the  conditions  of  Theeorem  8,  but  is  supported  on 
the  closed  arc  -8  <  0  <  8.  In  this  case,  a  Corollary  to  the  Theorem  may  be 
stated  as  follows; 

Corollary  1 ;  Let  f  be  a  function  satisfying  all  conditions  of  Theorem  8, 
except  that  its  region  of  support  is  the  closed  arc  -8  <  0  <  8.  Then  the 
Toeplitz  determinant. 
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has  the  asymptotic  formula  as  K-*-«= 


DR[f]  %  2  1/12  e  35  '  (  1} 


(cosf) 


E [ F ] 2K  4G[F]K 


(sin 


where  G  and  E  are  as  defined  in  Theorem  8  but 
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where  f  has  support  on  [-8,8]  .  With  the  substitution  9  =  <t>  -tt, 
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with  g(tO  =  f(<t  -  T-) .  Note  that  g  has  support  on  [  tt  -  8,  it  +  S]  ,  and  letting 
a  =  if  -  g,  g  has  support  on  [a,  2r,  -  a] .  Thus  if  f  meets  all  the  other 
requirements  of  Theorem  1  on  [-8,  3)  ,  then  g  satisfies  the  conditions  of 
Theorem  1. 

In  order  to  relate 
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the  following  theorem  may  be  used. 


Theorem  9:  Let  two  square  matrices,  A  and  B  be  defined  by 


A  =  a(p ,q) ) 

\  /  p=0, 1, . .  .K-l 

q=0, 1, . .  .K-l 


and 


B 


4 


(-1)P  q  a(p-q) 


p=0, 1, 
p=0, 1, 


.K-l 

.K-l 


Then  det  (A)  =  det  (B) 


Proof:  Let  M  be  a  diagonal  matrix 


M  =  diag  ((-DP) 
p=0, 1,  •  •  .K-l 

and  note  that 

M-1  =  diag  ((-D  P)  =M 

p=0 , 1 , . .  .K- 1 

Now,  the  matrix  B  can  be  written  in  terms  of  A  as 

B  =  MAM-1  =  MAM 
Therefore 

det[B  ]  =  detfMAM]  =  det2[M]  det  (A] 
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Using  this  theorem, 
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where  g(4>)  =f(4>  -tt)  is  supported  for  <j>e  [a,  2tt  -  a]  with  a  =  tt  -  6.  Further, 
g  meets  all  the  requirements  of  Theorem  8,  so  that  limit  of  the  Toeplitz 
determinant,  (D-91),  can  be  found.  The  result  is  that 
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where  if 
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but 
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Let  the  argument  of  f  in  (A-23)  be  x,  so 
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so  that 

x  =  2  sin  1  |sin  ^cos  <J>  J  (D-106) 

and 

F( cj>)  =  g  ^2  sin  _1jsin  cos  <pj  j  (D-107) 

This  Theorem  allows  the  determination  of  the  asymptote  of  the  ratio  of 
Toeplitz  determinants  as  follows: 

Corollary  2:  Let  f  be  a  function  satisfying  the  conditions  of  Corollary  1, 
and  let  D^ff]  be  the  Toeplitz  determinant 


D-32 


I 


}j(p  q)e  f(0)  d6\ 


p=0, 1, 
q=0, 1, 


.K-l 

.K-l 


(D-108) 


Then  as  K-*-» ,  the  ratio  D^+^[f] /D^[f]  has  the  asymptotic  formula 
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with  G  and  F  as  defined  in  Corollary  1. 
Proof:  By  Coilary  1, 


(D-109) 


1  /  \  -  ■  v  41 

DK+i[fl  2T2  e3  r(_1)  (cos  f)  4  E[F1  2  (K+1)  4GtplK+1  (sin  I  ) 
DK[f]  ^  F  ~  72 

212  e3^  '(_1)  (cos  |-)  E[F]2K  -1/4  G[F]K  (sin  |-) 


-1/4 

1  •  2K+1 
sin  f 

(D- 110) 

_/K+l\ 

G[F] 

-(-TT) 

V  z/ 

~G[F] 

(si")  (sin 

l)“ 

D-33 


md 


Values  of  m 


From  equation  (D-78), 


CO 


n=0 


B  (Z)  = 


n=l 


bnZ 


n 


Now,  denote 


c(m)(Z)  =dmC(Z) 


dZ 


m 


B(m)(Z)  =dmB(Z) 


dZ 


m 


and  note  that  for  m  >  0, 


C(m)(0)  =  m!  c  *  ,  B (m) (0)  =  m!b 
m  m 


(D-lll) 


(D-112) 


(D-113) 


(D-114) 


(D-115) 


D-34 


Using  (D-lll),  and  dropping  the  arguments,  Z,  for  convenience 


C(1)  =  exp [B]  B(1) 


C(2)  =  exp[B]  |  B(1)  +  B(2)  J 


C(3)  =  exp[B]  B(1)  !b(1)  +B(2) 


+  exp[B]  <2B(1)  B(2)  +  B(3) 


= exp[B]  B(1)  +  3B(1)  B(2)  +  B(3) 


C(4)  =  exp[B]  B(1)  |b(1)  +  3B(1)  B(2)  +  B(3) 


(1)  B(2)  +  3B(2)  +  3B(1)  B(3)  +  B(4)! 


+ exp[B]  |3B'‘'  B  +  3B  +  3BV/BV“'+B 
exp[B]  (b^1^  +  6B^  B^2^  +  4B(1^  B^ 


+  3B(2)  +  B(4) 


C(5)  =  exp[B]  |b(1)  +  10B(1)  B(2)+10B(1)  B(3) 


+  15B 


Q)  g ( 2)  +  10g(2)  g(3)  +  5B(1)  B(4^  +  B(5H 


Therefore,  setting  Z  =  0  and  using  (B-5)  and  (B-6), 


m 


m  <  - » 

■  s  a 


cm-k 


(D-116) 

(D-117) 


(D-U8) 


(D- 119) 


(D-120) 


(D-121) 


D-35 


I 


where,  from  (D-117) 


cQ  =  exp  (bQ/2) 


(D-122) 


Some  specific  c  's  are  then 
m 


C1  =  c0  bl  =  bl  exp  <V2) 


(D- 123) 


C2=Ci*bl+C0b2 


=  (C0  bl  +  c0  *2)  =  Cbl2  +  b2]  exp  [b0/2J 


4c  *  4c  4c 

C3  —  C2  bl  +  2cl  b2  +  c0  b3 


[bl2  +  b2]  bl  c0  +  2  bl  b2  C0  +  C0b3 


=  [b  1  +  3^  b2  +  bg]  exp  (bQ / 2) 


C4  =  c3bl  +  3c2b2  +  3clb3  +  cob4 


(D- 124) 


(D-125) 


-  ‘»13  +  ^  *  S1  bl  Co‘  *  3b2  'bl2  +  »2>  =02 


4  4 

+  3  b,  b_  cn  +  cn  b  . 
1  3  0  0  4 


(D-126) 


=  [bj4  +  6bj2  bg  +  4bj  bg  +  Sbg2  +  b4]  exp  [bp/2] 


D-36 


*  *  *  *  *  * 

cc  =  c .  b1  +  4e0  b0  +  6  c0  b0  +  4c,  b  4  +  crt  b- 
o  41  6  6  l  3  14  05 

=  bi  [b^  +  ebj’bjt^bjtsbj^b^i  cj 
*  4b2  [bj3  ♦  3b,  b2  *  b2]  c0*  *  6b3  [bj2  *  b2]  0* 

+  4b4bl°0tb5C0 

=  [b^  +  10b13  b2  +  lOb^  b3  +  15b:  b22  +  10b2  b3  +  5bx  b4 
+  b5]  exp  [bfl/2]  (D-127) 


From  (D-76) 


ni-  1 

Eco(m,u.)  =  y] 

k=0 


(D- 128) 
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APPENDIX  E;  CANCELLER  PERFORMANCE  IN  TERMS  OF 


The  power  spectral  density  of  the  optima]  canceller  output,  from 
equation  (G-17)  of  Appendix  C  is  given  by 


Ek(u;)  =  S^fco)  -  1(oj)s1(w) 


(E-I) 


Sn(oj)  =  PSD  of  the  primary 

Sg  ( co )  =  CSD  matrix  of  the  reference  hydrophone  outputs 

—1  (u) )  =  CSD  vector  between  the  primary  and  reference  hydrophone 
outputs 


For  the  far  field  extended  source  described  in  Appendix  A,  assuming 
a  uniformly  distributed,  narrow,  spatially  uncorrelated  source,  the  statistics 
of  the  hydrophone  outputs  are 


S , , (u) )  =  20  a  2(uj  )  +  a  2 
it  in 

.cod 

2 ,  .  J  COS  *o  „+ 

=  Cj  (co)  e  G  3^ 


where 


(E-2) 


(E-3) 


(E-4) 


s  Q  =  [so(o,w),  sQ  (1,0), ...  s  0(  K  - 1 ,  CO )  1 


sin  (J-pl^^O  sin  <J>  1 

Vp’u)  =  rr  77  r  jd  .  ~y 

2  (J'p)  ~rsmK_ 


(E-5) 


E-l 


and 


S 

e 


n 


'(u))I  +  o,  (  m  )G 


S  .  .  G 

eo(ui  ) 


(E-6) 


where 


Seo(“> 


(sin  (p-q)  $  sin  <pQ  \ 

|[(p-q)  ^sin  <fo]  / 


and 


G  =  diag 
P 


■  .  aid 

-]  P  COS  <fi 

e 


o 


p  =  0,1, ...K-l 
q  =  0,1, ...K-l 


(E-7) 


(E-8) 


The  inverse,  Sg^w),  can  be  obtained  in  terms  of  the  eigenvalues  and 
eigenvectors  of  Sgo(w ),  which  can  be  defined  in  terms  of  the  Discrete 
Prolate  Spheriodal  Sequences  (DPSS)  as  discussed  in  Slepian  [  13 3 . 

Note  that  (E-7)  can  be  written 


q  =  0, 1, . . .  ,K- 1 

with 

d  . 

Y  =  t  sin  d> 

'A  T0 

W  = 

\  =  wavelength  at  frequency  ou 


S  (  w )  = 
eo 


i  / gin  [2TrW(p-q)3  \ 

y  ^  *  (p-q)  / 


E-2 


Now,  the  DPSS  satisfy 


K-l 


Esin  [  2  7T  ( n  -  m)W  ]  y(k)  (K  w)  =  \  (K  ,W)V(k)  (K  ,  W)  (E-10) 

Tt  (n-m)  m  kv,/n 

m=0 

for  k=0, 1 _ K-l,  where  Jv^k)  (K,W)j  n=Q  +1  +2  is  the  k— DPSS  and 

where  X^(K,W)  is  the  associated  eigenvalue.  Equation  (E-10)  for  k=0,l,  ...K-l 
can  be  written  in  vector  form  as 


s^Vk(K,W)  =  Xk(K,W)  V^k)  (K,W) 


where 


?n  = 


sn(q) 


=  [Vo)’  sn(1>'  »n(k-»] 


sin  2tt  (n-q)W 
Tt(n-q) 


and 


Vk  (K,W) 


V^k)  (K,W), 


..,  V 


(k) 


t  (K , W) 


(E- 11) 


(E- 12) 


(E- 13) 


(E- 14) 


E-3 


Therefore 


”  + 

s 

V(k)(K  ,W) 

-o 

O 

+ 

s 

-1 

V^(K , W)  =  \k(K,W) 

V^k)(K.W) 

+ 

*k-l 

V<k_>  (K , W) 

k=0, 1,  . .  .K-l  (E - 15) 


Comparing  (E-15)  with  (E-7),  it  can  be  seen  that 


Seo<^  =  y 


S+(u>) 


^K(uj)-1 


if  W  =  y$ 


(E- 16) 


with  y  =  fd/c  sin<f>o  and  uj  =  2tt  f.  Therefore,  multiplying  both  sides  of  (E-15) 
by  1/y  gives 


s«a(k^)  ^  »k(K^*)  4  (“•*♦) 


(E- 17) 


In  this  form,  it  can  clearly  be  seen  that  the  set  of  eigenvectorsj^CK  ,W)  j 
k=0, 1,  . . . ,  K- 1,  of  S  (u> )  are  the  V,  (K  ,W)  given  by  (E-14)  with 

GO  K 

associated  eigenvalues  [-  \  (K,W)].  Noting  that  G  is  diagonal  and  that 

Y  k 

G  *  =  G+,  it  can  be  seen  that 


(E- 18) 


[G+So  G]  G+  V.  (K,W)  =  X  ( K  ,  W )  G+  V,  (K,W) 

GO  K  K 

so  that  the  G+  Sgo  G  has  the  same  eigenvalues  as  Sgo,  but  its  eigenvectors 
are  |g+V  (K,W)J  . n  , .  Consequently,  Se  has  the  same  eigen- 

|  ■■  K  |  K  U  i  1 1  •••  K  X 

vectors,  but  its  eigenvalues,  rk|(K,W)j  k_Q  ^  ^  are 

2 

2  °r 

rk  ( K , W )  =  +  -L-  X  (K,W)  (E - 19) 

Therefore, 


K  - 1 

SQl  =  Y,  N^'1  (K,W)  G+ Vk  (K,W)Vk+  (K,W)G  (E-20) 

k-o 

and 

K-i  2 

sX'1*,  =  T  rk_1  (K,W)  |s*  G+  Vk  (K,W)  |  (E-21) 

k-0 


Using  (E-3)  gives 


st  G+  V,.  (K , W) 


-1 


-k 


2  +  + 

°!  GG  Yk  6 


aid  T 

+]  - J  cos 


=  o^s^Vj^CK  ,W)e 


. ,  wd  T 

+1  —  J  cos4>0 


but 


(E-22) 


(K’W) 


K-l 


E 

m— n 


sin[27r  (J-m) y$  ] 
t(J-m)2TTY] 


v(k) 

m 


(K,W) 


E-5 


y 


*  £ 


-Sm-2(j(m)m)VV1  vmk)  (K*W)  =  “  Xk(K,W)VJk)(K,W)  <E'23> 


by  applying  (E-10)  again.  Substituting  this  in  (E-21)  yields 


,  +Q  ‘  1 

— lSe  -1 


K- 1  a  4X  2(K,W) 

£  2  2  k  } -  1VJk>  <K-W> 

k 


Therefore 


2  2  1  24»a2  A  2(K,W)  2( 

E  =  0n“  +  0I  P  y  X  - T  ~~ ~  2 -  |V^k)(K,W)  }(E-24) 

‘  2w°n  +  2*o2X(K,W)  J  * 


For  the  moment,  suppose  that  the  noise  field  is  totally  due  to  the  extended 


source,  on  =  o,  so 


E  =  Or*  2*  -  - 
1  Y 


*  £ 


A  (K.W)  j  Vj(k)(K,W) 


( E  -  2  5) 


A  good  check  on  this  result  is  to  let  J  =  Jq  for  some  0  Jq  <  K  -  1 ,  so 
that  one  reference  is  coincident  with  the  primary  hydrophone.  This  should 
give  perfect  cancellation.  In  that  case 


u  i 

Y.  Ak(K,W)V^k)  ( K ,W)Vjk)  (K,W) 


E-6 


(K,  W) 


K-t 

e 

^  ^  sin[  2n(  JQ 

-  2)W] 

tt( J  -  l) 

k=0 

to 

II 

o 

c 

K- 1 

E 

sin[  2 7T (  Jq  -  £)W] 

K  - 1 

E 

tt(J  -  l) 

2=0 

o 

k=0 

V(k)  (K,W)V(k) 
Jo 


V£(k)  (K  ,W)  Vjk)  (K,W) 
o 


(E-26) 


But  the  DPSS  are  orthogonal  when  index  limited  to  [0,K-1] 


K- 1 

y;  v£(k)  (K,W)Vj(k)  (K,W)  =  6U 

t _ n  O  1 


so 

K~ 1 

y  X^K.WlVj^5  (K,W)Vj(k)  (K  ,W)  =  2W  =  ^ 


(E-27) 


(E-28) 


Substituting  this  in  (E-26)  gives  E  =  0  indicating  perfect  cancellation  of  the 
source  as  it  should  when  a  reference  is  coincident  with  the  primary. 
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APPENDIX  F:  RESULTS  FOR  THE  UNIFORMLY  DISTRIBUTED  FAR  FIELD 
MODEL  BASED  UPON  TOEPLITZ  PROPERTIES  OF  THE  HYDROPHONE  CSD 

When  the  extended  source  is  narrow  in  the  sense  of  (A-24)  and  spatially 

uncorrelated,  substitution  of  (D-53)  into  (D-49)  yields  the  irreducible  canceller 

output  spectrum ,  as 


E  (l,o))  =  exp 


W  jf  ’  ln  d“ 


=  exp 


A />[*«>  ♦£ 


/  a+2Tik  +  6  \ 

I  \  2rvT"  °) 


Rect 


k=-°° 


0!+2TTk 


2ny<p 


da 


When  y<1>  i  i,  using  (D-54),  this  reduces  to 


Ew(l ,oi)  =  on(aj)  exp 


/ 


ln 


1  + 


2  . 

yo_(oj) 


Sj(w,  <j>  + 


d<j> 


(F-2) 


F-l 


In  the  case  when  extended  source  is  uniformly  distributed  on  [$  -4> , 
as  in  (A-40),  then  (F-l)  and  (F-2)  can  be  evaluated  in  closed  form. 
Since  the  uniform  source  distribution  is  a  reasonable  representation  for  a 
number  of  sonar  problems,  this  result  is  quite  important.  Let 

(  O?  (to)  ,  $£[$  "  $  ,  $  +  <5>1 

S ,  ( <(>  ,to)  =  ) 

(  0  ,  otherwise 

and  define  an  integer,  n,  such  that 

~<  Y*  <  ^  (F-3) 


The  positive  alias  in  (a, to)  as  defined  by  (F-l)  centered  at  k/y  extends 

e 

from  (k/y)  -  to  (k/y)  +  <I> ,  so  if 


k  .  .  1 

—  $  <  - 
y  2y 


th 


(F-4) 


lo"'[^-27j 


the  k  positive  alias  extends  through  the  interval  of  integration 
If  (F-4)  is  met,  the  kth  negative  alias  also  extends  through  the  interval  of 
integration.  Let  kQ  be  the  number  oi  aliases  extending  completely  through 
2y’2^  .  Using  the  above,  kQ  =  [yf  -  I/2J  where  [xj  is  the  largest  integer 

■*  **  XU 

<  X.  The  next,  (kQ  +  1)  ,  alias  extends  from  (kQ+l)/y  -  $  to  (kQ  +  l)/y  +  <{> 

First,  consider  n  (in  (F-3))  even,  so  that 


k0  =  h*  -  |J  =  j- 1 


(F-5) 


Therefore 


(F-6) 


(k  +1) 

-  -  $  =  -  $  <0 

Y  2y 


since,  by  (3-98) 


Y*  >  £ 

th  n  X 

Hence,  the  (kQ  +  1;  positive  alias  extends  from  2^  ~  ®  2y  ’  w^e  t^le 

negative  alias  extends  from  — p  to  -  ^  +  $  ,  overlapping  as  shown  in  Figure 

F-l. 


Figure  F-l.  Folded  Spectral  Density  for  Uniform  Source,  n  even 


Using  this,  there  are  2kQ  aliases  (kQ  positive  plus  kQ  negative)  extending 
through  [-1/2y,  1  / 2y ]  plus  the  primary  spectrum  and  the  sum  of  the  (kQ+l)th 
aliases.  Hence,  in  (F-l) 


F-3 


This  result  is  identical  to  that  obtained  by  Slepian  for  the  terrporal  linear 


prediction  of  a  perfectly  bandlimited  random  process  in  uncorrelated  noise. 

However,  Slepian  considered  the  problem  in  terms  of  the  Discrete  Prolate 

Spheroidal  Sequences  (DPSS).  The  DPSS  and  their  eigenvalues  are  associated 

with  the  hydrophone  CSD  matrix  produced  by  the  uniform  source.  This  results 

depended  upon  the  asymptotic  evaluation  of  several  complex  integrals  and  did  not 

demonstrate  the  effect  of  aliasing  of  the  wave  number-frequency  spectrum, 

H  (o,w),  on  the  linear  prediction  error. 

0  * 


If  n  is  odd  then 


and  the  <kQ  +l)th  positive  alias  extends  from  (ko+l)/V  -$  to  (kQ+l)/\  +$, 
where 


(ko+1)  -  <i>  =  n+1  .  (F-12) 

—  ~2Y  -*>0 

since  ~Y<fr<n+l .  Then  the  (k  +  11th  aliases  do  not  overlap,  as  shown  in 

n  ® 


Figure  3-2.  Therefore,  by  the  same  argumei '  :  used  the  even  case. 


Figure  F-2.  Folded 'Spectral  Density  for  Uniform  Source;  n  odd 
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F.l  Interference  Only  Case 

As  mentioned  above,  the  case  where  the  noise  field  consists  of  only  a  nar¬ 
row  extended  source  with  no  ambient  background  does  not  occur  in  sonar. 

2 

However,  consideration  of  the  ambient  noise  free,  0^(00)  =  0,  case  yields 
insight  into  the  behavior  of  the  LMS  optimal  canceller  in  the  field  produced 
by  an  extended  source .  Consider  the  irreducible  canceller  output  spectrum , 

E  (l.oo)  when 


a  ( w)  =  0 
n 


It  is  easily  shown  that  (F-9)  reduces  to 


E  (l.oo) 


2<J>  0.  (ou) 


|o*(w)  =  0 


i  i  i  r 

—  n  I1  n  ) 


(F-16) 


with  n  defined  as  the  integer  such  that 


n  .  t  .  n+1 

Y  $  <  — 


(F- 17) 


Therefore,  if  -y <t>  <  so  that  n=0, 


E  (l.oo) 


=  0,  Y$  <  g 


(F- 18) 


1  O^(oo)  =  0 


and  the  extended  source  can  be  cancelled  completely,  at  least  asymptotically 
in  K,  the  number  of  references.  On  the  other  hand,  if  >  1/2,  E,  (l,oo) 

K 

approaches  some  positive  value  as  K  becomes  large.  The  significance  of  this 


where 


(F-15) 


which  is  the  same  as  (F-8).  Therefore,  n  odd  yields  the  same  result. 


F-8 


Ml  3 


I 


fact  can  be  appreciated  by  observing  Figure  F-3,  which  shows  E^d.u:)  as  a 
function  of  for  y4>  Y.  1/2.  Clearly,  when  yP  >  1/2  the  cancellation  of  the 
extended  source  is  extremely  limited.  The  condition  yP  <  1/2  is  therefore  a 
spatial  Nyquist  condition,  with 

d  .  . 

Y  =  —  sm  p 
a  o 

the  spatial  sample  interval,  and  <t> ,  the  source  extent,  analogous  to  bandwidth 
in  the  temporal  case.  This  result  should  be  no  surprise  given  the  equiva¬ 
lence  of  the  temporal  linear  prediction  problem  and  the  spatial  cancellation 
problem  developed  in  Appendix  D. 


Figure  F-3.  Irreducible  Canceller  Output  Power 
vs.  W  for  Uniform  Extended  Source 
in  Noise  Free  Case 

Using  the  Strong  Szego  Limit  Theorem  for  Circular  Arcs  [12],  as  discussed 
in  Appendix  D,  it  is  possible  to  determine  the  asymptote  in  K  of  Ek  (1,  a>) 
for  the  uniform,  narrow,  spatially  uncorrelated  source.  For  this  source  in 
the  noise  free  case  the  hydrophone  CSD  is,  from  Appendix  A. 


Sk(w)  =  se<P.P+*-w> 

=  Oj(u>)  stolMZiSm 

I  TT(p-q)Y 


(F- 19) 
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But  from  (D-38) 


sk<”>  * 


.X  [ 

2tt  I 

•'-TT 


H  (a  ,<jj)  da 
e 


so  if  <  tt,  choosing 


Oj  (w) 

H  ( a  ,  to )  =  — — —  Rect 

G  I 


[&]• 


ac[-  tt,  tt] 


where 


Rect 


[*]■ 


[!•?] 


1,  te 
0,  elsewhere 


(F-20) 


yields  the  correct  S^Cw).  This  function,  He(a,w)  satisfies  the  conditions  of 
the  Strong  Szego  Limit  Theorem,  stated  as  Theorem  5a  in  Appendix  D,  so 
that 


Ek(l,u>)  '  G(w) 


£  sin  tty  ^  J 


2K+1 


where  '  denotes  asymptotic  in  K  and,  from  (D-58)t 


G(oj)  =  exp 


/TT 

-TT 


log 


-i -  Rect  »  2  sm 


[sin(-rfW)  cos  9] 
4ttW 


(F-21) 


de 

(F-22) 


F-10 


where 


|  2  sin  ^suUttW)  cos  Q]  |  _!l,  |2sin  1[sin(7TU')cos  9]  |  £  2ttW 
Kect  |  4ttW  1  )0,  otherwise  '  (F-23) 


Figure  F-4.  Relationship  Between  ttW 
and  sin  1[sin(Tr\V)  cos  9] 

But  from  Figure  F-4,  it  can  be  seen  that 


sin  *[sin(TiW)  cos  0]  4  |nw|. 


all  0£(-tt,tt) 


or 


2sin  sin  (itW)  cos  0]  < 


|2ttW|  ,  all  0e[-tt,  tt) 


which  implies  that 


G(co)  =  exp 


2tt 


2, 

0.  (w) 

log  -  d  0 

B  y 


2,  , 
ctj  (to) 

y 


(F-24) 


so  that 


Ek(1,uj) 


O?(to)  r  -]2K+1 

—  sin  TTW  I 


(F-25) 


F-ll 


Thus  for  large  K,  the  canceller  output  spectrum  goes  to  zero  exponentially  in 
K,  with  the  rate  dependent  upon  sin  (-nW). 


F-12 


1 

I 


F.2  Interference  Plus  Noise  Case 

When  ambient  noise  is  present,  the  change  in  the  irreducible  canceller 
output  spectrum,  (l,w),  as  W  transitions  from  W<g  to  W>j  is  not  quite  as 
dramatic,  but  still  supports  the  notion  of  W=i  as  a  spatial  Nyquist  criterion, 

particularly  at  high  interference  to  noise  ratios.  Figures  F-5,  F-6,  F-7,  and  F-8 
show 

Er(1,  u) )  =  (F-26) 

20<|%) 


as  a  function  of  W  for  interference  to  noise  ratios  of  40,  30,  20,  and  10  dB, 
respectively.  The  interference  to  noise  ratio  is  defined  as 


ENR  = 


2  $<rr2M 

<r  2(u>) 
n  ' 


(F-27) 


When  noise  is  present  it  can  be  seen  that  it  may  be  necessary  to  choose  W  signifi¬ 
cantly  less  than  1/2  to  achieve  acceptable  cancellation,  even  with  a  larte  number  of 

reference  hydrophones.  Of  course,  from  Appendix  D,  the  extended  source 

o 

cannot  be  cancelled  to  the  noise  floor  when  (r^  (w)  ^  0. 

| 

2  2 

For  high  interference-to -noise  ratio,  2$ o-j.  (w)/irn  (w)>>  1,  so  that  in 
(F-10)  for  n  >  1 


Sain 


and  (F-9)  is  approximately 


E«(1,  w)s 


2®(Tj  (w) 


nM) 


2ui-n 


(F  -28) 


(F-29) 
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FIGURE  F-5.  Irreducible  Canceller 
Output  Power  vs  W  for 
Sin  x/x  CSD 


Hence,  for  high  INR  andY$>l/2  (n>l),  the  interference  plus  noise  case  reduces 
approximately  to  the  noise  free  case. 


F-18 


APPENDIX  G:  CANCELLATION  OF  A  SOURCE  PRODUCING  AN 
EXPONENTIAL  HYDROPHONE  CROSS -SPECTRAL  DENSITY  MATRIX 


APPENDIX  G:  CANCELLATION  OF  A  SOURCE  PRODUCING  AN 
EXPONENTIAL  HYDROPHONE  CROSS-SPECTRAL  DENSITY  MATRIX 


REFERENCES 


Figure  G-l:  Array  Geometry  for  Cancellation 
Consider  a  line  array  of  K  reference  sensors  used  to  cancel,  in  the 
minimum  mean  square  sense ,  an  interference  from  the  output  of  a  primary 
hydrophone  which  lies  along  the  same  line  as  the  references.  Let  the 
references  be  spaced  d  feet  apart,  and  let  the  space  between  the  nearest 
reference  and  the  primary  be  md  feet,  as  shown  in  Figure  G-l.  Suppose 
that  the  ambient  noise  field  produces  zero  mean  noise  with  spectral  density, 

9 

on“(u)),  at  the  hydrophones,  uncorrelated  from  hydrophone  to  hydrophone. 
In  addition,  the  noise  field  contains  an  extended  interference  which 
produces  zero  mean  noise  at  each  hydrophone  with  the  cross-spectral 
density  between  hydrophones  spaced  kd  apart  given  by 


0lWB‘“>  I"  id 

Then  the  cross-spectral  density  (CSD)  between  the  p— -  and  q— 
hyrdophone  outputs  is 

se(p,q,cj)  =  on^(o))6pq  +  Oj^(w)e  lp  3 


(G-l) 


( G-2) 


G-l 


while  the  CSD  between  the  p —  reference  and  the  primary  is 


S  (p  ,  J  ,  uj  ) 
e 


2, 

r  ( u;)e 


j)  |p-J  Id 


( G-3) 


Appendix  C  developed  the  power  spectral  density  of  the  minimum  mean 
square  canceller  using  K  references  as 


det  [Sh(K,o))J 
EK(w)  =  det  [S  (K,u>)] 


where 


S  (K,co) 


(se < p , q ))  p=o, i, . .  .K-i 
p=0, 1, .  .  .K-l 


and 


Sh(K,u)) 


Se(K,io)  s^(K ,  u>) 
▼ 

s^+(K,w) 


( G-4) 


(G-5) 


(G-6) 


with 


s.iK.w)  =  [s  (o,J,ut),  s 
—lee 


se(K-l,J,w)] 


(G-7) 


bll  " 


n 


'(w)  +  a. 


2(oi) 


I 


(G-8) 


# 


G-2 


For  the  model  being  considered  here,  it  can  be  seen  that 


G.l  Special  Case:  m  =  1 


In  the  special  case  where  m  =  1,  both  Se(K,u)  and  S^(K,u))  are  Toeplitz 
matrices.  This  allows  use  of  the  results  of  Szego^  ^  ,  as  discussed  in 
Appendix  D,  to  determine  the  irreducible  canceller  output  power  for  the 
LMS  canceller  with  a  source  producing  an  exponential  CSD.  It  will  be 
convenient  in  this  Appendix  to  explicitly  indicate  the  value  of  m  in  the 
notation  for  the  canceller  output  power,  so  let 

EK(m,u))  =  Ek(oj) 

K  md  =  distance  between  reference  array  and  primary 
From  Appendix  D ,  let 

Sk  =  Se(p,p+k,uj)  any  p 

so  that 


Then  from  Appendix  D,  the  irreducible  canceller  output  power  at  frequency 
w  is 


G-4 


where 


H,(C.U)=  £  Vlk“ 

k  =  -oo 

In  the  special  case  of  an  exponential  CSD,  (G-2), 

=  66  (k)  +  p  I 

with 

II,  k=0 

6(k)  =  j 

l0,  otherwise 

then 

det  [Sh(K,u)>]  =  tOj2(o))] K+1det  [Se (K+l, u)] , 
det  [Se(K,o))]  =  [Oj2((jo)]K  det  [Se(K,w)]; 

and 

CO 

He(a,u)  =  B  +  ^  Pfk!ejka 

k=-a> 

i  2 

=  6  + - hP _ 

1-2  pcos  a+p2  (G-14) 

=  )  +  ( 1—  p  )  -  2  8  pcos  a 

l-2pcosa+p“ 


G-5 


therefore 


In  H  (a,cu)  -  -ln(  1+p2)  -  ln(  1  -  ■  — ,  cos  a) 

i  V 


(G-15) 


+  In  [  8(  1+p2)  +  (l-p2)]  +  ln  (1 - - s-  cosa) 

B(l+p  )  +  (l-p^) 


Fortunately , 


In  [  1-Acos  a]  da  =  2-rrtn 


1+ 

2 


so 


and 


(G-16) 


(G-17) 


2Bp 


B(  1+p2)  +  (1-p2) 


cos  a 


da  =  In 


28p 


B( 1+p  )  +  (1 


s) 


2 


r 


H(a)da  =  -  In  ( 1+p  )  -  In  - *- 


+  In  8( 1+p2)  +  (1-p2) 


)]  +  In 


1  +M  - 


2Bp  \ 

(1+p2)  +  (1-p2)/ 
2 


=  -  In  - 


(1+P2)  +  'Z ( 1+p2)2  -  ( 2 p ) 2 


+  In  - 


8(l+p2)  +  (1-p2)  +  >/[g(l+p2)  +  (1-P2)]2  -  (28p)2 


6(  1+p2)  +  (1-p2)  +  /Tb(1+p2)  +  (1-P2)]2  -  (2Bp)2 


(G-19) 


Therefore 


[a  2(u>)]K  1  det  [S  (K+l,w)] 

£_(!,«)  =  — 4 - * - - - 

[Oj  (w)]  det  IS  (K,u)] 


aI2(a))exp  ^  /  In  He(a,u>Xlcx| 


_  2,  x  8(  1+p2)  +  (1-p2)  +'/[B(l+p2)  +  (1-p2)]2  -  (Bp)2 
°I  - 2 - 


(G-20) 


Now  as  p  -*■  o  (for  example,  as  the  hydrophones  move  farther  apart, 


d  , 

lim  E^d.oj)  =  Oj 2 ( oj )  3+-— ^  "  Oj2(w)  +  on2(w)  (G-21) 

p  -*■  o 

which  shows  that  there  is  no  cancellation,  as  would  be  expected.  On  the 
other  hand,  as  p  -*•  1,  the  references  become  coincident  with  the  primary  and 

lim  E  0,0))  =  a  2(w)  (G-22) 

cw  n 

p  - 1 

which  indicates  perfect  cancellation. 

Next,  consider  8  -*■  o,  indicating  the  noise  free  (interference  only)  case. 
Then 

lim  E^d.uO  =  aj2(w)  (1  -  p2)  (G-23) 

8  ->  o 

Therefore,  even  in  the  absence  of  a  noise  floor,  the  interference  is  not 
cancelled  completely.  This  is  a  significant  difference  between  the  exponential 
cross-spectral  density  and  the  sin  x/x  CSD  considered  in  Appendix  F. 

G.2  General  Case;  m>l 

Return  to  equation  (G-10),  and  observe  that  if  the  Kth  column  of 
SH  (K,w)  is  multiplied  by  pm  and  subtracted  from  the  last  column  (which 
does  not  change  the  determinant),  it  becomes 


G-8 


2,  x 

Oj  (w) 


m+K-1 


S  ( K ,  oj  ) 

e 

2 .  , 

Oj  (w) 


m  „ t  o  2m 

.  p  e+p  -p 


(G-24) 


Repeating  this  on  the  K  row  and  the  last  row  gives 


2  . 
Oj  (u) 


Se(K,oj) 

0^(00) 


(G-25) 


a  m  ,  „  ^  o  2iti  _  2m) 
.  o  -Bp  (B  +  P  -  P  +  Bp 


Under  these  operations,  the  determinant  is  unchanged,  so  det  [S„(K,cj)]  is 

ri 

equal  to  the  determinant  of  (G-25),  which  is  easily  shown  to  be 
det[SH(K,w)]  =  [aI2(u>)]  IB(l+p2m)  +  (l-p2m)]  det  [Se(K,u)] 


-  [or2(w)]2  B2  p2m  det  (Se(K-l,a))] 


(G-26) 


det  tS^K.w)) 

EK(m,w)  =  det  [Se(K,w)]' 

det  [S  (K-l.u)] 

2,  \ro/ix  2m,  .  ,,  2m,,-,  r  2,  ..2.2  2m  - 

=  Oj  (a»)[3(l+P  )  +  (1-P  >1  tor  (a,)]  6  P  det  [S^K^)] 


(G-27) 


G-9 


The  limit  of  the  ratio  of  determinants,  using  (G-12)  and  (G-14)  is 


det  [S  (K-l.u)]  2-1  I  i  r 

^  di—Se(K, -)T-  =  [°I  (w)1  exp  "2?  J  ln  He(a’w)do 

-  7T 


=  [aI2(ui)]_1 


j^3(l+p2)  +  (1-p2)  +  v^Bd+p2)  (1-p2)]2-  (2gp)2  j 


(G-28) 


Therefore 


Em(m,io)  =  Oj  2(oj)  j  [  B(  l+p2m)  +  (1-P2m)] 


2B2  p2m 


6(l+p2)  +  (1-p2)  +  Vt6(l+p2)  +  ( 1-p2) J ^  -  (2Bp) 


When  m  =  1,  this  reduces  to  (20),  as  it  should. 


( G-29) 


Now,  when  p  -*  o,  this  goes  to 


lim  E^im.oj)  =  Oj2(w)  [3+1]  =  Oj2(o)  +  on2(u)) 
p  -*■  o 


(G-30) 


again  indicating  no  cancellation,  as  would  be  expected.  When  p  -*■  1, 


lim  E  (m,<jj)  =  o  (w) 

p  -  1  °°  n 


( G—  31) 


which  shows  that  the  interference  is  completely  cancelled  when  one  of  the 
references  is  coincident  with  the  primary.  Finally,  when  6  -*•  o,  the 
noise  free  case, 


lim  E  (m,u)  =  aT  2(u>)  [l-p2m] 

„  °°  l 

S  ■*  o 


(G-32) 


G-10 


9m 

The  fact  that  the  canceller  output  power  has  the  lower  bound  (1-p  m) 

2 

Oj  (uj)  in  the  noise  free  case  can  be  explained  in  terms  of  the  special  nature 

of  the  exponential  CSD.  The  sequence  of  reference  hydrophone  outputs  at 

a  single  frequency  can  be  shown  to  be  a  first  order  Markov  sequence.  As  a 

result,  in  the  noise  free  case,  the  optimal  canceller  uses  only  the  output  of 

2m  2 

the  reference  closest  to  the  primary  and  achieves  the  bound  (1-p  )  Oj  (to) 

with  that  single  reference.  Adding  further  references  does  not  reduce  the 
canceller  output  power. 

The  amount  of  cancellation  achieved  by  the  canceller  is  best  characterized 
by  the  minimum  relative  canceller  output  power,  given  by 

E^fm  ,co) 

ER(m,w)  =  - 2 -  (G-33) 

Oj  (co)  +  (to) 

2  2 

where  Oj  (to)  +  an  (to)  is  the  output  power  of  the  primary  hydrophone 
without  cancellation.  Noting  that 

ot2(co)  +  a  2(io)  =  o  2(to)  (1  +  6) 

Ini 


gives 


ER(m ,to)  =  j-Lg-  j  [  g( l+p2m)  +  (l-p2m)]- 


2B2  p2m 


(G-34) 


8(l+p2)  +  (1-p2)  +  \[ B( 1+p2)  +  (1-P2)]2  -  (26p)2 


G-ll 


Figures  G-2  through  G-6  plot  the  minimum  relative  error  given  by  (34)  as 
a  function  of  p  for  values  of  10  log  (1/6)  of  50,  40,  30,  20,  and  10  dB 
respectively.  It  can  be  seen  that  cancellation  to  near  the  noise  floor, 

9 

(w),  requires  that  the  primary  and  references  be  very  highly  correlated. 
This  is  again  due  to  the  Markov  nature  of  the  hydrophone  outputs  under  the 
exponential  CSD.  At  high  interference  to  noise  ratio,  the  optimal  canceller 
uses  only  the  reference  closest  to  the  primary  and  the  output  power  is 
proportional  to  (1-p  m).  Additional  references  have  no  effect,  so  the  can¬ 
celler  depends  completely  on  the  proximity  of  the  closest  reference  to  the 
primary  to  achieve  cancellation.  Hence,  p  must  be  very  close  to  unity  to 

produce  cancellation  to  near  the  noise  floor.  Again  the  curves  in  Figures  2 

2  2  2 

through  6  approach  8/(l+8)  =  on  (uo)/(Oj  (u))+on  (w))  as  p  -*•  1  because 
the  model  maintains  the  ambient  noise  as  uncorrelated  from  hydrophone  to 
hydrophone,  regardless  of  p. 

G.  3  Comparison  With  Results  for  Sin  X/X  Cross-Spectral  Density 

Appendix  F  considered  the  line  array  structure  of  Figure  G-l  with  an 
extended  source  producing  the  hydrophone  CSD, 


2  2  sinl(p-q)  4>  sin  4>q) 

se(p,q,u>)  =  an  (a))  6pq  4-  (a.)  r  _ ,d 


(G-35) 


2  l(P-q)  sin  *01 


This  represents  a  uniform  source  with  angular  extent  ±i>  about  the  angle  $  , 
in  the  same  plane  as  the  line  array.  For  convenience,  let 

Y  =  —  sin  <t>  (G-36) 

'  y  o 

and 

W  =  $y  (G-37) 


G-12 


1 1  rui 


so  that  G-35  can  be  written 


se(p.q.w) 


cn2(w) 


6  +  _v 2(.) sin  niAszsm 

pq  I  7T(p-q)v 


(G-38) 


The  minimum  canceller  output  spectrum  for  this  CSD  has  been  determined 

explicitly  only  for  the  case  m  =  1.  In  the  noise  free  case,  using  the  results 

of  Widom  ^  ,  it  has  been  shown  that  in  the  absence  of  ambient  noise, 

2 

a  (u»)  =  0,  the  canceller  output  spectrum,  Ev(1,oj),  aporoaches  zero 
n  iv 

exponentially  with  the  number  of  references,  K,  if  W  <  1/2.  Then  the 
minimum  canceller  output  spectrum  is 


E  (l,u>)  =  0;  o  2 ( w )  =  0,  W  <  1/2  (G-39) 

oo  n 

Note  that  VV  =  y<&  is  analogous  to  the  product  B(w)d  in  the  exponential  case. 
The  result  given  in  (G-39)  is  a  significant  difference  between  the  sin  x/x 
CSD  and  the  exponential,  which  gave 

E  (l.ia)  =  of2(u>)  (1-p2)  ;  a  2(w)  =  0  (G-40) 

<*>  i  n 


regardless  of  B('jj)d.  This  difference  is  due  to  the  Markov  property 
underlying  the  exponential  CSD. 

In  the  noise  free  case,  c  2(w)  =  0,  when  j  <  W  <  for  n  =  1,  2.  .... 
the  minimum  canceller  output  spectrum  is  given  by 

9  ,  2W-n 

E  (l.u)  =  2* 4).'  (u>)  n(l+i-)  (G-41) 

°°  in 


G-18 


The  primary  hydrophone  output  spectrum  without  cancellation  is 
2$og^(jj)  so  that  the  relative  canceller  output  power  for  the  sin  x/x  CSD  is 

.  2W-n 

Er(  1 , u))  =  n(l  +  i  )  ( G-42) 

This  is  plotted  in  Figure  G-7  for  n=l  through  n=5,  and  shows  that  cancella¬ 
tion  is  severely  limited  if  W  >  1/2.  This  sampling  criterion  in  the  noise  free 
case  does  not  have  an  analog  in  the  case  of  the  exponential  CSD.  On  the 
other  hand,  when  B(io)d  >  1/2,  analogous  to  W  >  1/2,  then  p  <  .606  and  the 
cancellation  for  the  exponential  CSD  is 

ER(  1 , oj)  >  (1-p2)  =  .632 

so  that  significant  cancellation  is  not  achieved  in  the  exponential  case 
either.  The  degradation  is  much  more  graceful  in  the  exponential  case, 
however . 


W 


Figure  G-7.  Minimum  Canceller  Output  Power 
vs.  W.  for  Sin  x/x  CSD,  m  =  1 
Interference-to-Noise  Ratio  =  +<*> 


G-19 


Of  course,  the  case  when  noise  is  present  is  of  much  greater  interest. 

It  is  shown  in  Appendix  F  that  for  m=l,  with  noise  present  the  CSD  given  in 


(G-38)  yields 

2$o2(w)  2W-n 

E»(1»u)  =  — 2W —  S(1  +  S) 

for 


j<W<^i,n  =  0,  1,  2  ... 


and 


2W  a  “(w) 
S  =  n  +  .  _  ft  - 

2$aI2(w) 


2  2 

For  high  interference  to  noise  ratio,  24» Oj  (w)/on  (w)  >>  1, 


S  =  n,  n  =  l,  2,  ... 

so  that  for  n  >  1,  this  reduces  approximately  to  the  noise  free  case, 
For  W  <  1/2 


k2W 

2  ,  2$  a“  (u») 

EJl.u)  =  a  *(u)  ll  +  - - -  J  ,  W  <  1/2 

n  1  2Wo  *(u>) 

n 


When  W  =  0, 


Eb(1,u)  =  an  (u>) 


(G-43) 


(G-44) 


(G-45) 


(G-42) . 


(G-46) 


(G-47) 


G-20 


indicating  cancellation  to  the  ambient  noise  floor,  and  when  W  =  1/2 

Er ( 1, jj)  =  on2(oj)  +  2$  aj2(oj)  (G-48) 

so  that  no  cancellation  occurs. 

It  appears  that  the  parameter  B(ai)d  in  the  exponential  case  and 
W  =  y  $  sin  (Ji  in  the  sin  x/x  case  play  a  similar  role,  relating  the  distance 
between  hydrophones  (spatial  sampling  interval)  to  the  source  extent  ($) 
or  the  CSD  bandwidth,  B(w).  Therefore,  in  order  to  compare  the  behavior 
of  the  optimal  canceller  with  the  two  CSD's  the  irreducible  relative 
canceller  output  power,  Er(1,<jj)  was  computed  as  a  function  of  B(w)d 
and  W,  respectively,  using  (G-34)  and  (G-43).  Figures  G-8  through  G-ll 
show  Er  (l,tu)  for  the  exponential  CSD  for  interference  to  noise  ratios  of 
40  dB,  30  dB,  20  dB ,  and  10  dB  for  B(w)d  varying  between  0.0  and  2.0. 

Figures  F-5  through  F-8  show  ER  (l,w)  for  the  sin  x/x  CSD  for  the  same 
range  of  INR  with  W  varying  from  0.0  to  2.0.  The  following  observations  may 
be  made  regarding  the  figures,  letting  W  =  B(w)d  for  convenience; 

(a)  If  W  >  1/2  in  the  sin  x/x  case,  virtually  no  cancellation  is 

achievable,  as  already  observed.  No  such  sampling 

criterion  is  apparent  for  the  exponential  case,  but  cancellation 
is  still  severely  limited  for  B(co)d  >  1/2.  Either  CSD  requires 

a  value  of  W  much  smaller  than  1/2  if  the  canceller  is  to  produce 
cancellation  in  the  10  to  40  dB  range. 

(b)  The  optimal  canceller  produces  better  results  for  the  exponential 
CSD  for  W  larger  than  approximately  0.4.  However,  in  the 
important  region  of  small  W,  producing  more  than  10  dB  of 
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cancellation,  the  canceller  is  better  able  to  reject  the  sin  x/x 
source.  This  is  due  to  the  Markov  property  of  the  exponential 
source,  as  discussed  above. 


G .  4  Conclusions 

The  minimum  canceller  output  power  has  been  obtained  at  a  single 
frequency,  u),  for  a  source  producing  an  exponential  CSD  along  the  line 
array  of  references.  In  the  case  of  m=l,  where  results  exist  for  the  sin  x/x 
CSD,  the  exponential  CSD  produces  significantly  different  results  than  the 
sin  x/x  case,  as  noted  in  (a)  and  (b)  above.  The  differences  between  the 
two  reflect  a  Markov  property  inherent  in  the  double  exponential  CSD.  It  is 
likely  that  the  differences  also  result  from  the  fact  that  the  sin  x/x  CSD  is 
produced  by  an  angle  limited  source,  while  the  exponential  case  results 
from  a  source  distributed  on  (-tt.tt).  Like  the  sin  x/x  CSD,  the  source 
producing  the  exponential  CSD  cannot  be  cancelled  to  an  ambient  noise 
floor  for  m=l.  Further,  the  minimum  canceller  output  power  is  shown  to 
increase  with  m,  just  as  the  numerical  results  indicated  for  the  sin  x/x. 
However,  the  parametric  behavior  of  the  canceller  with  the  two  CDS's 
differs  sufficiently  that  the  exponential  case  cannot  be  used  to  predict  the 
canceller  performance  with  the  sin  x/x  CSD. 
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APPENDIX  H  -  SPATIAL  RESPONSE  OF  THE  LMS  CANCELLER  IN  THE 
PRESENCE  OF  AN  EXTENDED  SOURCE 

Although  the  primary  emphasis  of  this  study  has  been  on  the 
cancellation  of  an  extended  interference,  without  regard  for  the  effect  of 
the  cancellation  process  on  a  signal  of  interest,  the  eventual  application 
of  the  canceller  will  involve  detection  of  a  plane  wave  signal  that  is  masked  by 
the  interference.  This  Appendix  therefore  considers  the  plane  wave  response 
of  the  cancelled  output  as  a  function  of  the  angle  of  a  plane  wave  signal.  In 
other  words ,  given  a  primary  hydrophone  with  unit  weighting  and  K  refer¬ 
ence  hydrophones  with  the  frequency  domain  LMS  weighting 

"se  1('^)  sj/w)  (H-l) 

what  is  the  spatial  response,  or  beam  pattern,  of  the  canceller  output. 

The  spatial  response  at  frequency,  ui,  of  an  array  of  K+l  hydrophones 
with  arbitrary  weighting 

W(o!)  =  [W0(u)f  W2(u),  ...  Wk(u))]T  (H-2) 

with  the  weight  applied  to  hydrophone  k  at  frequency  w  is  given 

by 
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w  s 
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and 


=  propagation  delay  to  the  kth  hydrophone. 

For  the  spatial  cancellation  problem  with  the  array  geometry  of 
Figure  H-l, 


k 

cos  A 

o 

k=0,  1,  ...  K-l 

J 

COS<() 

o 

k=K 

that  is,  the  Kth  hydrophone  (the  primary)  takes  the  Jth  position, 
the  first  paragraph  in  this  section. 


W(ui)  = 


(w) 


The  vector,  d,  can  be  written  as 
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J  COS<l 
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Then  (H-3)  yields 


w = 


.  —  J  cos$ 

]  C  s 


(u.) 


(H-9) 


Now,  assume  that  the  extended  source  is  the  uniformly  distributed, 
uncorrelated,  narrow  source  defined  in  Appendix  A,  satisfying  (A- 26), 
(A-36),  and  (A-40).  It  is  shown  in  Appendix  E  that  for  this  source, 
Se_1(o))  can  be  written  in  terms  of  its  eigenvalues  and  eigenvectors  as 
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where  y  =  d/X  sim{>o  and 

Vk(K,W)  =  [Vo(k)(K,W),  V1(k)(K,W),  ...,  VK_1(k)(K,W)]T,  (H-ll) 

{V  ^k\K  ,W)}  is  the  kth  Discrete  Prolate  Spheroidal 

m  m-0,  1,  . . . K,-i 

Sequence  (DPSS)  index  limited  to  [0,  K-l],  and  Ak(k,w)  is  the  associated 
eigenvalue.  Also 
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Using  (H-10)  gives 
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using  (E-18)  and  (E-19)  from  Appendix  E.  Note  that 


K-l 

Vk+  (K , jj)  G  d  =  ^  Vm(k)(K,,;) 
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.  —  m  (cos$g  -  cos$0) 


(H-14) 


is  just  a  Fourier  transform  of  the  DPSS. 

Using  expressions  (H-13)  and  (H-14),  the  spatial  response  of  the 
canceller  in  the  presence  of  a  uniform  extended  source  has  been  evaluated 
for  various  values  of  source  extent,  2$  ,  source  angular  position,  <j>Q,  and 
interference  to  noise  ratio,  INK.  Figures  H-l  through  H-8  show  the  spatial 
response  when  m=l,  so  that  the  reference  array  plus  primary  comprises  a 
uniform  array,  with  the  following  parameters: 
d/X  =  0.5 
INR  =  40  dB 
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2j  =  interference  source  extent  =  20° 

=  interference  angle  =  45° 

The  figures  show  the  response  for  the  number  of  references,  K,  equal 
to  1,  2,  ....  8  respectively.  In  Figure  H-l,  with  only  one  reference,  the 
array  response  only  has  a  single  null  which  is  steered  to  the  center  of  the 
interference  source.  As  more  references  are  added  (in  Figures  H-2 
through  H-8),  the  notch  in  the  spatial  response  at  45°  broadens  to  approx¬ 
imately  20°  wide,  the  extent  of  the  source.  For  K  <  7,  this  broadening 
of  the  notch  is  accompanied  by  a  considerable  attenuation  of  the  response 
in  other  directions.  For  example,  with  K=2,  there  is  a  15  dB  attenuation 
at  endfire  (<t>=0°),  even  though  the  interference  is  45°  away.  Gradually, 
as  more  references  are  added,  the  notch  sharpens,  becoming  more  nearly 
rectangular,  and  the  attenuation  at  other  angles  away  from  the  interference 
decreases  to  a  more  acceptable  value. 

|  Note  that  once  K=3  or  4,  additional  references  do  not  appear  to  deepen 

the  notch  very  much,  suggesting  that  further  interference  cancellation  will 
not  be  substantial.  For  the  parameters  given  above, 

W  =  4  j  sin  <f>0  =  .062 

Figure  K-21,  which  shows  the  canceller  output  spectrum  as  a  function  of 
K  for  u)=.06  and  INR=40  dB ,  verifies  that  little  further  cancellation  occurs 
beyond  K=3.  Therefore,  Figures  H-l  through  H-8  indicate  that  the  reason 
for  adding  references  may  be  to  assure  adequate  signal  response  in  some 
direction  other  than  the  interference  arrival  angles.  For  example,  in  Fig¬ 
ure  H-3,  with  K=3,  a  signal  arriving  at  60°,  15°  away  from  the  interference. 
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suffers  25  dB  of  attenuation,  while  when  K=8  in  Figure  H-8,  the  signal 
attenuation  is  only  12  dB . 

When  interpretting  these  figures  it  should  be  kept  in  mind  that  the 
primary  consists  of  a  single  omnidirectional  hydrophone,  rather  than  an 
array,  as  would  be  used  in  most  sonar  situations.  As  indicated  in 
Appendix  L,  when  a  primary  array  is  used,  the  responses  like  those  of 
Figures  H-l  through  H-6  would  be  multiplied  by  the  response  pattern  of 
the  primary  array.  This  would  greatly  attenuate  (by  the  primary  array 
sidelobe  level)  those  regions  outside  the  mainlobe  of  the  primary  array. 
Assuming,  then,  that  the  reverberation  that  limits  the  active  detection 
performance  arrives  on  or  near  the  main  lobe,  the  signal  response  in  the 
vicinity  of  the  notch  is  of  primary  interest.  The  reason  for  adding 
references  would  then  be  to  increase  the  slope  of  the  sides  of  the  notch. 

It  is  observed  in  Appendix  K  that  as  the  primary  is  moved  farther 
from  the  references,  the  cancellation  of  the  extended  source  degrades. 
Figures  H-9,  H-10,  and  H-ll,  which  show  the  spatial  response  for  the 
same  case  as  H-l  through  H-8  but  with  m=5,  10,  and  20,  respectively 
confirm  this.  The  depth  of  the  notch  is  40  dB,  32  dB,  and  23  dB  for 
m=5,  10  and  20,  respectively  while  m=l  gives  a  notch  approximately  47  dB 
deep.  In  addition  to  the  reduction  in  cancellation  as  in  increases. 

Figures  H-9  through  H-ll  show  that  there  is  a  degradation  in  the 

array  response  to  signals  arriving  from  directions  other  than  the  angular 
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FIGURE  H- 1 4 .  Spatial  Response  for  20°  Uniform  Extended  Source  at  45 
d/X  =  .5,  m-1  ,  K=*l  ,  INR  =  20  dB 
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APPENDIX  I  -  USE  OF  THE  OUTPUT  OF  A  STEERED  LINE  ARRAY  AS  THE 
PRIMARY 

In  most  situations,  the  primary  input  to  the  optimal  canceller  is  not 
the  output  of  a  single  omnidirectional  hydrophone,  but  that  of  an  array  of 
hydrophones  steered  in  some  direction  of  interest.  The  directional  res¬ 
ponse  of  the  array  will  reduce  the  effect  of  reverberation  arriving  outside 
the  mainlobe  by  an  amount  equal  to  the  sidelobe  level,  but  reverberation 
arriving  in  the  mainlobe  will  remain  unattenuated  at  the  beam  output. 
Further,  situations  occur  in  the  sonar  environment  in  which  reverberation 
arriving  in  the  sidelobes  has  sufficient  power  to  prevent  detection  of  a 
much  weaker  signal  in  the  mainlobe. 

This  Appendix  considers  the  case  when  the  primary  input  to  the 
canceller  is  the  output  of  a  line  array  of  hydrophones,  colinear  with  the 
reference  array  as  shown  in  Figure  1-1,  and  steered  in  a  direction  <j>^,  as 
shown.  Let  the  N  hydrophones  in  the  primary  array  be  uniformly  spaced 
dp  feet  apart,  and  let  the  distance  between  the  primary  and  the  reference 
array  be  Ld  feet,  as  shown  in  Figure  1-1.  The  plane  wave  assumption 
made  in  Appendix  A  is  assumed  to  be  valid  over  both  the  reference  and 
primary  arrays. 

The  outputs  of  the  N  hydrophones  in  the  primary  array  are  denoted 
yn(t)  for  n=0,  1,  ...  N-l.  A  beamformer  steered  in  the  direction  <j>d 
computes  its  output  as 

N-i  . 

r — 1  d 

z(t)  =  2,  an  yn  (t  ~  fN~n_1)  cos  <)>d)  (1-1) 

n=0 


1-1 


Further,  the  nn(t)  are  uncorrelated  with  the  ambient  noise  components  of 
the  reference  hydrophone  outputs, 


E  [r^Ct)  nk(t+x)] 
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Since  z(t)  replaces  e^(t)  as  the  primary  input  to  the  canceller  struc¬ 
ture,  (A- 14)  becomes 
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Therefore  (A-15)  is  replaced  by 
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The  bracketted  term  is  just  the  response  of  the  primary  array,  steered  in 
direction  <|>d,  to  a  plane  wave  from  direction,  $  Denote 
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Figure  1-1.  CJse  of  a  Line  Array  as  the  Primary  Input 
where  the  an  are  real  shading  coefficients.  This  beamformer  output 
replaces  ed(t)  as  the  primary  input  to  the  canceller.  Extending  equations 
(A-6)  and  (A-7)  to  the  primary  array  gives  the  hydrophone  outputs  as 


YnU)  = 


nn(t) 


/'■ 


[t~£  [ (L+K-l)d  + 


nV 


cos  ]  d<j> 


(1-2) 


where  nn ( t )  is  the  ambient  noise  component  of  the  output  of  primary  hydro¬ 
phone  n.  As  before,  the  nn(t)  are  assumed  to  be  zero  mean,  stationary 
random  processes  with 
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so  that 


f 

/  <f>2,  w)B(cos  $1  -  cos  (t>d> 

"it  -n 

exp  {-j  [(L+K-l;d  cos  $  ^  -  kd  cos  <j>2)  ^d^ 2 
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Comparison  of  this  result  with  (A-15)  shows  that  the  case  with  a 
primary  array  can  be  treated  by  placing  a  single  primary  hydrophone  a 
distance  Ld  from  the  nearest  reference  hydrophone  if  the  hydrophone  has 
the  directional  response 

H(4>)  =  B(cos  —cos  $d)  exp  {-j  ^  (N-l)dp  cos  $d}  (1-9) 

This  equivalent  model  is  illustrated  in  Figure  1-2. 

Similarly,  it  can  be  seen  that  (A-16)  becomes 

N-l 

ru(T)  =  E  ( z ( t )  z(t+x)]  =  an2  Rn(r) 

n=0 

N-l  N-l 

*2  I  an*. 
n=0  m=0 

d 

+  -£■  m  (cos  <j>2  -  cos  <t>d) ,  <j>1,  <j,2}  d^d^  (I- 10) 


If 


RT  - 


-2-  n  (cos  d> . 
c  y  1 


cos  4>d) 


s,  (k,uo) 


/ 


1-5 


N=1 

2% 

n=0 


.IT  „TT 


a  (go)  + 
n 


/  /  W 


<t>2>  w)B  (cos  4^  -  cos  <j>d> 


-ir  -it 


B *(cos  <J>2  -  cos  4>d)  d^d^ 


(1-11) 


This  result  shows  that  in  the  equivalent  model  of  Figure  1-2,  the  spectral 
density  of  the  ambient  noise  component  of  the  primary  hydrophone  output 
must  be 
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Next,  consider  this  model  in  the  special  cases  discussed  in  Appdenix 
A.l.  Under  the  narrow  source  approximation, 
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Figure  1-2.  Equivalent  Model  for  Primary  Hydrophone  Array 


and  (A- 17)  is  replaced  by 
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Equations  (1-13)  and  (1-14)  correspond  to  (A-27)  and  (A-28),  respectively. 

When  the  source  is  additionally  assumed  to  be  spatially  uncorrelated 
(1-13)  and  (1-14)  further  simplify  to 
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In  Appendix  A.  2,  the  case  when  the  extended  source  is  uniformly 
distributed  was  considered.  Under  this  assumption,  the  CSD,  (k,  uj), 
and  the  PSD,  (a)  assumed  the  simple  forms  given  by  (A-40)  and 
(A-41).  When  an  array  is  used  at  the  primary  input  and  the  source  is 
uniform  as  defined  in  (A-38),  (1-13)  and  (1-16)  become 
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^  (K+L- 
c 


1-k)  cos 


'(u) 


(  d 

[B  (cos ( +  <j>)  -  cos  ij>d)  exp  {-j  (N-l)  cos  4>d  ] 

exp  {-j  —  (K+L-l-k)  <p  sin  4>Q}d<f  (1-17) 


and 


8U(.)  ' 


o  (w)  + 


t 


2.  . 

a.  (u) 


B  [cos(^r.  +  <(•)-  cos  <f> .  1 


d<J> 


(1-18) 


Hence  s^(k,u)  does  not  reduce  to  the  sin  x/x  form  given  in  (A-40)  unless 
B  [cos(<(|q  +  <fr)  -  cos  $^1  is  constant  for  4>e  [-$,  $].  That  is,  the  form  of 
(A-40)  only  occurs  if  the  directional  response  of  the  primary  array  is 
constant,  or  at  least  nearly  constant,  over  the  extent  of  the  interference 
source. 
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APPENDIX  J  -  THE  USE  OF  AN  ADAPTIVE  LMS  CANCELLER  IN  THE 
NON-STATIONARY  REVERBERATION  ENVIRONMENT 


In  practice,  the  second  order  statistics  of  the  reverberation  are  not 
known,  and  the  canceller  is  implemented  using  a  multiple  LMS  adaptive  can¬ 
celler  structure  as  discussed  in  the  Introduction.  However,  most  of  the 
work  in  this  report  has  considered  the  performance  of  an  optimal  (in  the 
least  mean  square  sense)  canceller  with  stationary  extended  source.  Since 
the  LMS  adaptive  filter  converges  in  the  mean  to  the  discrete  Wiener  filter, 
the  performance  of  the  optimal  filter  can  be  regarded  as  representative  of 
that  of  the  adaptive  canceller  if  algorithm  noise,  sampling,  and  finite  filter 
length  effects  are  negligible.  The  assumption  of  stationarity  is  necessary 
for  the  formulation  of  the  problem  in  terms  of  the  Wiener  (LMS  optimal) 
filter . 

In  a  passive  sonar  environment,  the  assumption  of  stationarity  on  the 
canceller  inputs  reasonably  models  most  sonar  problems.  On  the  other 
hand,  the  active  sonar  environment,  taken  as  a  whole,  is  markedly  non¬ 
stationary.  A  transmission  cycle  will  consist  of  a  transmission,  followed  by 
periods  in  which  reverberation  from  various  sources  (surface,  bottom,  con¬ 
vergence  zone,  etc.)  dominates  the  noise  field.  Between  reverberation 
periods,  the  ambient  noise  will  be  the  primary  limitation  on  sonar  perform¬ 
ance.  It  is  therefore  necessary  to  consider  the  effects  of  this  non- 
stationary  noise  field  on  the  performance  of  the  adaptive  canceller  and  any 
modifications  to  its  structure  necessary  to  accommodate  this  environment. 

The  canceller  structure  being  discussed  here  is  directed  specifically 
toward  suppression  of  surface  reverberation  from  the  convergence  zone. 

If  the  filter  is  adapting  during  periods  which  the  noise  is  due  to  the 


J-l 


ambient  field  or  to  reverberation  from  other  sources,  then  the  filter 
weights  may  not  be  reflective  of  the  convergence  zone  reverberation.  In 
order  to  allow  detection  of  targets  appearing  during  the  onset  of  surface 
reverberation  during  any  given  transmission  cycle,  the  canceller  must  have 
acquired  enough  information  about  the  reverberation  statistics  during  the 
preceding  cycles,  and  retained  enough  of  that  information  between  the 
reverberation  periods. 

It  is  well  known ^  that  the  LMS  adaptive  algorithm  adapts  much  more 
rapidly  in  an  environment  producing  correlation  between  its  inputs  (such 
as  reverberation  from  any  source)  than  to  one  producing  uncorrelated 
inputs  (such  as  ambient  noise  field  would).  That  is,  the  filter  "learns" 
the  correlation  properties  of  an  input  rapidly  in  comparison  to  the  way  it 
"forgets"  those  properties  once  the  correlation  disappears.  Therefore,  in 
order  to  maintain  the  properties  of  the  convergence  zone  reverberation  in 
the  filter  weights ,  it  is  first  essential  that  all  other  forms  of  reverberation 
be  excluded  from  the  adaptation  process.  Because  of  the  long  range  to  the 
convergence  zone,  this  can  be  done  by  freezing  the  adaptive  weights  fol¬ 
lowing  transmission  until  all  other  forms  of  reverberation  have  died  out, 
then  starting  the  adaptation  process.  Since  the  two  way  propagation  time 
to  the  convergence  zone  is  so  long,  the  adaptation  could  be  started  some 
fixed  time  interval  after  transmission. 

The  question  then  arises  as  to  the  effect  of  the  periods  in  which  only 
ambient  noise  is  present  on  the  filter  weights.  This  can  be  analyzed  for 
a  single  reference  structure  using  a  simplified,  approximate  model,  as 
follows.  The  LMS  algorithm  for  real  weights  is  given  by  [  3  ] 

W(N+1)  =  [I  -  u  X  (nlX^n)]  W(n)  +  yd(n)  X(n)  (J-l) 


J-2 


where 


W(n)  =  weight  vector  at  time  n 
X(n)  =  reference  hydrophone  output  data  vector 
=  [x(n),  x(n-l),...x(n-M)]T 
d(n)  =  primary  hydrophone  output  sequence 
u  =  adaptation  constant 

Using  the  assumption  that  the  present  weight  vector  is  uncorrelated  with 
the  present  data  vector  and  averaging  (J-l)  yields  a  difference  equation 
for  the  mean  weights, 

E  [W(n+l)l  =  [I  -  uRxx(n>]  E  [W(n)]  +vRdx(n)  (J-2) 

where 

Rxx(n)  =  E  [X(n)XT(n)] 

R^In)  =  E  [d(n)  X(n)] 

If  Rxx(n)  and  Rdx^n^  are  indePen<^ent  the  n>  then 

lim  E[W(n)]  =  R  'V,..  (J-3) 

n->“ 

which  is  the  discrete  Wiener  filter  as  indicated  above. 

Now,  assume  that  a  pulse  is  transmitted  starting  every  T  seconds,  and 

that  the  weights  are  frozen  during  the  period  nT  +  T  ,  that  is,  for  the 

w 

Tw  seconds  following  the  start  of  each  transmission  as  shown  in 
Figure  J-l.  This  presumably  excludes  all  reverberation  except  for 
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convergence  zone  from  the  filter  surface  reverberation  during  adaptation. 

From  the  point  of  view  of  the  adaptive  weights,  then,  the  input  consists  of 
the  intervals  [nT  +  T  ,  (n+l)T]  for  n  =  1,2,...  concatenated  together  as 
shown  in  Figure  J-l  to  form  a  continuous  input  sequence.  Suppose  that  this 
concatenated  input  is  sampled  every  T  seconds  to  produce  the  inputs  x(n) 
and  d(n) . 

First,  consider  the  primary  input  to  the  canceller.  Let  the  gate  func¬ 
tion,  G(n),  shown  in  Figure  J-l  characterize  the  period  during  which  the 
reverberation  is  present.  Let  dQ(n)  be  the  portion  of  the  desired  input 
due  to  ambient  noise  and  let  r^(n )  be  the  reverberation  as  seen  by  the 
primary.  Then 

d(n)  =  dQ(n)  +  G(n)rd(n)  (J-4) 

Similarly,  let  x  (n)  be  the  ambient  noise  at  the  reference  input  and  r  (n) 
be  the  reverberation  as  seen  by  the  reference.  It  is  now  assumed  that  the 
gate  function  seen  by  the  reference  is  just  a  delayed  version  of  G(n),  i.e., 

x(n)  =  xQ(n)  +  G(n-A)  rx(n)  (j-5) 

with  A  an  integer.  Then 

X(n)  =  X^n)  +  HG(n-A)rx(n)  (j-6) 

where 

Xq  =  [xQ(n),  xQ(n-l),  ...  xo(n-M)]T 

-x(n)  =  (rx(n)’  rx(n’1>-  •••»  rx(n-M))T 

J-5 


(J-7) 

(J-8) 


and 


HG(n-A)  =  diag  [Ct(n-m-A)] 
m  =  0, 1,  . .  .M 


Consequently 

Rxx(n)  =  PnX  +  HG(n'A)RRR  HG  (n‘A) 


(J-9) 


(J-10) 


and 


5dx(n)  =  HG(n'A)Rdx  G(n'A) 


(J-ll) 


where  it  has  been  assumed  that  the  ambient  noise  is  uncorrelated  in  time 
with  power,  Pn,  and  is  also  uncorrelated  between  hydrophones.  RRR  is 
the  reverberation  return  covariance  matrix,  and  Rdx  the  reverberation 
return  cross -correlation  vector.  Note  that  the  reverberation  has  been 
assumed  stationary  during  the  period  when  G(n)  is  unity. 

Now,  consider  HG(n-A)  for  o  n  <  k  as  shown  in  Figure  J-l. 

HG(n-A)  is  an  identity  matrix  if 

o^n-m  -  A  <_  p  ( J- 12) 


for  m  =  0 ,  1 ,  . . .  M ,  hence  ifM  +  A<_n_^p  + 


HG(n-A)  I>||Vl  +  A£nip  +  A 


(J- 13) 


Thus,  for  p-M  iterations  of  the  filter  during  the  first  p  iterations, 
HG(n-A)  =  I.  Similarly,  for  o  £  n  <.  k, 


HQ(n-A)  =  0 


(J- 14) 
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if  p<n-m-A<k  for  m  =  0,1,...  M.  This  occurs  if  p+M+A< n< k+A  ,  so  that 
HqOi-A)  =  0  for  k-p-M-2  iterations  during  the  k-p  iterations  from  n  =  p+1 
to  n  =  k-1.  These  arguments  can  easily  be  extended  to  the  successive 
transmission  cycles. 

It  now  is  assumed  that  the  length  of  the  filter,  M,  is  much  less  than 
the  reverberation  duration,  p,  and  the  duration  of  the  reverberation-free 
period,  k-p.  Then  "most  of  the  time",  G(n-A)  will  be  either  I  or  O,  and 
it  seems  reasonable  that  the  short  periods  when  this  is  not  true  will  not 
significantly  affect  the  filter  response.  Then,  approximately 

HQ(n-A)  =  G(n)I  ( J- 15) 

and 

R*x(n)  =  Po'  +  G(n)  RRR 

R^Cn)  =  G(n)  G(n-a)  Rdj(  (J-1J) 

If  it  is  additionally  assumed  that  the  duration  of  the  reverberation,  p,  is 
long  in  comparison  with  the  inter-hydrophone  delay,  A,  then  (J-ll)  can  be 
approximately  written 

R^x(n)  =  G(n)  Rdx  (J-18) 

without  significantly  affecting  the  results. 

Using  equations  (J-16)  and  (J-18)  in  (J-2)  yeilds 

E[W(n+l)]  =  {I(l-uPn)  -  uG(n)RRR)  E[W(n)]  +  PG(n)Rdx  (J-19) 
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Since  ROD  is  a  positive  definite  symmetric  real  matrix,  it  can  always 

KK 

be  written  as  normal  form  as 

Rdd  =  PT  A  p  =  P'1  A  p  (J-20) 

KK 

where  P  is  an  ortho  normal  matrix  consisting  of  the  eigenvectors  of  RRR’ 
and  A  =  Diag  [X  *2>  ...XM) 

=  eigenvalues  of  RRR 

Hence  (J-19)  can  be  written  as 

E[W(n+l)]  =  (l(l-yPn)  +mPT  G(n) AP}  E[W(n)]  +  PG(n)Rdx  (J-21) 

Multiplying  (J-21)  from  the  left  by  P  and  letting  Q(n)  =  P  E[W^(n)]  yields 
Q(n+1)  =  {(l-yPn)I  -  yG(n)A)Q(n)  +  uG(n)PRdx  (J-22) 

Since  the  matrices  inside  the  brackets  are  diagonal,  each  component  of 
Q(n)  can  be  solved  for  separately 

Q.(n+1)  =  1  -  p[PN+G(n)X.]}  Q.(n)  +  pG(n)(PRdx).  (J-23) 

The  above  is  a  piece-wise  constant  parameter  linear  system  with  on-off 
input  and  on-off  system  parameters.  Then  (J-23)  can  be  contrasted  with 
the  approach  in  [11]  where  the  eigenvalues  are  constant  but  the  eigenvec¬ 
tors  are  not.  Here  the  eigenvalues  of  (J-19)  (not  the  eigenvectors)  are 
the  functions  of  time.  Thus  transforming  to  the  eigenvectors  of  the  data 
covariance  matrix  allows  closed  form  solution  for  the  mean  weights. 
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Now  (J-23)  can  be  solved  in  a  piece-wise  manner  as  follows.  Letting 


(PR^x).  =  a.,  during  the  time  Of_n£p  (reverberation  return  is  present) 

n 

Q.(n)  =  Q.(o)  g!1  +  mi  ^  gj1'"1 

m=l 

where 

gi  =  1  '  y(Pn+Xi)  (J‘24) 


gj  is  the  discrete  time  impulse  response  of  the  system  described  by 
(J-23)  when  G(n)  =  1.  The  sum  can  be  written  in  closed  form  so  that 


Q^n)  =  Q.(o)  g*1 


+  jja. 


0  £n 


( J-25) 


Assume,  for  simplicity  that  Q.(o)  =  O  all  i.  Then,  for  p<  n<  k, 
G(n)  =  o  and 


Qj(n)  =  Q.(p)  bn  p ,  b  =  1  -  uPn. 


(J-26) 


That  is,  when  the  reverberation  disappears,  the  driving  function  of  (J-19) 
is  zero. 

The  system  decays  from  its  value  when  the  input  disappears  with 
impulse  response  bn.  After  one  complete  period  (n=k), 


Q.(k)  =  Q.(p)  bk~p 


(J-27) 
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If  one  repeats  the  same  arguments  during  the  successive  periods  of  G(n), 
the  following  solution  results 


uoUi-gi 
Q  [mk+p]  =  — r— - 
gi 

Q.[(m+l)k]  =  bkp 
at  the  switching  times. 
Q,(n>  =  bn-(mk+p) 


p.  m 

'  V  r  P  k-p.r 

-  2,  s  1 

r=o 

Q.[mk+p] 

In  between  switching  times, 

Q.[mk+p] 


(J-28) 


for  mk+p  <  n<  (m+l)k 

(J-29) 


and  for  mk  <  n  <  mk  +  p 

,.  n-mk. 

QjCn)  =  gj11'^  Qjlmk)  +  ua.  ^ -  (J-30) 

Equations  (J-29)  and  (J-30)  have  simple  physical  interpretations. 
Equation  (J-29)  implies  that,  between  switching  times  when  no  input  is 
applied,  the  system  decays  with  a  rate  b  from  its  value  when  the  input 
was  last  applied,  while  (J-30)  implies  that  the  system  response  consists  of 
(a)  a  transient  response  term  g.n  mk  Q^(mk)  due  to  the  initial  condition 
(Q.(mk))  when  the  input  is  turned  on,  with  decay  rate  g^,  and  (b)  a 
response  to  the  input  ucc.  with  the  step  response  of  first  order  system  with 
time  constant  g.,  ( 1- gjn  mk)/(l-g.). 


As  the  number  of  iterations  increases,  the  Qj(n)  becomes  periodic 


with 


lira  Q  (n)  =  bn'(mk+p)  Qi(mk"p)  mk+pmi(m+l)k  (J-31) 

m-*-® 

ua. 

lim  Q.(n)  -  y-T 
m-*-®  6i 

As  a  check,  letting  m-*-®  in  (J-28)  and  n=mk+p  in  ( J- 32)  yields  equality. 
Also,  letting  n=(m+l)  and  m-*-®,  (J-31)  and  (J-28)  yield  agreement.  Also 

pa. 

lim  Q.(n)  =  y— p ,  P  =  k  (no  dead  time) 
m-*-®  1  ®i 

lim  Q.(n)  =  0  ,  p  =  o  (no  input)  (J-33) 

m-*-® 

This  result  is  analogous  to  the  response  of  an  RC  filter  to  a  periodic 
rectangular  pulse  train  as  shown  in  Figure  J-2.  There  is  an  initial  tran¬ 
sient  buildup  that  tries  approach  a  steady  value  that  would  occur  with  no 
pulsing.  Because  of  the  pulsing,  a  decay  is  superimposed  upon  this 
transient.  If  the  response,  when  the  pulse  is  on,  is  large  in  comparison 
to  the  response  when  the  pulse  is  off,  the  steady  state  value  is  quickly 
achieved.  If  the  reverse  is  true,  then  the  weig!  ts  do  not  change  signifi¬ 
cantly.  Two  factors  contribute  to  the  response;  the  percentage  of  time  in 
each  mode  (duty  cycle)  and  the  rapidity  of  the  response  in  each  mode. 

Now  the  rapidity  of  the  response  in  each  mode  is  given  by  gj  and  b  with 

g.  =  1  -  u(P  ♦*  )  <  1  -  upn  =  b  (J-34) 


g™(l-OgK-p 


- i -  +  (l-g. 

l.g.Pbk‘P 


n-mk 


mk<n<mk+p  (J-32) 


J-ll 
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Figure  J-2.  Weight  Response  Analogy  to  RC  Filter  Response 
Using  (J-28),  as 


ua 


i  > 


Qitmk+p]  =  (1.gj)  d.g.Pt.K-P) 

N 

Equation  (J-35)  can  be  viewed  as  consisting  of  three  factors, 

ua , 

steady-state  response  with  no  pulsing  =  j-  -1 

®i 

transient  response  during  one  pulse  =  1-g^ 

k-  ~1 

long  term  transient  response  =  (l-g.PbK 


(J-35) 


(See  Eq.  (J-33)) 
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Imbedded  in  the  last  factor  is  the  time  constant  to  achieve  a  periodic 
steady-state.  From  (J-28),  the  factor  in  the  sum  can  be  written  as 
follows : 


m  km 

(g?  bk~p^  =  ^gp/k  b1_k^  (J-36 

However,  km  =  total  no.  of  algorithm  iterations  (ignoring  relatively  few 
iterations  during  one  period  for  large  k).  Hence,  the  factor  gjP^k  b1  p^k 
is  the  effective  time  constant  of  the  pulsed  system.  Now 


g.P/k  bl-p/k  = 


l-U(Pn+A.) 


,  P 
k 


1-uP 


a-.pn) 


1  - 


U  A. 


1-uP 


’I 


n 


(l-PPn) 


(J-37) 


But  p/k  is  the  duty  cycle.  If  ypn  «1  and  -X.  «  1,  all  i,  then 


gp/k  b1_p/k  «  (1-yE  X.)  (l-upn)  *  1  -  u(Pn+Xi  £>  (J-38) 

Hence  the  response  of  the  pulsed  system  is  determined  by  the  product  of 
u  and  the  sum  of  the  noise  power  plus  the  eigenvalue  times  the  duty  cycle 
(sort  of  the  time  average  eigenvalue).  This  result  is  very  physically 
satisfying. 

Alternatively,  in  equation  (J-37),  the  factor  ( l-g.p) /( l-g.p  bkp)  can 
be  viewed  as  the  steady- state  amplitude  loss. 
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One  may  now  transform  back  to  the  original  coordinate  system  to 
obtain  the  behavior  of  E[W(n)].  Assuming  steady-state  periodic  behavior 
and  using  equations  ( J- 31)  and  (J-32),  for  mk+p  <  n<  (m+l)k 


lim  E[W(n)]  =  bn'(mk+P>  £[W(mk+p)] 
n-*-  * 


_  bn-(mk+p) 


T 

uP  Diag 
i=0, 1, . . . 

M-l 


l-giP  bk'P 


-dx 


(J-39) 


and  for  mk^n^mk+p 


lim  E[W(n)] 

n-*-°° 


^PT  Diag 


n-mk 


1-gj P 

1_gi  1-gP  bk‘P 


,k-p  + 


i-gi 


n-mk  1 


i-gi 


PR 


dx 

(J-40J 


It  is  clear  from  equations  (J-39)  and  (J-40)  that  the  mean  weights  are 
time  varying,  even  in  steady-state.  The  mean  weights  are  periodic,  how¬ 
ever,  to  determine  the  difference  between  the  mean  weight  and  the  Wiener 
weights  requires  specific  knowledge  of  P.  Note,  as  a  check,  for  p=k, 
(J-40)  becomes 


lim  E[W(n)]  =  uP1  Diag 


n-+oo 


i-e,- 


PR  . 
-dx 


T 

uP  Diag 


u(  X  .+P  ) 
i  n' 


PR  , 

-  dx 


-1 


=  Rxx  Rcjx  =  Wiener  weights 


\ 


(J-41) 
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Of  primary  interest  are  the  weights  during  the  time  the  reverberation 


is  present,  given  by  (J-39).  Note  that  (1-g.)  =  1  -  [  1- u(Pn+  > .)  ]  = 
u(Pn+  X.),  so  (J-39)  can  be  written 


lim  E[W(n)]  =  PT  diag 
n-<-«>  i 


b 


k-p 


n-mk  ,  k-p  ,  , ,  n-mk 
b  *  +  (1-g. 


diag 

i 


P  +X. 
n  1 


PR  . 
-dx 


PT  tPn+  A]'1  PBdx 


+  PT  diag 


„  n-mk] 
gj  diag 

-*  i 


l-b1^ 


1-g.P  bk'P 


-1 


PfPn+  P*dx 


5dx  -  (1-bk'P)  pT  Din'mk  D2  PIRxx  ?dx>  (J-42) 


where 

Dx  =  diag  [g.]  =  diag  [l-v(P  +*.)] 
i  i 


(J-43) 


and 


D„  =  diag 


i-Bip  bk-e 


diag 

i 


l-tl-KPn+^)]P  [l-uPn]k'p 


( J- 44) 
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where  R  is  the  reference  covariance  that  would  occur  if  the 
xx 

reverberation  were  always  present. 

Then  the  first  term  is  just  the  Wiener  weights  that  would  occur  if  the 
reverberation  were  always  present,  while  the  second  term  is  a  time  varying 
deviation  of  the  weights  due  to  the  pulsed  nature  of  the  reverberation. 

Note  that  since  m  simply  defines  the  particular  transmit  cycle  in  which  n 
lies,  ( J - 42)  can  be  rewritten 

lim  E[W(n)]  =  Rxx  "  ‘5dx  *  <l-bk-P)PT  Dj"  Dj  P(RXX  (J-45) 

n-*  °° 

with  n'  the  time  since  the  onset  of  reverberation  in  the  current  cycle. 

At  the  beginning  of  the  reverberation  (n'=0) 

lim  E[W(p)l  =  Rxx  ‘‘r  .  (l-bk’P)PT  D2  P(Rxx  ‘‘r^I  (J-46) 

n-*-* 

while  at  the  end  of  the  reverberation  (n'=p) 

lim  E[W(n) )  =  Rxx  .  <l-bk-P)PT  D1PD1  P(Rxx  " ^Rdx>  (J-47) 

n-*-°» 

The  conditions  under  which  the  weights  converge  to  Rxx  *Rdx  are  that 
either 

(l-bk'p)  =  0  =>k  =  p  (J=48) 

or 

D*'  D2  =  0  (J-49) 


J-16 


But 


D^n  D2  =  diag 
i 

For  a  filter  with  acceptable  algorithm  noise,  [Pn+Aj]«  so 

D^1  D0  >  0  ( J - 51) 

Hence  the  Wiener  weights,  Rxx  Rdx>  are  not  reached  in  steady  state 

unless  p=k,  so  that  the  reverberation  is  always  present.  This  means  that 

the  canceller  performance  is  degraded  fromthat  which  would  be  achieved  if 

reverberation  were  alwyas  present  at  the  filter  input.  The  degradation 

K“D 

will  depend  upon  the  factor  (1-b  *)  and  upon  the  time  varying  maxtnx, 

Dj11  D^,  given  by  (J-50). 

Of  most  interest  here  is  the  mean  square  error  of  the  adaptive  can¬ 
celler,  which  represents  the  power  in  the  output.  This  is  given  by 

E[E2(n)]  =  E[d(n)-WT(n)  X(n)]2 

=  E[d2(n)]  -  2E[WT(n)  d(n)  X(n)] 

+  E[WT(n)  XT(n)  W(n)]  ( J- 52) 


['•"iv'i']  ” 

'  P  r  -I  k-p 

1  -  [H[Pn+,i,J  t1_UPn. 


(J-50) 


J-17 


r 


Just  as  in  the  derivation  of  the  mean  weights,  (J-2),  it  is  quite 
common  to  assume  that  the  weights  are  uncorrelated  with  the  present 
data,  so  that 

Ew,  A [E2(n) ]  =  E[d2(n)]  -  2W(n)  R,  (n)  +  WT(n)  R  (n)  W(n)  (J-53) 
wen)  —  —ax  —  xx  — 

where  E^n^  denotes  the  expectation  conditioned  on  W(n).  In  the  actual 
adaptive  implementation,  the  weights  are  random  variable  which  can  be 
written  as 

W(n)  =  E[W(n)]  +  Qw(n)  (J-54) 

with  Qw(n)  the  zero  mean,  fluctuating  part  of  W(n).  Ignoring  these 
fluctuations  is  equivalent  to  neglecting  algorithm  noise,  which  reduces 
(J-53)  to 

E[E2(n)]  =  E[d2(n)]  -  2E[W(n)]  Rdx(n) 

+  E[WT(n)]  Rxx(n)  E[W(n) ]  (J-55) 

If  the  reverberation  were  always  present,  then  the  inputs  are 
stationary,  and  (  )  yields  the  Wiener  solution  in  steady  state, 

lim  E[E2(n>  =  R  d  -  R*,  Rd  (J-56) 

n+  = 

2 

with  Rdd  =  E[d  (n)]  under  the  condition  that  the  reverberation  is  always 


present.  However  in  the  actual  case,  it  has  been  shown  above  that  the 
mean  weights  are  periodic  as  n-*-=°,  and  are  given  by  equation  (J-42). 


This  can  be  rewritten 


lim  E[W(n)]  =  [I-A]  R 
n->“ 


-1 

xx 


( J- 57) 


where 


=  M-hk~P)pT  d  n  mk  D„P 
1 


A  =  ( 1-b 


( J- 58) 


The  matrix,  A,  displays  the  time  varying  nature  of  the  man  weights. 

Note  that  when  k=p,  so  that  reverberation  is  on  all  the  time,  A=0,  and 
equation  (J-57)  reduces  to  (J-56).  On  the  other  hand,  if  p=o,  so  that 
reverberation  is  always  present,  A=l,  and  the  mean  weights  are  zero,  as 
they  should  be  in  an  uncorrelated  environment.  Substituting  (J-57)  into 
(J-55)  under  the  condition  that  the  reverberation  is  present  gives 


lim  [E2(n)J  =  Rdd 
n-*-00 


-  T  -  -1 
5dx  Rxx 


5dx  +  ?dx 


XX 


A«*x  Akxi  5dx  <J-59> 


The  first  two  terms  are  just  the  mean  square  error  that  would  be  achieved 
if  the  reverberation  were  always  present,  from  (J-56),  while  the  second 
term  is  the  increase  due  to  the  "pulsed"  character  of  the  inputs. 

Now,  denote  the  second  term  in  (J-59)  as  y  and  use  (J-58  to  write 

Y  =  (l-bk~P)2  Rdx  PT  D2(n‘mk)  D2  (Pn+A)'1  P  Rdx  (J-60 

t  Vi  - 

Let  a.  be  the  i—  element  of  p5dx’  projection  of  Rdx  on  the 

i—  eigenvector.  Then  (J-60)  becomes 


J-19 


(J-61) 


-  a-bk'P)2  £ 


2  D  2(n~mk)  d2 
K-l  af  li  2i 


i=0 


Pn  +  9 

n  2, 


where  and  D^.  are  the  i—  diagonal  components  of  and  Dg, 

IK  __ 

respectively,  and  where  ^  is  the  i—  eigenvalue  of  Rxx-  This  can  be 
written  as 


K-l 

Y=(l-bk-P)2  £ 


ctj[l-u(pn+  X.)]2(n_mk) 


i=0  (Pn+Xi)  [1  '  d-u(Pn-x.))P  d-uP„)k’P]2 


(J-62) 


n 


In  order  to  evaluate  the  increase  in  he  canceller  output  power  due  to  the 
"pulsed"  character  of  the  reverberation,  it  is  necessary  to  know  the 
eigenvalues  of  Rxx  and  the  projections  of  Rrfx  on  the  eigenvectors  of  Rxx- 
Further,  the  increase  in  output  power  can  be  seen  to  depend  upon  y,  but 
this  dependence  is  quite  complex,  as  (J-62)  shows. 

J.  1  Special  Case  -  Narrowband  Interference 

Suppose  that  the  interference  consists  of  a  single  frequency  sine  wave 
of  frequency,  wo,  and  that  the  array  has  the  geometry  of  Figure  A-l.  Then 

kxx  ( ^1  d_  <1+  •+•  d*  d^)  (J-63) 

2 

where 


[I,e',uo(At), 


e-jwo2(M)i 


e-jwQ(K-l)(At)]T 


(J-64) 


J-20 


+  T 

where  d  d  =  K,  and  for  simplicity  d  d  =  0.  Also 

At  =  spacing  between  filter  taps 

and 

=  PI  ,  ~  j cli  A  ,  ,  jo.  A  , j,. 

-ux  —  (e  •>  o  d  +  eJ  o  d*) 

This  special  case  gives 


(J-65) 


( J-66) 


P  =  [(d  +  d*),  a*,  a*,  ...  a*R_2JT 


( J-67) 


where  the  a.  are  any  set  of  orthonormal  vectors  or  thogonal  to  d  +  d*. 
Then 


ai~ 


o. 


x  Pr  Cos  od  A  ,  i  =  1 

L  1  0 

otherwise 


Substituting  (J-68)  into  (J-62)  gives 


(5-68) 


[1  -  ( l-uPn)k_P] 2  MPj2  [1  -  u(Pn+MPj)]2(n'mk) 
(Pn+MPj)  [1  -  (l-y(Pn+MPI))P  (l-uPn)k‘p] 


(J-69) 


J-21 


J-22 


APPENDIX  K:  NUMERICAL  EVALUATIONS  OF  THE  LMS 
CANCELLER  PERFORMANCE  WITH  A  UNIFORMLY  DISTRIBUTED, 
SPATIALLY  UNCORRELATED  NARROW  SOURCE 


I 

I 
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APPENDIX  K  -  NUMERICAL  EVALUATIONS  OF  THE  LMS  CANCELLER 
PERFORMANCE  WITH  A  UNIFORMLY  DISTRIBUTED,  SPATIALLY  UNCOR- 
RELATED  NARROW  SOURCE 

Using  the  results  of  Appendix  E,  the  canceller  output  spectrum  for 

K  reference  hydrophones,  E^(<j),  can  be  evaluated  on  the  computer. 

Similarly,  using  the  results  of  Appendix  H,  the  spatial  response  of  the 

optimal  canceller  can  be  determined.  Evaluation  of  (E-20)  for  Ek(^>  and 

(H-13)  for  the  spatial  response,  requires  determination  of  the 

(k) 

DPSS,  tVj  i  K_^  for  JilK,  and  the  associated  eigenvalues, 

xlc(K,W).  The  DPSS,  { Vn^(K ,W)^ k_Q  ^  for  0£niK-l  are  computed 

by  numerically  solving  the  eigenvalue  problem 

S(K,u>)  Vk(K,W)  =  *k(K,W)  Vk(K,W)  (K-l) 


where 


SCK.O  =  (K-2 

V  ”<p-q)  >  p=o,  1,  ...  K-l 

q=0,  1,  ...  K-l 

The  resulting  Vk(K,W)  are  given  by 

Vk(K,W)  =  [V0(k)(K,W),  V^^CK.W),  ...,  VK1(k)(K,W)]T 

and  the  Ak(K,W)  are  the  desired  eigenvalues.  The  v/k^(K,W)  for  JiK 
can  then  be  determined  by  evaluating  equation  (E-10), 


K-l 

I 


sin  [27rW(n-m>] 
ir(n-m) 


(k)(K,W) 


(K,W)  =  Ak(K,W)  Vn 


(K-4) 


In  the  development  of  the  narrow  source  model  in  Appendix  A,  it  was 
assumed  that  $  was  sufficiently  that 

sin  $  $ 

(K-5) 

COS  a  1 

where 

W  =  <t>  y  sin  (K-6) 

This  restriction  will  be  interpreted  as  requiring  that  $*.2  radians 
(approximately  11. 5°) ,  mainly  to  assure  that  cos  4>=*1.  Therefore 

W  i  .  2  j  sin  'fg  <  .2  j  (K-7) 

Usually  in  array  designs, 

.5  <j  £1  (K-8) 

If  d IX  <  .5,  the  ambient  noise  components  of  the  hydrophone  outputs  will 
be  correlated,  limiting  their  usefulness  for  increasing  array  gain.  If  the 
upper  limit  is  exceeded,  then  the  array  response  will  include  grating  lobes. 
These  considerations  suggest  that  the  main  interest  will  be  in  values  of 
W  <  .2.  Certainly,  given  the  sampling  criterion,  W  £  .5,  developed  in 
Appendix  F,  values  of  W  ^.5  will  be  of  no  practical  interest.  Several 
cases  with  W  >  .5  will  be  shown,  however,  to  illustrate  the  importance  of 
the  sampling  criterion. 


K-2 


The  plots  of  the  numerical  results  included  here  are  summarized  in 

(W) 

Table  K-l.  Each  plot  shows  the  canceller  output  spectrum,  ,  at  a 

single  frequency,  W,  as  a  function  of  the  number  of  reference  hydrophones 
used  in  the  canceller  for  a  particular  value  of  W.  The  results  are  given 
for  a  range  of  values  of  the  parameter  m,  where  the  reference  array  is 
a  distance,  md,  from  the  primary  hydrophone.  To  determine  the  cancella¬ 
tion  for  a  particular  case  of  interest,  the  hydrophone  spacing  is  computed 
in  wavelengths  for  the  frequency  of  interest,  Since  cancellation  will 
improve  as  W  decreases,  the  largest  value  of  ^  (the  angle  to  the  inter¬ 
ference,  always  £  n/2)  should  be  vised  to  compute  (K-6),  so  that  the 
resulting  cancellation  performance  will  be  an  lower  bound  on  achievable 
cancellation.  The  distance  from  the  primary  to  the  references  is  then 
determined  in  units  of  d  to  provide  the  parameter,  m. 

Figures  K-l  through  K-5  were  plotted  using  values  of  W  that  violate 
the  sampling  criterion,  W  >  1/2,  to  further  illustrate  just  how  limited 
cancellation  is  if  this  condition  is  not  met.  These  curves  clearly  support 
Figure  F-3,  which  showed  the  irreducible  canceller  output  spectrum,  E^il.w) 
for  the  case  m=l.  Since  it  was  also  shown  in  Appendix  D  that 

f 

E^fm.co)  >  E^(l,aj),m  >  1  (K-9) 

it  would  be  expected  that  cancellation  is  also  severely  limited  for  mM  when 

W  >.5. 

The  remaining  Figures  illustrate  the  properties  of  the  canceller 
structure  derived  in  the  preceding  sections.  In  the  presence  of  noise, 
cancellation  to  the  noise  floor  is  never  achieved,  as  predicted,  and  in  fact 


the  output  spectrum  can  be  significantly  above  the  noise  floor.  Appendix  D 
derived  the  cancellation  floor,  E^Cl,*),  for  the  m=l  case  with  W<l/2  as 


E  (1,*)  =  2$oIV)(lNR)  1  [1  +  ^p)2W  (K-10) 

9 

where  24a^  («)  is  the  interference  power  received  at  any  of  the  hydro¬ 
phones  and  where  INR  is  the  interference  to  ambient  noise  ratio, 


INR 


2<*>0j2(u)) 

a  (w) 
n  v  ’ 


(K-ll) 


2 

The  power  at  the  primary  output  without  use  of  the  canceller  is  2 4^  0-0  + 

9 

cn  (w),  so  that  the  maximum  achievable  cancellation,  0^(1,^),  is 


C(l,w)  = 


HU  1,0 


2<5>Oj2(u)  )  +  on2(u) 


1  +  INR 


fl  +  *M]2W 
11  2W  1 


(K-12) 


This  has  been  computed  and  shown  on  Figures  K-6  through  K-62. 

Figures  K-6  through  K-62  clearly  illustrate  the  diminishing  returns  in 
cancellation  performance  as  references  are  added.  If  the  cancellation  floor 
is  significantly  below  the  uncancelled  primary  output,  that  is,  if 

E“< w)| m=i  =  E~(1’w)  <c  2$0!2^>  +  °n2(U))  (K-13) 


then  the  performance  improves  dramatically  with  the  first  few  references. 
However,  the  figures  show  a  well  defined  transition  point  in  K  beyond 
which  the  additional  interference  rejection  afforded  by  adding  another 


K-5 


Figure  K-2;  Canceller  Output  Power 
vs.  No.  of  References  when  Sampling 
Criterion  is  Exceeded 
No  Ambient  Noise 
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Figure  K-3;  Canceller  Output  Power 
vs.  No.  of  References  when  Sampling 
Criterion  is  Exceeded 
No  Ambient  Noise 
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reference  is  quite  modest.  For  example,  with  INR=40  dB  and  VV = 1 0 
(Figure  K-30),  the  first  reference  gives  28  dB  of  rejection,  the  second 
adds  6  dB  more,  but  the  third  gives  only  about  2  dB  additional  cancella¬ 
tion.  The  hardware  required  to  incorporate  each  new  reference  into  the 
structure  is  the  same,  and  includes  the  computational  load  of  computing  the 
weights  in  an  adaptive  implementation.  Further,  in  a  adaptive  canceller, 
the  adaptive  loop  associated  with  each  reference  adds  algorithm  noise  to 
the  system.  This  algorithm  noise  can  actually  offset  the  improvement 
due  to  cancellation  improvement.  The  choice  of  the  number  of  references 
must  trade  off  the  rejection  added  by  a  reference  against  the  additional 
algorithm  noise  and  computations. 

In  Appendix  F,  it  was  shown  that  the  canceller  output  spectrum  for  K 

references  and  m=l,  that  is,  Ev(l,tu)  in  the  noise  free  case  (o2  (w)=0) 

x\  n 

is  asymptotic  in  K  to 


A(K,W) 


2<foI2 

2VT 


(^) 


[sin  w  WJ 


2K+1 


(K-14) 


Figure  K-63  shows  EK(1,W)  for  various  W  plotted  along  with  A(K,W)  for 
2 

an  ( i»0=0.  It  can  be  seen  that  A(K,W)  is  virtually  indistinguishable  from 
2 

Ev(l,<*»)  when  o  (<*))=0  even  for  small  K,  so  that 

IS.  II 


Ek(l,<*0 


a 

n 


2(<*0=0 


2$aI2(u)) 

~ 2W - 


[sin  w  W]2K+1 


(K-15) 
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When  the  interference-to-noise  ratio  is  high,  INR">1,  then  (K-17)  reduces 
approximately  to 

^  1/2  |  [log  (sin  ttW)]-1  [(2W-1)  log  (^)]-l  }  ,  W<l/2  (K-18) 

Since  the  selected  number  of  references  must  be  an  integer,  and  taking 
into  account  that  the  intersection  given  by  (K-18)  is  slightly  lower  than 
the  transition  point  in  the  figures,  let 

KQ  =  fKjl  (K-19) 

where  TXl  is  the  smallest  integer  greater  than  or  equal  to  x.  Choosing 
the  number  of  references,  K,  equal  to  Kq  for  m=l  does  not  mean  that  fur¬ 
ther  cancellation  is  not  possible,  but  that  further  cancellation  is  small  relative 
to  the  increase  in  computations  and  algorithm  noise.  As  an  example,  in  Fig¬ 
ure  K-18  with  INR=40  dB  and  W  =  .08,  (K-19)  gives  Kq=3  and  K=3  produces 
27  dB  of  cancellation.  Use  of  five  additional  references  (K=8)  provides 
only  2.5  dB  more  rejection  of  the  interference,  which  would  probably  not 
justify  the  additional  computational  cost  or  make  up  for  additional  algorithm 
noise  of  using  five  more  references. 

The  number,  Kj,  provides  insight  as  to  how  the  number  of  references 
required  in  a  given  situation  changes  with  W  and  with  INR.  From  (K-18) 
it  can  be  seen  that  Kj  is  linear  in  10  log  (INR),  which  is  the  interference- 
to-noise  ratio  in  dB.  The  slope  of  Kj  with  respect  to  10  log  (INR)  is 

^  ^1  1  2W-1 

d(  10  log  (INR))  =  20  log  (sin  ^W)  (K- 20) 


K-  71 


so  that  approximately  (2W-l)/(20  log(sin  ”W)  references  are  required  per 

dB  of  INR.  Figure  K-64  shows  Kj  plotted  as  a  function  of  10  log  (INR) 

for  various  W .  The  dependence  on  the  interference- to-noise  ratio  is  quite 

weak  when  W1.01,  that  is,  when  the  interferences  are  very  narrow. 

The  dependence  of  Kj  or  W  is  more  complicated  as  shown  in  Figure 

K-65,  which  shows  Kj  plotted  as  a  function  of  W  for  INR=20,  30,  and  40  dB. 

-2 

For  small  W  (say,  W<10  )  the  dependence  of  on  W  is  very  weak, 

-2 

regardless  of  INR.  As  W  increases  above  10  ,  the  value  of  Kj  begins 

to  increase  rapidly  with  W,  going  asymptotically  to  infinity  at  W=l/2. 
Recalling  from  Section  3. 3. 1.1  that  no  cancellation  is  possible  when 
W=l/2,  it  would  be  expected  that  Kj'*®  at  this  point. 

Now,  suppose  the  criterion  K=KQ  is  used  to  select  the  number  of 
references  used  in  a  given  situation.  Figures  K-66  shows  that  the  can¬ 
cellation  achieved  is  within  several  dB  of  the  cancellation  floor,  (K-12)  for 


this  choice  of  K. 


In  fact,  it  can  be  seen  that  the  value  of  E^gll.u))  is 


very  nearby  log  (INR)  dB  above  the  floor,  so  that 


ER  (l.u)  «  (INR)’1  B.d,*)  =  2ioj2(a))(INR)  '9  [1  +  2W  (K-21) 

provided  that  ( INR )’*"  E„(l,w)  ±  +  o  20").  This  condition  just 

assures  that  the  approximation  is  not  used  when  it  would  produce  an 
apparent  increase  in  the  interference  power  relative  to  the  uncancelled 
primary  output.  If  INR  ?>1  and  W<1,  as  it  will  be  in  all  cases  of  interest, 
(K-21)  can  be  approximated  as 


Ek-  (l,u>)  *  2*o  2(<u)  (INR)2W_’9  ( 2W) 
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Figure  K-65.  Change  in  K.  as  a  Function  of  W 


-21V 

For  0  <  VV  <  1 ,  ( 2W )  is  within  about  1.6  dB  of  unity ,  so 


Ew  (1,*)  =  2*o.Z(u>)  (INR)2W  *9 
K0  1 


(K-22) 


The  cancellation  achieved  is  then  approximately 


CK0(1>a))  “ 


E„  (l.u) 
K0 


on"(u)  +  24oi20)  1  +  INR 


2W+  1 

(HLR)_Z— -1  S  (INR)2W--9  (k-23) 


for  INR>>1.  Expressed  in  dB ,  this  yields 


10  log  CK  (l,w)  =  -.9  [10  log  (INR)]  +  2W  [10  log  INR]  (K-24) 

0 


The  interference  rejection  in  dB  varies  approximately  linearly  with  both 
the  interference  to  noise  ratio  in  dB  and  the  parameter  W.  This  result  is 
valid  over  the  range  10  4iWf..45  and  102£INR£104. 

Further  inspection  of  Figures  K-6  to  K-62  shows  that  for  any  given 
values  of  W  and  INR,  the  curves  of  Eg(m,u)  have  roughly  the  same  shape 
as  that  of  EK  ( 1 .  <*>) ,  and ,  in  particular ,  have  approximately  the  same 
transition  point  in  K.  This  suggessts  that  Kg  may  be  a  suitable  guide  for 
the  selection  of  the  number  of  references,  K,  even  when  m>l.  It  can  be 
seen  from  the  figures  that  for  m»l,  E^Cm.oo)  continues  to  fall  off  some¬ 
what  faster  than  E^(1,uj)  for  K>Kq.  For  example,  in  Figure  K-31,  with 
W  =  .01  and  INR  =  103,  KQ  =  2  which  yields  10  log  C2(l,u>)  =  -25.5  dB. 

Six  additional  references  give  10  log  Cg(l,u))=-27.5  dB  of  cancellation,  so 
the  choice  of  K=Kq  seems  justified.  However,  for  m=5,  K=Kg=2  gives 
10  log  C2(5»w)=16  dB  of  cancellation,  while  K=3  yields  19  dB  and  K=5, 
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22  dB  of  cancellation-  Of  course,  there  will  always  be  some  ambiguity  as 
to  how  much  more  rejection  must  be  provided  by  a  new  reference  to  justify 


I 

I 


its  use.  Certainly,  in  any  given  design  situation  can  be  used  as  a 
design  guide,  with  final  choices  made  by  consulting  the  curves  in 
Figures  K-6  through  K-62. 
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